Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,930,597 C→A 6.4% V391V (GTG→GTT)  apeB ← M18 family aminopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,930,5970CA 6.4% 195.9 / 2.8 63V391V (GTG→GTT) apeBM18 family aminopeptidase
Reads supporting (aligned to +/- strand):  ref base C (38/21);  new base A (2/2);  total (40/23)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

GAACCGCCGCGAACAGGCCGCGCCTGCTTACATCAAGGCAGCTCGCGGCTACGGTAGAACGCGGTCAGTACCTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGCACAAGTGGCGGAAGAACCCGGCCGTCTCGCTGTTGGTG  >  NC_002947/1930448‑1930741
                                                                                                                                                     |                                                                                                                                                
gAACCGCCGCGAACAGGCCGCGCCTGCTTACATCAAGGCAGCTCGCGGCTACGGTAGAACGCGGTCAGTACCTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATcc                                                                                                                                                  >  1:1007522/1‑150 (MQ=255)
   ccgccgCGAACAGGCCGCGCCTGCTTACATCAAGGCAGCTCGCGGCTACGGTAGAACGCGGTCAGTACCTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCAcc                                                                                                                                               >  1:154614/1‑150 (MQ=255)
       cgcgAACAGGCCGCGCCTGCTTACATCAAGGCAGCTCGCGGCTACGGTAGAACGCGGTCAGTACCTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTgc                                                                                                                                           >  2:226358/1‑150 (MQ=255)
                           gtACATCAAGGCAGCTCGCGGCTACGGTAGAACGCGGTCAGTACCTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCg                                                                                                                       <  1:430227/149‑1 (MQ=255)
                             aCATCAAGGCAGCTCGCGGCTACGGTAGAACGCGGTCAGTACCTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCGGGGGGCTGGCggg                                                                                                                     >  1:166138/1‑149 (MQ=255)
                               agcAAGGCAGCTCGCGGCTACGGTAGAACGCGGTCAGTACCTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGa                                                                                                                   >  2:863307/3‑150 (MQ=255)
                               aTCAAGGCAGCTCGCGGCTACGGTAGAACGCGGTCAGTACCTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGa                                                                                                                   >  2:1686225/1‑150 (MQ=255)
                               aTCAAGGCAGCTCGCGGCTACGGTAGAACGCGGTCAGTACCTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGa                                                                                                                   >  1:10632/1‑150 (MQ=255)
                                tCACGGCAGCACGAGGCTACGGTAGAACGCGGGCAGTACCTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGAt                                                                                                                  >  2:1553127/1‑150 (MQ=255)
                                tCAAGGCAGCTCGCGGCTACGGTAGAACGCGGTCAGTACCTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGAt                                                                                                                  >  2:1166161/1‑150 (MQ=255)
                                tCAAGGCAGCTCGCGGCTACGGTAGAACGCGGTCAGTACCTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGAt                                                                                                                  >  1:745023/1‑150 (MQ=255)
                                          tCGCGGCTACGGTAGAACGCGGTCAGTACCTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATg                                                                                                        >  1:1902383/1‑150 (MQ=255)
                                              ggCTACGGTAGAACGCGGTCAGTACCTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGGCCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGtgc                                                                                                    >  2:659988/1‑148 (MQ=255)
                                                 tACGGTAGAACGCGGTCAGTACCTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGc                                                                                                 <  1:1166161/150‑1 (MQ=255)
                                                         aaCGCGGTCAGTACCTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCAGCTTCTTGACACCCACCGTGCGCTCACCCAGGTCGCTGGCGAAGGTCTGCCCGATTGTCGAGCCACAGccc                                                                                         >  1:34265/1‑150 (MQ=255)
                                                          aCGCGGTCAGTACCTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCa                                                                                        <  1:1143868/150‑1 (MQ=255)
                                                          aCGCGGTCAGTACCTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCa                                                                                        <  1:1143873/150‑1 (MQ=255)
                                                                  aGTACCTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGAGGCTGGCGGTGATCGGCCCCCTGGTCGATCCACCGCCCCCGTCGCTg                                                                                >  2:771518/1‑150 (MQ=255)
                                                                    tACCTTCACCAGGTGCGCCAGGTCATGGCGGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTgcg                                                                              >  2:25407/1‑150 (MQ=255)
                                                                       cTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTgcgccc                                                                           >  1:344392/1‑148 (MQ=255)
                                                                       cTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGcac                                                                           >  2:239922/1‑150 (MQ=255)
                                                                           aCCAGGTGCGCCAGGTCATGGCTCCCACACAGCTCGCGAATGGAGGGCATGGCAAAGGTCGGCAAACCTATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACcaca                                                                       >  2:951222/1‑150 (MQ=255)
                                                                                gTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGct                                                                  >  2:12529/1‑150 (MQ=255)
                                                                                    gcCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGc                                                              >  1:1336355/1‑150 (MQ=255)
                                                                                     ccAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCa                                                             >  1:1778562/1‑150 (MQ=255)
                                                                                       aGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCAcc                                                           <  2:546253/150‑1 (MQ=255)
                                                                                       aGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTAGGCAAACCGATATCAACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATgg                                                                                                       <  1:776298/106‑1 (MQ=255)
                                                                                       aGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTAGGCAAACCGATATCAACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATgg                                                                                                       >  2:1748240/1‑106 (MQ=255)
                                                                                       aGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTAGGCAAACCGATATCAACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATgg                                                                                                       >  2:776298/1‑106 (MQ=255)
                                                                                       aGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTAGGCAAACCGATATCAACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATgg                                                                                                       <  1:1748240/106‑1 (MQ=255)
                                                                                         gTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCgg                                                         <  1:239922/150‑1 (MQ=255)
                                                                                          tCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGt                                                        >  1:1378046/1‑150 (MQ=255)
                                                                                          tCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGt                                                        >  2:1236907/1‑150 (MQ=255)
                                                                                            aTGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGctc                                                                 >  2:905655/1‑139 (MQ=255)
                                                                                            aTGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGctc                                                                 <  1:905655/139‑1 (MQ=255)
                                                                                                     acacaGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGt                                                                                                                     <  1:1632021/78‑1 (MQ=255)
                                                                                                     acacaGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGt                                                                                                                     >  2:1632021/1‑78 (MQ=255)
                                                                                                     acacaGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGc                                             >  1:1222510/1‑150 (MQ=255)
                                                                                                       acaGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCa                                           >  1:1157668/1‑150 (MQ=255)
                                                                                                        caGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACCCCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCAt                                          >  1:357587/1‑150 (MQ=255)
                                                                                                         aGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATg                                         <  2:10632/150‑1 (MQ=255)
                                                                                                            tCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGcac                                      >  2:434236/1‑150 (MQ=255)
                                                                                                              gcgAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGCACaa                                    >  2:957282/1‑150 (MQ=255)
                                                                                                               cgAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGCACAAg                                   <  2:34265/150‑1 (MQ=255)
                                                                                                                gAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTg                                                               >  2:1586639/1‑121 (MQ=255)
                                                                                                                gAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTg                                                               <  1:1586639/121‑1 (MQ=255)
                                                                                                                gAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGCACAAGt                                  <  2:1902383/150‑1 (MQ=255)
                                                                                                                gAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGCACAAGt                                  <  1:1686225/150‑1 (MQ=255)
                                                                                                                acctGGAGTGCATGGCAAAGGTCGGCGGATCAATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGCACAAGt                                  <  1:863307/147‑1 (MQ=255)
                                                                                                                    ggAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGCACAAGTGGCg                              <  1:12529/150‑1 (MQ=255)
                                                                                                                       gTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGCACAAGTGGCGgaa                           <  1:25407/150‑1 (MQ=255)
                                                                                                                       gTGCAAGGCAAGGGCCGTCAAGCCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGCAAGATCGTAAAGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGCACAAGTGGCGgaa                           <  1:1553127/150‑1 (MQ=255)
                                                                                                                           aTGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGCACAAGTGGCGGAAGAAc                       >  2:1728943/1‑150 (MQ=255)
                                                                                                                             ggCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGCACAAGTGGCGGAAGAAccc                     <  1:1728943/150‑1 (MQ=255)
                                                                                                                                aaaGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGCACAAGTGGCGGAAGAACCCGGc                  >  1:1043241/1‑150 (MQ=255)
                                                                                                                                aaaGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGCACAAGTGGCGGAAGAACCCGGc                  >  2:767885/1‑150 (MQ=255)
                                                                                                                                aaaGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGCACAAGTGGCGGAAGAACCCGGc                  >  1:332231/1‑150 (MQ=255)
                                                                                                                                        gCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGCACAAGt                                  >  1:1702859/1‑126 (MQ=255)
                                                                                                                                        gCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGCACAAGt                                  <  2:1702859/126‑1 (MQ=255)
                                                                                                                                        gCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGCACAAGTGGCGGAAGAACCCGGCCGTCTCGc          >  1:324366/1‑150 (MQ=255)
                                                                                                                                          aaaCCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAg                                                                                            <  2:650627/66‑1 (MQ=255)
                                                                                                                                          aaaCCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAg                                                                                            >  1:650627/1‑66 (MQ=255)
                                                                                                                                           aaCCGATATCCACCGTGCGCACACCCAGGTGGCTAGCGGTGACCCGCCCGCTGCTCCAGCCACAGCCCATGTCACTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGCACAAGTGGCGGAAGAACCCGGCCGTCTCGCTGt       >  2:946921/1‑150 (MQ=255)
                                                                                                                                                atatCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGCACAAGTGGCGGAAGAACCCGGCCGTCTCGCTGTtggtg  <  2:745023/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GAACCGCCGCGAACAGGCCGCGCCTGCTTACATCAAGGCAGCTCGCGGCTACGGTAGAACGCGGTCAGTACCTTCACCAGGTGCGCCAGGTCATGGCTGCCACACAGCTCGCGAATGGAGTGCATGGCAAAGGTCGGCAAACCGATATCCACCGTGCGCACACCCAGGTGGCTGGCGGTGATCGGCCCGATGGTCGAGCCACAGCCCATGTCGCTGCGCACCACAAAGCTCTGCACCGGTACTTCCTCGGCCATGCACAAGTGGCGGAAGAACCCGGCCGTCTCGCTGTTGGTG  >  NC_002947/1930448‑1930741

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