Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,656,533 T→C 5.1% Q93R (CAA→CGA)  PP_1449 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,656,5330TC 5.1% 357.7 / 4.8 115Q93R (CAA→CGA) PP_1449hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (56/53);  new base C (3/3);  total (59/56)
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.12e-01

Only 100 of 120 total aligned reads displayed.

TGGGGTTGGCCAGTATGTAGTCGGCCGGGCGCCCGAAGATTTCCTGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGATGACCGGCACCCCGTGGCCGTTATGGAT  >  NC_002947/1656384‑1656672
                                                                                                                                                     |                                                                                                                                           
tGGGGTTGGCCAGTATGTAGTCGGCCGGGCGCCCGAAGATTTCCTGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAAGGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGTGCGCCGGCCAGtt                                                                                                                                             <  1:1498745/150‑1 (MQ=255)
  gggTTGGCCAGTATGTAGTCGGCCGGGCGCCCGAAGATTTCCTGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGGGCTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGc                                                                                                                                           >  1:1455626/1‑150 (MQ=255)
  gggTTGGCCAGTATGTAGTCGGCCGGGCGCCCGAAGATTTCCTGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGc                                                                                                                                           <  1:778061/150‑1 (MQ=255)
     ttGGCCAGTATGTAGTCGGCCGGGCGCCCGAAGATTTCCTGCGGCCCCTCGATCAGCGGGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGAc                                                                                                                                        >  1:321150/1‑150 (MQ=255)
         ccAGTATGTAGTCGGCCGGGCGCCCGAAGATTTCCTGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGcc                                                                                                                                    >  1:2338562/1‑150 (MQ=255)
          cAGTATGTAGTCGGCCGGGCGCCCGAAGATTTCCTGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGccc                                                                                                                                   >  1:1828425/1‑150 (MQ=255)
          cAGTATGTAGTCGGCCGGGCGCCCGAAGATTTCCTGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTCGCCCTGGAACTGCGGGTTAGCCGCCAGGGCGCCGGCCAGTTGCGACTGccc                                                                                                                                   >  2:260013/1‑150 (MQ=255)
            gTATGTAGTCGGCCGGGCGCCCGAAGATTTCCTGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCgg                                                                                                                                 <  1:1903306/150‑1 (MQ=255)
              aTGTAGTCGGCCGGGCGCCCGAAGATTCCCTGCGGCCCCTTGATCAGCGAGGCGTTGCGGCTGATCACTTAGTTAAGAATGGTCGAGGCGCCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCt                                                                                                                               <  1:1916499/150‑1 (MQ=255)
                   gTCGGCCGGGCGCCCGAAGATTTCCTGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGt                                                                                                                          >  1:1848972/1‑150 (MQ=255)
                     cGGCCGGGCGCCCGAAGATTTCCTGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTTGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTgg                                                                                                                        <  2:1148363/150‑1 (MQ=255)
                        ccGGGCGCCCGAAGATTTCCTGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGGCAGTCGCAACTGGCCGGc                                                                                                                                <  1:1518247/139‑1 (MQ=255)
                        ccGGGCGCCCGAAGATTTCCTGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGGCAGTCGCAACTGGCCGGc                                                                                                                                >  2:1518247/1‑139 (MQ=255)
                         cGGGCGCCCGAAGATTTCCTGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCgtt                                                                                                                    <  1:1017327/150‑1 (MQ=255)
                          gggCGCCCGAAGATTTCCTGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCgttg                                                                                                                   >  2:2071772/1‑150 (MQ=255)
                                ccGAAGATTTCCTGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGGCAGTCGCAACTGGCCGGc                                                                                                                                <  1:1046595/131‑1 (MQ=255)
                                    aGATTTCCTGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCa                                                                                                         <  2:1617293/150‑1 (MQ=255)
                                     gATTTCCTGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAg                                                                                                        <  2:1475900/150‑1 (MQ=255)
                                     gATTTCCTGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCAGCTGCGGTGGCGTTGTTCAGCACCAg                                                                                                        <  2:664943/150‑1 (MQ=255)
                                        ttCCTGCGGCCCCTCGATCAGCGAGGCGTTCCGGCACATCCCATGGTCTGTACAGGTCGAGGCGGCCTGGCCTTTGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGcc                                                                                                     >  1:1023605/1‑150 (MQ=255)
                                           cTGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCgg                                                                                                  >  1:1557543/1‑150 (MQ=255)
                                            tGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGcc                                                                                                                                    <  2:995415/115‑1 (MQ=255)
                                            tGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGcc                                                                                                                                    >  1:995415/1‑115 (MQ=255)
                                            tGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCggg                                                                                                 <  2:2338562/150‑1 (MQ=255)
                                                 ccccTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGc                                                                                            >  1:1557713/1‑150 (MQ=255)
                                                  cccTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGCTCTAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCa                                                                                           >  1:1757493/1‑150 (MQ=255)
                                                   ccTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCGCGGCTGCGGGGGCGGTGTTCAGCCCCAGGCCCGGGGTGGCaa                                                                                          >  2:351083/1‑150 (MQ=255)
                                                    cTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAAc                                                                                         >  1:965240/1‑150 (MQ=255)
                                                    cTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAAc                                                                                         >  1:2092507/1‑150 (MQ=255)
                                                          cAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCACCGTTGTa                                                                                   <  2:2402344/150‑1 (MQ=255)
                                                          cAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTa                                                                                   <  2:321150/150‑1 (MQ=255)
                                                          cAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTa                                                                                   <  2:2092507/150‑1 (MQ=255)
                                                                ggCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGAt                                                                             <  2:715057/150‑1 (MQ=255)
                                                                ggCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGAt                                                                             <  1:378834/150‑1 (MQ=255)
                                                                       cGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTgg                                                                      <  2:1557713/150‑1 (MQ=255)
                                                                       cGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTgg                                                                      <  1:167848/150‑1 (MQ=255)
                                                                           tGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGtt                                                                                                                                             >  2:1578846/1‑75 (MQ=255)
                                                                           tGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGtt                                                                                                                                             <  1:1578846/75‑1 (MQ=255)
                                                                           tGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTtgt                                                                  <  1:1097012/150‑1 (MQ=255)
                                                                             aTCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTtgtgt                                                                >  1:511916/1‑150 (MQ=255)
                                                                              tCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTtgtgtg                                                               >  1:175524/1‑150 (MQ=255)
                                                                               cACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGa                                                              >  2:136399/1‑150 (MQ=255)
                                                                                aCTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGAc                                                             >  2:1971093/1‑150 (MQ=255)
                                                                                 cTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGt                                                                                               <  2:1737548/115‑1 (MQ=255)
                                                                                  ttCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGc                                                                                                                             >  1:1203399/1‑84 (MQ=255)
                                                                                  ttCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGc                                                                                                                             <  2:1203399/84‑1 (MQ=255)
                                                                                  ttCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTtg                                                                   <  2:1106092/142‑1 (MQ=255)
                                                                                    cGTTAAGAATGGTCGAGGCGGCCTGGCACTGGAAATGCGGGTTGGCGGCCCTTGCTCCGGCCAGTTGCGACTTCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGcc                                                         >  1:1023479/1‑150 (MQ=255)
                                                                                     gTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGc                                                                                                                                     <  1:789886/73‑1 (MQ=255)
                                                                                     gTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGc                                                                                                                                     >  2:789886/1‑73 (MQ=255)
                                                                                     gTTAACAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGCCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGccc                                                        <  2:1023479/150‑1 (MQ=255)
                                                                                      ttAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGt                                                                                 <  2:1873667/124‑1 (MQ=255)
                                                                                         aGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTg                                                    <  1:118761/150‑1 (MQ=255)
                                                                                          gAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTtgtgt                                                                >  2:2176800/1‑137 (MQ=255)
                                                                                           aaTGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGa                                                                              >  1:1280440/1‑122 (MQ=255)
                                                                                           aaTGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGa                                                                              <  2:1280440/122‑1 (MQ=255)
                                                                                                tCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCggg                                                                                                 >  1:2137785/1‑98 (MQ=255)
                                                                                                tCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCggg                                                                                                 <  2:2137785/98‑1 (MQ=255)
                                                                                                tCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTtgtg                                                                 >  1:749384/1‑130 (MQ=255)
                                                                                                tCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTtgtg                                                                 <  2:749384/130‑1 (MQ=255)
                                                                                                tCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTgg                                             >  2:1722648/1‑150 (MQ=255)
                                                                                                 cGAGGCGGCCTGGCCCTAGAACTGCGGGTTGGCCTCCAGGGCGCCGGCCAGTTGCCACTGCCCGGCTGCGGTGGCGTTTTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTggg                                            <  1:2160985/150‑1 (MQ=255)
                                                                                                  gAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCg                                                                                                                            <  1:501461/69‑1 (MQ=255)
                                                                                                  gAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCg                                                                                                                            >  2:501461/1‑69 (MQ=255)
                                                                                                  gAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGtt                                                                    >  2:1944132/1‑125 (MQ=255)
                                                                                                  gAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGtt                                                                    <  1:1944132/125‑1 (MQ=255)
                                                                                                      cggcCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCa                                       <  1:2071772/150‑1 (MQ=255)
                                                                                                        gcCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGtt                                                                                      >  2:1498328/1‑101 (MQ=255)
                                                                                                        gcCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGtt                                                                                      <  1:1498328/101‑1 (MQ=255)
                                                                                                               ccTGGAACTGCGGGTTTGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGa                              <  1:1674834/150‑1 (MQ=255)
                                                                                                               ccTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACa                                                            <  2:2238420/120‑1 (MQ=255)
                                                                                                               ccTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACa                                                            >  1:2238420/1‑120 (MQ=255)
                                                                                                               ccTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGa                              >  2:1590659/1‑150 (MQ=255)
                                                                                                                cTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCgttgtt                                                                                                                 >  1:2094284/1‑66 (MQ=255)
                                                                                                                cTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCgttgtt                                                                                                                 <  2:2094284/66‑1 (MQ=255)
                                                                                                                cTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGAt                             >  2:1314699/1‑150 (MQ=255)
                                                                                                                     aCTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGATGACCg                        <  1:598973/150‑1 (MQ=255)
                                                                                                                      cTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAAc                                                                         <  1:1294628/100‑1 (MQ=255)
                                                                                                                      cTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGATGACCgg                       >  1:2346996/1‑150 (MQ=255)
                                                                                                                       tGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAAc                                                                                         <  2:2294956/83‑1 (MQ=255)
                                                                                                                       tGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAAc                                                                                         >  1:2294956/1‑83 (MQ=255)
                                                                                                                       tGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGATGACCGGc                      >  1:2270472/1‑150 (MQ=255)
                                                                                                                           ggTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGATGACCGGCAccc                  >  1:1992397/1‑150 (MQ=255)
                                                                                                                               ggCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGATGACCGGCACCCCGTg              <  2:1757493/150‑1 (MQ=255)
                                                                                                                               ggCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGATGACCGGCACCCCGTg              <  2:175524/150‑1 (MQ=255)
                                                                                                                                gccgccAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGATGACCGGCACCCCGTgg             <  1:1501998/150‑1 (MQ=255)
                                                                                                                                 ccgccAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGATGACCGGCACCCCGTGGc            >  1:1228896/1‑150 (MQ=255)
                                                                                                                                 ccgccAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGATGACCGGCACCCCGTGGc            <  2:237851/150‑1 (MQ=255)
                                                                                                                                   gccAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGATGACCGGCACCCCGTGGCCg          >  1:2285107/1‑150 (MQ=255)
                                                                                                                                    ccAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGATGACCGGCACCCCGTGGCCGt         >  2:268344/1‑150 (MQ=255)
                                                                                                                                    ccAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGATGACCGGCACCCCGTGGCCGt         <  1:250000/150‑1 (MQ=255)
                                                                                                                                      aGGGCGCCGGCGAGTTGGGACTGCCCGGCCGGGGTGGCGTGCTGCACGTCCTGTCCCGGGGTGGCAACGTTGTAGTGGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGATGACCGGCACCCCGTGGCCGTTa       <  2:1129141/150‑1 (MQ=255)
                                                                                                                                       gggCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGTCAGCCCCGTGGCGTTGGGCGGCACGATGTCGATGACCGGCACCCCGTGGCCGTTAt      <  2:1557543/150‑1 (MQ=255)
                                                                                                                                       gggCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGATGACCGGCACCCCGTGGCCGTTAt      >  1:433952/1‑150 (MQ=255)
                                                                                                                                       ccctggcGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGATGACCGGCACCCCGTGGCCGTTAt      >  1:39809/7‑150 (MQ=255)
                                                                                                                                           gcCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGATGACCGGCACCCCGTGGCCGTTATGGAt  >  2:2010374/1‑150 (MQ=255)
                                                                                                                                            ccGGGCAGTCGCAACTGGCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGc                                                  <  2:1812989/101‑1 (MQ=255)
                                                                                                                                            ccGGGCAGTCGCAACTGGCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGc                                                  >  1:1812989/1‑101 (MQ=255)
                                                                                                                                                    ttGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGt                                                                     >  2:1930829/1‑74 (MQ=255)
                                                                                                                                                    ttGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGt                                                                     <  1:1930829/74‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
TGGGGTTGGCCAGTATGTAGTCGGCCGGGCGCCCGAAGATTTCCTGCGGCCCCTCGATCAGCGAGGCGTTGCGGCTGATCACTTCGTTAAGAATGGTCGAGGCGGCCTGGCCCTGGAACTGCGGGTTGGCCGCCAGGGCGCCGGCCAGTTGCGACTGCCCGGCTGCGGTGGCGTTGTTCAGCACCAGGCCCGGGGTGGCAACGTTGTAGTCGATGAACTGGTTGTGTGACAGCCCCGTGGCGTTGGGCGGCACGATGTCGATGACCGGCACCCCGTGGCCGTTATGGAT  >  NC_002947/1656384‑1656672

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