Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,560,480 T→C K109E (AAA→GAA)  PP_1366 ← transcriptional regulator MvaT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,560,4800TC99.4% 1139.1 / NA 323K109E (AAA→GAA) PP_1366transcriptional regulator MvaT
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (188/132);  minor base A (2/0);  total (191/132)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 325 total aligned reads displayed.

GAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAGTTTCGACTCTG  >  NC_002947/1560333‑1560626
                                                                                                                                                   |                                                                                                                                                  
gAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTGTCAACCACGTCGCCGCCCCAAGTGGCTTccc                                                                                                                                                  >  2:5655146/1‑150 (MQ=255)
  aCGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCAc                                                                                                                                                <  2:2544371/150‑1 (MQ=255)
   cGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCAct                                                                                                                                               <  1:527774/150‑1 (MQ=255)
     cTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCActct                                                                                                                                             >  2:1530361/1‑150 (MQ=255)
      tATTCAATCGGCTGATGGCAGATAAACAATAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTACAGGAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCtt                                                                                                                                            >  1:5967486/1‑150 (MQ=255)
      tATTCAATCGGCTGATGGCAGATAAACAAAATCGCCGGCAATGTGCCGGCGTTTTTTTTGTGCGAAGGAAGAGGGCGTATGGGGTTAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCtt                                                                                                                                            >  1:1986704/1‑150 (MQ=255)
           aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCACTGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCACGTGGCTTCCCACTCTTTCAgg                                                                                                                                       >  2:5195430/1‑150 (MQ=255)
           aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAgg                                                                                                                                       >  1:4475347/1‑150 (MQ=255)
           aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTAAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAgg                                                                                                                                       >  1:4522079/1‑150 (MQ=255)
               ggCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCt                                                                                                                                   >  2:3852594/1‑150 (MQ=255)
                  tGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTtgt                                                                                                                                >  1:888676/1‑150 (MQ=255)
                   gATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTtgtg                                                                                                                               >  2:2217319/1‑150 (MQ=255)
                     tGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGt                                                                                                                             <  1:143135/150‑1 (MQ=255)
                       gCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCGTGTGGTTg                                                                                                                           >  2:3111622/1‑150 (MQ=255)
                         aGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCa                                                                                                                                                 <  1:4063218/126‑1 (MQ=255)
                             aaaCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCt                                                                                                                     >  1:4933513/1‑150 (MQ=255)
                              aaCAAAAACGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCtt                                                                                                                    >  2:1635721/1‑150 (MQ=255)
                               aCAAAATCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCttt                                                                                                                   >  2:1830709/1‑150 (MQ=255)
                               aCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCttt                                                                                                                   <  1:420919/150‑1 (MQ=255)
                                cAAAAACGCCGGCACTGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTg                                                                                                                  <  2:3904461/150‑1 (MQ=255)
                                cAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTg                                                                                                                  <  2:4989025/150‑1 (MQ=255)
                                  aaaaCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGt                                                                                                                >  2:480208/1‑150 (MQ=255)
                                  aaaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGt                                                                                                                <  1:3639846/150‑1 (MQ=255)
                                  aaaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGt                                                                                                                >  1:1892414/1‑150 (MQ=255)
                                    aaCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGttt                                                                                                              >  1:1260925/1‑150 (MQ=255)
                                    aaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGttt                                                                                                              <  1:3503799/150‑1 (MQ=255)
                                    aaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGttt                                                                                                              <  2:4731008/150‑1 (MQ=255)
                                    aaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGttt                                                                                                              <  1:2574265/150‑1 (MQ=255)
                                    aaCGCCGGCAATGTGCCGGCGTTTTTGTGGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCCGCTTTCAACCACGTCGCCGCCCCACGTGGCTTCCCACTCGGTCAGGGGGTTGTGGTTGCCGCCTTTGGgtt                                                                                                              >  2:5065402/1‑150 (MQ=255)
                                     aCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGtt                                                                                                                            <  1:693883/135‑1 (MQ=255)
                                     aCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGtt                                                                                                                            >  2:693883/1‑135 (MQ=255)
                                     aCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTc                                                                                                             >  2:5873198/1‑150 (MQ=255)
                                     aCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTc                                                                                                             >  2:5451326/1‑150 (MQ=255)
                                           gCAATGTGCCGGCGTTTTTGTTGTGCGCAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCAc                                                                                                       <  2:3713664/150‑1 (MQ=255)
                                                    cGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGtt                                                                                              >  2:4866007/1‑150 (MQ=255)
                                                        gTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGc                                                                                          <  2:265964/150‑1 (MQ=255)
                                                              gttgtGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTc                                                                                    >  2:4546451/1‑150 (MQ=255)
                                                              gttgtGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTc                                                                                    <  2:4264615/150‑1 (MQ=255)
                                                                  tgCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGt                                                                                >  1:1221886/1‑150 (MQ=255)
                                                                  tgCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTGTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGt                                                                                >  2:4080680/1‑150 (MQ=255)
                                                                   gCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTa                                                                               >  1:3187965/1‑150 (MQ=255)
                                                                     gAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCActc                                                                                                                                              >  1:75957/1‑85 (MQ=255)
                                                                     gAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTAtt                                                                             >  1:4161850/1‑150 (MQ=255)
                                                                     gAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTAtt                                                                             <  1:5462736/150‑1 (MQ=255)
                                                                     gAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTAtt                                                                             <  1:369235/150‑1 (MQ=255)
                                                                     gAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTAtt                                                                             >  1:3384879/1‑150 (MQ=255)
                                                                         cAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCtt                                                                         <  1:1588381/150‑1 (MQ=255)
                                                                           agagGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGa                                                                       >  2:2903667/1‑150 (MQ=255)
                                                                           agagGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGa                                                                       <  1:3778692/150‑1 (MQ=255)
                                                                                 cGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTtgc                                                                 <  2:2814340/150‑1 (MQ=255)
                                                                                   tATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTtgcgt                                                               <  1:243579/150‑1 (MQ=255)
                                                                                      gCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTtgcgtgcg                                                            >  1:801954/1‑150 (MQ=255)
                                                                                       cGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTgcgc                                                           <  2:2282219/150‑1 (MQ=255)
                                                                                        ggATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTgcgcg                                                          >  1:5860331/1‑150 (MQ=255)
                                                                                            aaTCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTg                                                      <  1:3418304/150‑1 (MQ=255)
                                                                                             aTCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTgg                                                     <  1:4597589/150‑1 (MQ=255)
                                                                                              tCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGt                                                    >  1:2908858/1‑150 (MQ=255)
                                                                                                aGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGa                                                  >  1:2651851/1‑150 (MQ=255)
                                                                                                    cagcagGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGt                                                                                      <  2:639179/110‑1 (MQ=255)
                                                                                                    cagcagGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGt                                              >  1:480557/1‑150 (MQ=255)
                                                                                                     agcagGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTa                                             >  2:3777505/1‑150 (MQ=255)
                                                                                                     agcagGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTa                                             <  1:3890561/150‑1 (MQ=255)
                                                                                                     agcagGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTa                                             <  2:1986704/150‑1 (MQ=255)
                                                                                                        agGGTTGCCCAGCTTTCAACCACGTCGCCGACCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTc                                          >  2:3520353/1‑150 (MQ=255)
                                                                                                         gggTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCt                                         <  2:4463686/150‑1 (MQ=255)
                                                                                                          ggTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCtt                                        >  2:4606216/1‑150 (MQ=255)
                                                                                                            ttGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTTCGCGCTTGGTGACGGTAGTCTTGa                                      <  1:5551215/150‑1 (MQ=255)
                                                                                                            ttGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGa                                      <  2:888805/150‑1 (MQ=255)
                                                                                                            ttGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGa                                      <  1:5379739/150‑1 (MQ=255)
                                                                                                               cccAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCg                                   >  2:5794928/1‑150 (MQ=255)
                                                                                                                ccAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGa                                  >  1:3236928/1‑150 (MQ=255)
                                                                                                                ccAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGa                                  >  1:4263950/1‑150 (MQ=255)
                                                                                                                 cAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCtt                                        >  1:2736752/1‑143 (MQ=255)
                                                                                                                 cAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGaa                                 >  1:2241037/1‑150 (MQ=255)
                                                                                                                 cAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGaa                                 >  1:4341088/1‑150 (MQ=255)
                                                                                                                  aGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAAc                                >  1:1169808/1‑150 (MQ=255)
                                                                                                                     tttCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCAc                             <  2:801954/150‑1 (MQ=255)
                                                                                                                     tttCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCAc                             >  2:347103/1‑150 (MQ=255)
                                                                                                                     tttCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCAc                             >  1:3058517/1‑150 (MQ=255)
                                                                                                                     tttCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCAc                             >  1:4859957/1‑150 (MQ=255)
                                                                                                                         aaCCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAgg                         >  1:3105860/1‑150 (MQ=255)
                                                                                                                            cACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGc                      <  1:2596452/150‑1 (MQ=255)
                                                                                                                             aCGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTgcgcgc                                                         >  2:3981427/1‑114 (MQ=255)
                                                                                                                             aCGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGcc                     >  1:1024291/1‑150 (MQ=255)
                                                                                                                             aCGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTAGTGCTGGACCGTGCGTGCGCGCTGGGTGACGGTAGTCTTGACCGAACCACGAGGTGcc                     >  1:3790900/1‑150 (MQ=255)
                                                                                                                                gcgccgcCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTg                  <  2:1649525/149‑1 (MQ=255)
                                                                                                                                  gccgccCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCt                <  2:2389915/150‑1 (MQ=255)
                                                                                                                                      ccccAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCtt                                                       >  2:1858210/1‑107 (MQ=255)
                                                                                                                                      ccccAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAGt            >  1:5049902/1‑150 (MQ=255)
                                                                                                                                      ccccAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAGt            <  2:3656196/150‑1 (MQ=255)
                                                                                                                                      ccccAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAGt            >  2:3228946/1‑150 (MQ=255)
                                                                                                                                      acccAAGTGGCTTCCCACTATTTCAGTGTCTTGAGGTTGCCGCCTTTGGTTTCGATCACTTCACAGTTGTGCTTGTTCTTGTATTGCTTGACATTTCGTGCGCTCTTGGTGACGGTATTCTTGACCGAACCACGAGGTGCCTTGCTCAGt            >  2:12679/2‑150 (MQ=255)
                                                                                                                                        ccATGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTGCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTTACGTTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAGttt          <  2:5298668/150‑1 (MQ=255)
                                                                                                                                        ccAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAGttt          >  1:4374448/1‑150 (MQ=255)
                                                                                                                                        ccAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAGttt          <  1:1830709/150‑1 (MQ=255)
                                                                                                                                        ccAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAGttt          <  1:1530361/150‑1 (MQ=255)
                                                                                                                                         cAAGTGGCTTCCCACTCATTCAGTGTCTTGTGGTTGCAGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGAACTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAGTTTc         <  1:5148999/150‑1 (MQ=255)
                                                                                                                                          aaGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAGTTTCg        >  2:1438590/1‑150 (MQ=255)
                                                                                                                                          aaGTGGCTTCCCACTCTTTCAGGGTCTAGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTGGTATTGCTTGGCCTTGCGTGCGCGCTTGGTGACGGTAGTATTGACCGAACCACGAGGTGCCTTGCGCAGGTTCg        >  2:1124213/1‑150 (MQ=255)
                                                                                                                                           aGTGGCTTCCCACTCTTTCAGGGTCTTGTGATTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAGTTTCGa       <  1:4014398/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAGTTTCGACTCTG  >  NC_002947/1560333‑1560626

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: