Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,674,958 A→T 5.9% intergenic (+11/+14) tqsA → / ← pntB pheromone AI‑2 transporter/pyridine nucleotide transhydrogenase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,674,9580AT 5.9% 193.7 / 2.4 68intergenic (+11/+14)tqsA/pntBpheromone AI‑2 transporter/pyridine nucleotide transhydrogenase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (34/29);  new base T (2/2);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

GATTTTTTGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAATCACG  >  NC_000913/1674809‑1675104
                                                                                                                                                     |                                                                                                                                                  
gATTTTTTGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTcagca                                                                                                                                                    >  1:1569507/1‑150 (MQ=255)
 aTTTTTTGGGGATGGTTGTTAGGCCCGGTGGGTATGCTGCTTTCAGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCATAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTcagcag                                                                                                                                                   >  1:1535922/1‑150 (MQ=255)
 aTTTTTTGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTGCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTcagcag                                                                                                                                                   <  1:505611/150‑1 (MQ=255)
   tttttGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCagag                                                                                                                                                 >  1:391904/1‑150 (MQ=255)
   tttttGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTGTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCagag                                                                                                                                                 <  2:720302/150‑1 (MQ=255)
     tttGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGGCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGc                                                                                                                                               <  1:917647/150‑1 (MQ=255)
           gATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCTGCTGAGGCCGTCac                                                                                                                                         <  1:1659845/150‑2 (MQ=255)
           gATGGGTGGTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCCTGAACAAACAGCGGGAGGGCAAAGCATCGCCGATCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCTGCTGAGGCCGTCac                                                                                                                                         <  1:1038766/150‑2 (MQ=255)
             tGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAggg                                                                                                                                       >  1:811211/1‑150 (MQ=255)
             tGGTTGTTAGGACCGGGGGGTATGCTGCGTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAggg                                                                                                                                       <  1:869281/150‑1 (MQ=255)
               gttgttAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGtt                                                                                                                                     >  2:1616644/1‑150 (MQ=255)
               gttgttAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGtt                                                                                                                                     >  2:959005/1‑150 (MQ=255)
               gttgttAGGACCGGTGGGTATGCTGCTTTCCGAGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGtt                                                                                                                                     <  2:184996/150‑1 (MQ=255)
                    tAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCATTAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACaga                                                                                                                                >  1:1397887/1‑150 (MQ=255)
                    tAGGACCGGTGGCTAAGCTGCTTTCCTTGCCGTTGACAATTATTGTAAAAATTGCGCTTGAACAAACAGCGGGAGGTCACAGCATCGACGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGACTGGGCCGTCAGGGTTACaca                                                                                                                                >  2:483737/1‑148 (MQ=255)
                     aGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACagag                                                                                                                               >  1:1680048/1‑150 (MQ=255)
                        aCCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCtt                                                                                                                            <  1:13232/150‑1 (MQ=255)
                               ggTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGAt                                                                                                                     >  2:377255/1‑150 (MQ=255)
                                 tATGCTGCTTTCCGTGTCGTTTACAATTTTTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGGATCGCCGTTCTGTTAAGGTATCTCAATATAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTg                                                                                                                   <  2:1701977/150‑1 (MQ=255)
                                 tATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGATGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTg                                                                                                                   <  2:1729286/150‑1 (MQ=255)
                                 tATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTg                                                                                                                   <  1:1568295/150‑1 (MQ=255)
                                 tATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTg                                                                                                                   <  2:1056923/150‑1 (MQ=255)
                                   tgctgcTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGGGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCa                                                                                                                 <  2:627709/150‑1 (MQ=255)
                                     ctgctTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCCGGATTGCATc                                                                                                               >  2:1235521/1‑150 (MQ=255)
                                     ctgctTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCATGATTGCATc                                                                                                               >  2:1469555/1‑150 (MQ=255)
                                     ctgctTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATc                                                                                                               >  2:1137611/1‑150 (MQ=255)
                                      tgctTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGTGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATcc                                                                                                              <  1:1616807/150‑1 (MQ=255)
                                      tgctTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGTGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATcc                                                                                                              <  1:1602987/150‑1 (MQ=255)
                                             cGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTgg                                                                                                       >  2:746741/1‑150 (MQ=255)
                                                   gtcgACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTgg                                                                                                 >  1:547672/4‑150 (MQ=255)
                                                   gtagACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTgg                                                                                                 >  1:547668/4‑150 (MQ=255)
                                                   gTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGCGCTTTCAGGATTGCATCCACGCTGGCTTTgg                                                                                                 >  2:1146516/1‑150 (MQ=255)
                                                    ttGACAATTATTGTCAAAATTGCTCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTTTTAAGCGATCTCAATAATGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGc                                                                                                <  1:1400638/150‑1 (MQ=255)
                                                        cAATTATTGTCAAAATTGCGCTTGAACAAACAGCGCGAGGTCAAAGCATGGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCACAGGCCGGCAGCGATACAGCGCGTACAGGATTGCCTACCCCCTGGCTGTGGCGTCa                                                                                            >  1:299526/1‑150 (MQ=255)
                                                              ttGTCTAAATTGCGCTTGAACTAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAAc                                                                                      <  1:420199/150‑1 (MQ=255)
                                                                gTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCATTAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAg                                                                                    >  2:829828/1‑150 (MQ=255)
                                                                gTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAg                                                                                    >  2:1555546/1‑150 (MQ=255)
                                                                     aaTTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCAtgt                                                                               <  2:1371113/150‑1 (MQ=255)
                                                                      aTTGCGCTTGAAAAAACAGCGGGATGTCTAAGAATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCTCCAGATGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCAtgtg                                                                              <  1:864905/150‑1 (MQ=255)
                                                                        tGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTggg                                                                            >  1:1492949/1‑150 (MQ=255)
                                                                              tGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGtttt                                                                      >  2:689144/1‑150 (MQ=255)
                                                                               gAACAAACCGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCCGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGGGGGTGTTTTc                                                                     >  2:134178/1‑150 (MQ=255)
                                                                                aacaaacaGCGGGAGGTCAAAGCATCGCCGTTCTGATAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGCTGTTTTcc                                                                    >  1:599543/1‑150 (MQ=255)
                                                                                 acaaacaGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCt                                                                   <  1:1616644/150‑1 (MQ=255)
                                                                                 acaaacaGCGGGAGGTCAAAGCAGCGCCGAGCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCt                                                                   <  1:173005/150‑1 (MQ=255)
                                                                                      caGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTgaag                                                              <  1:407539/150‑1 (MQ=255)
                                                                                              ggTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTAGAAGAACAGCgg                                                      >  1:703064/1‑150 (MQ=255)
                                                                                               gTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCggg                                                     >  2:1474176/1‑150 (MQ=255)
                                                                                                 cAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGtt                                                   >  1:1146009/1‑150 (MQ=255)
                                                                                                 cAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGAATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTACTTGAAGAACAGCGGGtt                                                   >  2:529478/1‑150 (MQ=255)
                                                                                                    cgCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTATCCGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTg                                                <  2:351588/149‑1 (MQ=255)
                                                                                                     gCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCGGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGc                                               >  2:278458/1‑150 (MQ=255)
                                                                                                      cATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGGCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGACGAACAGCGGGTTTTGca                                              <  1:1055579/150‑1 (MQ=255)
                                                                                                       aTCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGcac                                             <  1:171673/150‑1 (MQ=255)
                                                                                                           ccGTTCTGTTAAGAGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCa                                         <  1:36807/150‑1 (MQ=255)
                                                                                                            cGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAg                                        >  1:1400511/1‑150 (MQ=255)
                                                                                                               tCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCAt                                     >  2:1444600/1‑150 (MQ=255)
                                                                                                               tCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCAt                                     >  1:585354/1‑150 (MQ=255)
                                                                                                                 tGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAg                                   <  1:959005/150‑1 (MQ=255)
                                                                                                                    tAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACCGCGGGTTTTGCACACCAGCATCGCCa                                >  1:133250/1‑150 (MQ=255)
                                                                                                                    tAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGACAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCTAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCa                                >  1:127999/1‑150 (MQ=255)
                                                                                                                     aaGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAg                               <  2:587799/150‑1 (MQ=255)
                                                                                                                        cGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAgtgt                            >  2:556284/1‑150 (MQ=255)
                                                                                                                           tctcAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCa                         <  2:586240/150‑1 (MQ=255)
                                                                                                                             tcATTAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATc                       <  1:829828/150‑1 (MQ=255)
                                                                                                                              cAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGATCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAAAAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCg                      >  2:222765/1‑150 (MQ=255)
                                                                                                                                     ccctGACGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGttt               >  2:223161/4‑150 (MQ=255)
                                                                                                                                      agTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTa              <  1:762225/150‑1 (MQ=255)
                                                                                                                                         gACGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAg           >  1:1410266/1‑150 (MQ=255)
                                                                                                                                         gACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAg           <  2:391904/150‑1 (MQ=255)
                                                                                                                                         gACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAg           >  1:970388/1‑150 (MQ=255)
                                                                                                                                           cGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGGTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGATCAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGAc         <  2:1397218/150‑1 (MQ=255)
                                                                                                                                           cGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCTGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGAc         >  1:903626/1‑150 (MQ=255)
                                                                                                                                                  agcagAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAATCACg  <  2:1400511/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GATTTTTTGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAATCACG  >  NC_000913/1674809‑1675104

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