Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,388,285 G→A 7.3% A281V (GCT→GTT) ‡ aaeA ← p‑hydroxybenzoic acid efflux system component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,388,2850GA 7.3% 257.5 / 7.6 83A281V (GCT→GTT) ‡aaeAp‑hydroxybenzoic acid efflux system component
Reads supporting (aligned to +/- strand):  ref base G (35/41);  new base A (2/4);  total (37/46)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

CCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTGGTGACCCC  >  NC_000913/3388137‑3388433
                                                                                                                                                     |                                                                                                                                                    
ccAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAg                                                                                                                                                      >  3:101130/1‑149 (MQ=255)
ccAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAg                                                                                                                                                      >  2:136196/1‑149 (MQ=255)
 cAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGc                                                                                                                                                     >  1:385722/1‑149 (MQ=255)
  aGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCa                                                                                                                                                    <  1:673482/149‑1 (MQ=255)
     tttACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAgtg                                                                                                                                                 <  2:668865/149‑1 (MQ=255)
     tttACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAgtg                                                                                                                                                 <  4:2391/149‑1 (MQ=255)
        aCCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGCTA‑CAgtggtg                                                                                                                                              <  3:314176/149‑1 (MQ=255)
        aCCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑CTACAgtggtg                                                                                                                                              <  2:536167/149‑1 (MQ=255)
         ccGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑CTACAgtggtgg                                                                                                                                             <  3:505239/149‑2 (MQ=255)
         ccGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑CTACAgtggtgg                                                                                                                                             <  3:460338/149‑2 (MQ=255)
               aGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAg                                                                                                                                       <  4:163422/149‑1 (MQ=255)
                  gAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGcc                                                                                                                                    >  3:198306/1‑149 (MQ=255)
                  gAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGcc                                                                                                                                    >  1:120423/1‑149 (MQ=255)
                   aaTATGTTGGTTAGCAATCGAGAAAATACCCAGCGTGATTAACCAAAGTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAAACCTGGCTTTCGTCGCGGTCTTGTTTGCCAGTGACCACGGCTG‑TAGGAGTGGTGCCCGCAGGCCa                                                                                                                                   <  4:629944/149‑1 (MQ=255)
                     tatGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGa                                                                                                                                 <  4:198306/149‑1 (MQ=255)
                       tGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATg                                                                                                                               >  4:263087/1‑149 (MQ=255)
                               aGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTg                                                                                                                       <  2:345384/149‑1 (MQ=255)
                                    tGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGt                                                                                                                  >  3:131314/1‑149 (MQ=255)
                                     ggAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGtt                                                                                                                 <  1:823802/149‑1 (MQ=255)
                                     ggAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGtt                                                                                                                 <  2:435154/149‑1 (MQ=255)
                                      gagaAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTg                                                                                                                >  2:62852/1‑149 (MQ=255)
                                        gaAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTc                                                                                                              <  3:231019/149‑1 (MQ=255)
                                          aaaTACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCga                                                                                                            <  3:630474/149‑1 (MQ=255)
                                             tACCCACCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGAc                                                                                                         <  2:390602/149‑1 (MQ=255)
                                              aCCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACg                                                                                                        >  4:328826/1‑149 (MQ=255)
                                                ccATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGaa                                                                                                      <  2:747077/149‑1 (MQ=255)
                                                ccATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGaa                                                                                                      >  1:650375/1‑149 (MQ=255)
                                                 cATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTCCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAAt                                                                                                     >  1:160553/1‑149 (MQ=255)
                                                   tCGTGCTTAACCAAACTCACGCAGGCGATGAGCCATGTTACGGAAGCACGAATCCTTCCTTTCGCCGCGATCTTGTTTGTCAGTGATAACCACTG‑TAGCAGGGGCGTTCGCAAGCCAGATGTTCTCTTGCTCGTTGTCGGGACGAATAc                                                                                                   >  4:837128/1‑149 (MQ=255)
                                                   tCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGTTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAc                                                                                                   <  4:694844/149‑1 (MQ=255)
                                                            aCCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGaa                                                                                          >  2:141970/1‑149 (MQ=255)
                                                             ccAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAac                                                                                         <  2:264093/149‑1 (MQ=255)
                                                             ccAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAac                                                                                         >  4:200050/1‑149 (MQ=255)
                                                                    cACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTgc                                                                                  <  4:77408/149‑1 (MQ=255)
                                                                        cAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCaa                                                                              <  4:640389/149‑1 (MQ=255)
                                                                         aGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAAc                                                                                               <  1:285746/131‑1 (MQ=255)
                                                                         aGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAg                                                                             >  1:428908/1‑149 (MQ=255)
                                                                         aGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTAGTTGTCGAGACGAATACGAAc                                                                                               >  2:285746/1‑131 (MQ=255)
                                                                          ggCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGa                                                                            >  1:242529/1‑149 (MQ=255)
                                                                             gATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGc                                                                         >  1:565070/1‑149 (MQ=255)
                                                                             gATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGc                                                                         >  4:445919/1‑149 (MQ=255)
                                                                                  gCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGa                                                                                                                                 >  3:222844/1‑88 (MQ=255)
                                                                                  gCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGa                                                                                                                                 <  4:222844/88‑1 (MQ=255)
                                                                                      ttttACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCa                                                                <  3:285837/149‑1 (MQ=255)
                                                                                      ttttACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCa                                                                >  1:512692/1‑149 (MQ=255)
                                                                                       tttACGGAAGCACGAATCCTGGCTTTTGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCaa                                                               <  1:641894/149‑1 (MQ=255)
                                                                                        ttACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGTCGCACCCATTCAAg                                                              <  1:671013/149‑1 (MQ=255)
                                                                                        ttACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAg                                                              <  3:121885/149‑1 (MQ=255)
                                                                                        ttACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAg                                                              <  2:120423/149‑1 (MQ=255)
                                                                                            ggAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTa                                                          >  3:86057/1‑149 (MQ=255)
                                                                                            ggAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTa                                                          >  3:477172/1‑149 (MQ=255)
                                                                                              acgaaCGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTaga                                                        <  1:436350/147‑1 (MQ=255)
                                                                                              aagaaCGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTaga                                                        <  4:532248/149‑1 (MQ=255)
                                                                                              aagaaCGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTaga                                                        <  2:497042/149‑1 (MQ=255)
                                                                                                gaaCGAATCCTTGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGt                                                      >  3:704957/1‑149 (MQ=255)
                                                                                                gaaCGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGt                                                      >  1:314412/1‑149 (MQ=255)
                                                                                                     aaTCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCtata                                                 <  3:660407/149‑1 (MQ=255)
                                                                                                     aaTCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCtata                                                 <  2:650375/149‑1 (MQ=255)
                                                                                                       tCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGaa                                                                                                      >  2:691267/1‑94 (MQ=255)
                                                                                                       tCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGaa                                                                                                      <  1:691267/94‑1 (MQ=255)
                                                                                                       tCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGt                                               >  4:661812/1‑149 (MQ=255)
                                                                                                        ccTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTc                                              <  3:661812/149‑1 (MQ=255)
                                                                                                        ccTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTc                                              <  3:263087/149‑1 (MQ=255)
                                                                                                              tttCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTtc                                                                                                                            >  4:158055/1‑65 (MQ=255)
                                                                                                              tttCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTtc                                                                                                                            <  3:158055/65‑1 (MQ=255)
                                                                                                              tttCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATc                                        >  4:510213/1‑149 (MQ=255)
                                                                                                              tttCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATc                                        >  1:843602/1‑149 (MQ=255)
                                                                                                                  gtcgCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCtt                                    <  3:445919/149‑1 (MQ=255)
                                                                                                                  gtcgCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCtt                                    <  4:337993/149‑1 (MQ=255)
                                                                                                                   tcgCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCttt                                   <  4:86057/149‑1 (MQ=255)
                                                                                                                   tcgCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCttt                                   <  4:392484/149‑1 (MQ=255)
                                                                                                                   tcgCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCttt                                   <  1:578310/149‑1 (MQ=255)
                                                                                                                      cgATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTgtc                                >  1:271110/1‑149 (MQ=255)
                                                                                                                      cgATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTgtc                                <  4:84543/149‑1 (MQ=255)
                                                                                                                       gATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTgtcg                               >  4:531804/1‑149 (MQ=255)
                                                                                                                        aTCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTgtcgt                              >  4:691566/1‑149 (MQ=255)
                                                                                                                         tCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTgtcgtc                             <  3:604317/149‑1 (MQ=255)
                                                                                                                                 gCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGtgc                     >  1:823709/1‑149 (MQ=255)
                                                                                                                                        ggcaccacTG‑CTACAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGc              >  1:862796/3‑149 (MQ=255)
                                                                                                                                        ggcaccacTG‑CTACAGTGGTGCCCGCAGGCCAGAGGTTCTCTTGCTGGTTGGCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGc              >  2:703992/3‑149 (MQ=255)
                                                                                                                                          caccacTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGt            >  4:398185/1‑149 (MQ=255)
                                                                                                                                           accacTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGtt           >  3:59428/1‑149 (MQ=255)
                                                                                                                                           accacTG‑CTACAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGtt           <  2:862796/149‑1 (MQ=255)
                                                                                                                                            ccacTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTg          <  4:450118/149‑1 (MQ=255)
                                                                                                                                             cacTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTgg         >  3:210439/1‑149 (MQ=255)
                                                                                                                                                     gCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTGGTGAcccc  >  2:106085/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG‑TAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTGGTGACCCC  >  NC_000913/3388137‑3388433

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑