Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,972,881 A→T 5.0% F223I (TTT→ATT)  flgL ← flagellar hook‑associated protein FlgL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,972,8810AT 5.0% 247.2 / 2.0 79F223I (TTT→ATT) flgLflagellar hook‑associated protein FlgL
Reads supporting (aligned to +/- strand):  ref base A (34/40);  new base T (2/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

GCGAGATCATCTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGGTGTCGTTGGTCGC  >  NC_002947/4972737‑4973028
                                                                                                                                                |                                                                                                                                                    
gCGAGATCATCTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGa                                                                                                                                                  <  3:117086/149‑1 (MQ=255)
    gaTCATCTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGc                                                                                                                                              <  6:320680/149‑1 (MQ=255)
     atcatcTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGcc                                                                                                                                             >  8:192321/1‑149 (MQ=255)
     atcatcTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGAACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGcc                                                                                                                                             <  4:341849/149‑1 (MQ=255)
        atcTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAATCAGTTTGCCgct                                                                                                                                          <  8:15862/149‑3 (MQ=255)
        atcTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCg                                                                                                                                            >  7:277077/1‑147 (MQ=255)
          cTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGc                                                                                                                                              <  6:251263/143‑1 (MQ=255)
          cTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGc                                                                                                                                              >  5:251263/1‑143 (MQ=255)
          cTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAgg                                                                                                                                        >  4:320921/1‑149 (MQ=255)
           ttGCCGCCTTCGGTATTGACGTCGAACTAGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAggt                                                                                                                                       >  3:314076/1‑149 (MQ=255)
             gccgccTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAggtg                                                                                                                                      <  3:340963/148‑1 (MQ=255)
             gccgccTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAg                                                                                                                                         <  7:192321/145‑1 (MQ=255)
              ccgacTTCGGTATTCACGTCGAACTTGCCCTTTCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGAACTGCTGAAGCTGACAGAATACTGCTGACCACCCGCAAAATCATTGCCGAAggtggt                                                                                                                                    >  4:206855/1‑149 (MQ=255)
               cgccTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGGTGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGgtggtg                                                                                                                                   >  4:101970/1‑149 (MQ=255)
                gccTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGt                                                                                                                                           <  8:123048/140‑1 (MQ=255)
                gccTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGt                                                                                                                                           >  7:123048/1‑140 (MQ=255)
                  cTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGcc                                                                                                                                >  3:133263/1‑149 (MQ=255)
                     cGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCg                                                                                                                               <  1:104058/147‑1 (MQ=255)
                      ggTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAg                                                                                                                             <  4:133263/148‑1 (MQ=255)
                       gTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGGTGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGt                                                                                                                           <  3:101970/149‑1 (MQ=255)
                       gTATTGACGTCGAACTTGCCCTTGCCGGAGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAgg                                                                                                                            >  3:172864/1‑148 (MQ=255)
                          ttGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCaa                                                                                                                        <  4:172864/149‑1 (MQ=255)
                               gTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGgc                                                                                                                  >  7:277856/1‑149 (MQ=255)
                                  gAACTTGGCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCTGCGGCATCGCTCTCGCTGTAGTCTGTGGCTCTGGTGAACGAGACAGTATACGGCTGACCAGCGGCAAACTGATTGCCGTAg                                                                                                                                         <  5:199316/124‑1 (MQ=255)
                                  gAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGACGGCGTCATCCCTAACGCTGTACTCGTTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCCGCAAACTGATTGCCGTAg                                                                                                                                         >  6:199316/1‑124 (MQ=255)
                                   aaCTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGgccgc                                                                                                               >  3:3439/1‑148 (MQ=255)
                                   aaCTTGCCCTTGCCGGCGGTTTCATTGGTGACGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGgccgcc                                                                                                              <  4:123837/149‑1 (MQ=255)
                                     ctagcccttgccGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAg                                                                                                            <  4:314076/146‑1 (MQ=255)
                                      ttgcccttgccGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGa                                                                                                           >  1:266332/1‑149 (MQ=255)
                                         cccttgccGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAg                                                                                                        <  3:206855/149‑1 (MQ=255)
                                                  gcggTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCtt                                                                                               <  7:314119/149‑1 (MQ=255)
                                                   cggTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTg                                                                                              <  5:345469/149‑1 (MQ=255)
                                                     gTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAg                                                                                                        <  8:357287/137‑1 (MQ=255)
                                                     gTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAg                                                                                                        >  7:357287/1‑137 (MQ=255)
                                                     gTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGcc                                                                                            <  8:277077/149‑1 (MQ=255)
                                                     gTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGcc                                                                                            <  7:352573/149‑1 (MQ=255)
                                                         cATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTc                                                                                         <  2:266332/148‑1 (MQ=255)
                                                          aTTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGACGTCAt                                                                                       >  3:248228/1‑149 (MQ=255)
                                                               tGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGc                                                                                  <  4:3439/149‑1 (MQ=255)
                                                               tGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACATTGCCGTAATCCGGc                                                                                  >  4:251821/1‑149 (MQ=255)
                                                                gATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGc                                                                                  >  6:293126/1‑148 (MQ=255)
                                                                    tCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACgg                                                                             <  3:107483/149‑1 (MQ=255)
                                                                    tCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACgg                                                                             >  7:220361/1‑149 (MQ=255)
                                                                     cATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGt                                                                            >  5:288387/1‑149 (MQ=255)
                                                                      aTTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTg                                                                           >  5:140753/1‑149 (MQ=255)
                                                                          ccGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTAc                                                                       <  3:320921/149‑1 (MQ=255)
                                                                                 cATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACATGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTc                                                                <  4:389274/149‑1 (MQ=255)
                                                                                  aTCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACgg                                                                             >  2:207124/1‑135 (MQ=255)
                                                                                  aTCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACgg                                                                             <  1:207124/135‑1 (MQ=255)
                                                                                     ccTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACt                                                            >  7:122859/1‑149 (MQ=255)
                                                                                     ccTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACt                                                            >  1:147319/1‑149 (MQ=255)
                                                                                     ccTCACGCTGTACTCGGTGGCACTGGTGAAGGTAACAGAATACGGCTGACCAGCGGCAACCTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACt                                                            >  1:378302/1‑149 (MQ=255)
                                                                                       tCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACt                                                            <  1:123291/147‑1 (MQ=255)
                                                                                       tCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACt                                                            >  2:123291/1‑147 (MQ=255)
                                                                                         aCGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGt                                                        <  1:83644/149‑1 (MQ=255)
                                                                                         aCGCTGTACTCGGTGGCACTGGGGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTggg                                                         <  6:140753/148‑1 (MQ=255)
                                                                                                  tCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCtgtt                                               >  3:34273/1‑149 (MQ=255)
                                                                                                        gCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTtggt                                         <  8:220361/149‑1 (MQ=255)
                                                                                                          aCTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCAt                                                                                       <  5:138978/101‑1 (MQ=255)
                                                                                                          aCTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCAt                                                                                       >  6:138978/1‑101 (MQ=255)
                                                                                                          aCTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTtggtgg                                       >  6:344593/1‑149 (MQ=255)
                                                                                                              gTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGcc                                     <  2:147319/147‑1 (MQ=255)
                                                                                                              gTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTg                                   >  5:256876/1‑149 (MQ=255)
                                                                                                                    gTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGa                              >  2:283256/1‑148 (MQ=255)
                                                                                                                       aCAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGaaaa                           <  1:283256/148‑1 (MQ=255)
                                                                                                                        cAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATg                         <  6:256876/149‑1 (MQ=255)
                                                                                                                         aGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCgg                                                    <  5:68223/121‑1 (MQ=255)
                                                                                                                         aGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCgg                                                    >  6:68223/1‑121 (MQ=255)
                                                                                                                          gAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATg                         <  3:219039/147‑1 (MQ=255)
                                                                                                                          gAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATg                         >  4:219039/1‑147 (MQ=255)
                                                                                                                          gAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCt                       <  6:288387/149‑1 (MQ=255)
                                                                                                                          gAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCt                       <  5:344593/149‑1 (MQ=255)
                                                                                                                           aaTACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTg                      >  2:39725/1‑149 (MQ=255)
                                                                                                                                 gCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCg                <  8:277856/149‑1 (MQ=255)
                                                                                                                                   ccacctacGGCAATCAGTTTGCCGTAGGTGGTGGCCGAGGTCAACAGGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCgg               >  1:89526/8‑149 (MQ=255)
                                                                                                                                      ccAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGAc                                                                               <  2:353690/81‑1 (MQ=255)
                                                                                                                                      ccAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGAc                                                                               >  1:353690/1‑81 (MQ=255)
                                                                                                                                       ctacGGCAATCAGTTTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCt                                                                         <  4:338605/83‑1 (MQ=255)
                                                                                                                                       ctacGGCAATCAGTTTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCt                                                                         >  3:338605/4‑86 (MQ=255)
                                                                                                                                               aaCTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGGTGTCGTTGGTCGc  >  8:377397/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GCGAGATCATCTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGGTGTCGTTGGTCGC  >  NC_002947/4972737‑4973028

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: