Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,895,071 A→G 7.6% L11P (CTA→CCA)  yegH ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,895,0710AG 7.6% 229.6 / 6.5 79L11P (CTA→CCA) yegHputative transporter
Reads supporting (aligned to +/- strand):  ref base A (34/39);  new base G (3/3);  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.49e-01

CCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTGATAACAAATTGGGAGAAC  >  NC_002947/3894925‑3895219
                                                                                                                                                  |                                                                                                                                                    
ccAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTGGct                                                                                                                                                   >  2:732537/1‑149 (MQ=255)
       cATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGt                                                                                                                                            <  1:27422/150‑1 (MQ=255)
        atgatgAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTGGCTAGGcctgt                                                                                                                                           <  1:732537/150‑4 (MQ=255)
         tgatgaGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTgg                                                                                                                                          <  2:1271109/150‑1 (MQ=255)
         tgatgaGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTgg                                                                                                                                          >  2:519909/1‑150 (MQ=255)
                  ccAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTa                                                                                                                                 <  1:141497/150‑1 (MQ=255)
                      gCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGcca                                                                                                                             <  2:1208804/150‑1 (MQ=255)
                       cTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATAAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATGAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGccat                                                                                                                            >  2:1180416/1‑150 (MQ=255)
                        tCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGccatt                                                                                                                           >  1:865609/1‑150 (MQ=255)
                         cAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGccattc                                                                                                                          <  2:1323653/150‑1 (MQ=255)
                         cAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGccattc                                                                                                                          <  2:720823/150‑1 (MQ=255)
                           ggCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTGGCTAGGcctgt                                                                                                                                           <  1:293486/131‑4 (MQ=255)
                           ggCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTGGCTAGGcctgt                                                                                                                                           >  2:293486/1‑128 (MQ=255)
                              cAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTa                                                                                                                     >  1:167813/1‑150 (MQ=255)
                              cAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTa                                                                                                                     <  1:1631463/150‑1 (MQ=255)
                              cAATGACCCGCGCGCGGTCGAGATGATGGGGCTGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTa                                                                                                                     <  2:1614297/150‑1 (MQ=255)
                                aTGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGt                                                                                                                   <  2:167813/150‑1 (MQ=255)
                                  gACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTgg                                                                                                                 <  1:88436/150‑1 (MQ=255)
                                        cgcgcgGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTc                                                                                                           <  2:54461/150‑1 (MQ=255)
                                         gcgcgGTCGCGCTGATGGGGAGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATACGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCa                                                                                                          >  1:269317/1‑150 (MQ=255)
                                          cgcgGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCAc                                                                                                         <  1:591072/150‑1 (MQ=255)
                                                     tGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATGAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGt                                                                                              <  1:1180416/150‑1 (MQ=255)
                                                       aTGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTca                                                                                            >  2:109751/1‑150 (MQ=255)
                                                         ggggCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGt                                                                                                                   <  1:1552728/125‑1 (MQ=255)
                                                         ggggCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGt                                                                                                                   >  2:1552728/1‑125 (MQ=255)
                                                              ggcAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGccattc                                                                                                                          >  1:320261/1‑113 (MQ=255)
                                                              ggcAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGccattc                                                                                                                          <  2:320261/113‑1 (MQ=255)
                                                                 aaCTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACggg                                                                                  <  1:111269/150‑1 (MQ=255)
                                                                         ggCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGAt                                                                          <  1:519909/150‑1 (MQ=255)
                                                                          gCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGACTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATg                                                                         >  1:845626/1‑150 (MQ=255)
                                                                               gATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAg                                                                                                                                                 <  1:293700/73‑1 (MQ=255)
                                                                               gATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAg                                                                                                                                                 >  2:293700/1‑73 (MQ=255)
                                                                                  ggCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGaaa                                                                 <  1:4076/150‑1 (MQ=255)
                                                                                     gATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTgg                                                                                                     >  2:1424154/1‑111 (MQ=255)
                                                                                     gATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTgg                                                                                                     <  1:1424154/111‑1 (MQ=255)
                                                                                            acCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGACAGAAAGGCCGAAGAc                                                       >  1:1486288/1‑150 (MQ=255)
                                                                                            acCAGGTTGTCGATAACCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGAc                                                       >  1:1365944/1‑150 (MQ=255)
                                                                                             ccAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACt                                                      >  2:1377328/1‑150 (MQ=255)
                                                                                              cAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTg                                                     >  1:1555740/1‑150 (MQ=255)
                                                                                                    gTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGc                                               >  2:1422903/1‑150 (MQ=255)
                                                                                                    gTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGc                                               >  1:719986/1‑150 (MQ=255)
                                                                                                     tCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCAAGACCGGCAGTACACTGGGGGGcc                                              >  2:1583507/1‑149 (MQ=255)
                                                                                                       gATGCCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGGCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTcc                                            >  2:884008/1‑150 (MQ=255)
                                                                                                        aTACCCAGCACCAGCTCGAGGGCGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACgggg                                                                                 <  1:373863/112‑1 (MQ=255)
                                                                                                        aTACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACgggg                                                                                 >  2:373863/1‑112 (MQ=255)
                                                                                                        aTACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGaaa                                                                 >  1:788332/1‑128 (MQ=255)
                                                                                                        aTACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGaaa                                                                 <  2:788332/128‑1 (MQ=255)
                                                                                                         tACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGGGGTATGGCCAGAAAGGCCGAAGACTGGGGGGCTCcgc                                          >  2:1620635/1‑150 (MQ=255)
                                                                                                            ccAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGaaa                                                                 >  2:649868/1‑124 (MQ=255)
                                                                                                            ccAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGaaa                                                                 <  1:649868/124‑1 (MQ=255)
                                                                                                            ccAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGaa                                       >  1:1127233/1‑150 (MQ=255)
                                                                                                                 aCCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACggg                                                                                  <  2:1237931/102‑1 (MQ=255)
                                                                                                                 aCCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACggg                                                                                  >  1:1237931/1‑102 (MQ=255)
                                                                                                                   cAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCCTTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTc                                <  1:884008/150‑1 (MQ=255)
                                                                                                                    aGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCa                               <  1:1422903/150‑1 (MQ=255)
                                                                                                                      cTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGTCGAAGACTGGGGGGCTCCGCGAAGGTGCTCata                             <  1:1620635/150‑1 (MQ=255)
                                                                                                                       tCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCatat                            >  2:1462139/1‑150 (MQ=255)
                                                                                                                          aGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGa                                                        >  1:409828/1‑119 (MQ=255)
                                                                                                                          aGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGa                                                        <  2:409828/119‑1 (MQ=255)
                                                                                                                           ggACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGAc                        <  1:1462139/150‑1 (MQ=255)
                                                                                                                              cGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCtt                     <  2:1555740/150‑1 (MQ=255)
                                                                                                                               gATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTg                    <  2:538164/150‑1 (MQ=255)
                                                                                                                                 tAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTc                                                                                             <  1:535697/75‑1 (MQ=255)
                                                                                                                                 tAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTc                                                                                             >  2:535697/1‑75 (MQ=255)
                                                                                                                                 tAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCata                             >  1:923514/1‑139 (MQ=255)
                                                                                                                                 tAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCata                             <  2:923514/139‑1 (MQ=255)
                                                                                                                                   aGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTGATaa                >  2:390989/1‑150 (MQ=255)
                                                                                                                                     cGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTGATAACa              >  2:938850/1‑150 (MQ=255)
                                                                                                                                       tCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGaa                                       >  1:888670/1‑123 (MQ=255)
                                                                                                                                       tCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGaa                                       <  2:888670/123‑1 (MQ=255)
                                                                                                                                          cacggCCTGGCTAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGaa                                                                              <  2:310032/78‑1 (MQ=255)
                                                                                                                                          cacggCCTGGCTAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGaa                                                                              >  1:310032/4‑81 (MQ=255)
                                                                                                                                          aCAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTGATAACAAAttt         >  2:302377/1‑149 (MQ=255)
                                                                                                                                           cAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTGATAACAAATTgg        <  1:1377328/150‑1 (MQ=255)
                                                                                                                                             ggCCTAGCCAGGCCGTGGGGTCGACTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTGATAACAAATTGGga      <  2:845626/150‑1 (MQ=255)
                                                                                                                                                 tAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCtt                     >  2:19735/1‑131 (MQ=255)
                                                                                                                                                 tAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCtt                     <  1:19735/131‑1 (MQ=255)
                                                                                                                                                 tAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTGATAACAAATTGGGAGAAc  <  1:938850/150‑1 (MQ=255)
                                                                                                                                                 tAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTGATAACAAATTGGGAGAAc  <  1:109751/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTGATAACAAATTGGGAGAAC  >  NC_002947/3894925‑3895219

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