Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,870,260 C→G 6.7% intergenic (+22/‑125) yeaH → / → cdgI DUF444 domain‑containing protein YeaH/putative c‑di‑GMP binding protein CdgI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,870,2600CG 6.7% 488.4 / 18.2 147intergenic (+22/‑125)yeaH/cdgIDUF444 domain‑containing protein YeaH/putative c‑di‑GMP binding protein CdgI
Reads supporting (aligned to +/- strand):  ref base C (79/58);  new base G (7/3);  total (86/61)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

Only 100 of 149 total aligned reads displayed.

GGCGAGAATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTC  >  NC_000913/1870111‑1870407
                                                                                                                                                     |                                                                                                                                                   
ggCGAGAATATGANCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATc                                                                                                                                                     >  1:1030244/1‑150 (MQ=255)
ggCGAGAATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATc                                                                                                                                                     >  2:589171/1‑150 (MQ=255)
 gCGAGAATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCg                                                                                                                                                    >  2:1444741/1‑150 (MQ=255)
  cGAGAATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCCGCCAGGTCATTATCGc                                                                                                                                                   >  2:703487/1‑150 (MQ=255)
  cGAGAATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGc                                                                                                                                                   >  2:1080845/1‑150 (MQ=255)
    agaATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCt                                                                                                                                                 <  2:287120/150‑1 (MQ=255)
      aaTATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTgg                                                                                                                                               <  1:1566574/150‑1 (MQ=255)
        tatGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGCGATAATGACCTGGCt                                                                                                                                             <  2:1386451/150‑1 (MQ=255)
        tatGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGCGATAATGACCTGGCt                                                                                                                                             <  2:1379484/150‑1 (MQ=255)
          tGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGa                                                                                                                                           >  2:330330/1‑150 (MQ=255)
                tCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTa                                                                                                                                     <  1:986441/150‑1 (MQ=255)
                 cTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAg                                                                                                                                    <  2:155573/150‑1 (MQ=255)
                  tGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGc                                                                                                                                   >  2:1693259/1‑150 (MQ=255)
                  tGCAATATACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGATAAAACTATCAGAAAGGTCATTATCGCCTGGCTGATTTTTAga                                                                                                                                   >  2:764224/1‑149 (MQ=255)
                   gCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCt                                                                                                                                  >  2:1331649/1‑150 (MQ=255)
                         tACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGt                                                                                                                            >  2:448754/1‑150 (MQ=255)
                          aCTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTa                                                                                                                           >  2:1896156/1‑150 (MQ=255)
                          aCTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTa                                                                                                                           >  1:213819/1‑150 (MQ=255)
                           cTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATACGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTaa                                                                                                                          >  2:474584/1‑150 (MQ=255)
                            tttCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTg                                                                                                                                            >  2:1339501/1‑131 (MQ=255)
                            tttCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTg                                                                                                                                            <  1:1339501/131‑1 (MQ=255)
                             ttCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAAt                                                                                                                        <  1:1444741/150‑1 (MQ=255)
                             ttCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAAt                                                                                                                        >  2:1172802/1‑150 (MQ=255)
                              tCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAAtt                                                                                                                       <  1:1911053/150‑1 (MQ=255)
                              tCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAAtt                                                                                                                       >  1:501438/1‑150 (MQ=255)
                                gACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTAt                                                                                                                     <  1:1080845/150‑1 (MQ=255)
                                 aCAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTGATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTAtc                                                                                                                    <  2:758534/150‑1 (MQ=255)
                                 aCAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTAtc                                                                                                                    <  2:626153/150‑1 (MQ=255)
                                    aCTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGACTGGCTGATTTTTAGCTTACTGTAAATTATCTCa                                                                                                                 >  1:974835/1‑150 (MQ=255)
                                        tGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATtta                                                                                                             <  1:575256/150‑1 (MQ=255)
                                         gCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATttat                                                                                                            >  1:1200740/1‑150 (MQ=255)
                                          cGATGCAGCACATCCGCGACCAGCCTGATATTTATCCGGTGTTCCGTGAACTGTTTCATAACCAAAATGCAACAGCTAAAGGCTACAACTATCAGCCAGGTCATTACCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATttatt                                                                                                           <  2:1419182/150‑1 (MQ=255)
                                              gcagcaCATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATtacat                                                                                                       <  2:1994781/150‑1 (MQ=255)
                                                 gcaCATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTATCTTACTGTAAATTATCTCATTTATtacataca                                                                                                    >  2:823318/1‑150 (MQ=255)
                                                 gcaCATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATtacataca                                                                                                    >  1:1979947/1‑150 (MQ=255)
                                                  cacaTCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTacatacat                                                                                                   >  2:1524393/1‑150 (MQ=255)
                                                  aacaTCCTCGACCAGGATGATATTTATCCTGTGTTACGTGAACTGTTTCATAAACAAAATGAAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTacatacat                                                                                                   >  2:1197987/2‑150 (MQ=255)
                                                     aTCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAg                                                                                                <  1:1829898/150‑1 (MQ=255)
                                                         gcgACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTa                                                                                            <  2:659255/150‑1 (MQ=255)
                                                          cgACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGGCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTAc                                                                                           <  1:756108/150‑1 (MQ=255)
                                                             ccAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAAt                                                                                        <  1:811956/150‑1 (MQ=255)
                                                                ggATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGc                                                                                     >  2:1890729/1‑150 (MQ=255)
                                                                ggATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGc                                                                                     >  2:1408832/1‑150 (MQ=255)
                                                                   tgatATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCttt                                                                                  <  2:526/150‑1 (MQ=255)
                                                                    gatATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTa                                                                                 >  1:449293/1‑150 (MQ=255)
                                                                     atatTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTa                                                                                                 >  1:1311135/1‑133 (MQ=255)
                                                                     atatTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTa                                                                                                 <  2:1311135/133‑1 (MQ=255)
                                                                     atatTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTaa                                                                                >  2:1508023/1‑150 (MQ=255)
                                                                      tatTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTaaa                                                                               >  1:645999/1‑150 (MQ=255)
                                                                       atTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAat                                                                              >  1:1474384/1‑150 (MQ=255)
                                                                          tATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATg                                                                           <  2:355847/150‑1 (MQ=255)
                                                                              cGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAg                                                                       >  2:347726/1‑150 (MQ=255)
                                                                                 tgtTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAg                                                                       <  1:1328247/147‑1 (MQ=255)
                                                                                     ccGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAAcc                                                                <  1:823318/150‑1 (MQ=255)
                                                                                            cTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACAt                                                         <  2:1979947/150‑1 (MQ=255)
                                                                                            cTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACAt                                                         <  1:1614736/150‑1 (MQ=255)
                                                                                              gTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATaa                                                       >  1:1942620/1‑150 (MQ=255)
                                                                                                ttCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAAtt                                                     >  1:1025375/1‑150 (MQ=255)
                                                                                                ttCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAAtt                                                     >  1:762477/1‑150 (MQ=255)
                                                                                                ttCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAAtt                                                     >  2:1709493/1‑150 (MQ=255)
                                                                                                 tCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAAttt                                                    >  1:1016564/1‑150 (MQ=255)
                                                                                                  cATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGAc                                                                         <  1:633409/128‑1 (MQ=255)
                                                                                                  cATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGAc                                                                         >  2:633409/1‑128 (MQ=255)
                                                                                                   aTAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAg                                                  >  1:1613545/1‑150 (MQ=255)
                                                                                                    tAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGt                                                 >  2:399947/1‑150 (MQ=255)
                                                                                                    tAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGt                                                 >  2:477744/1‑150 (MQ=255)
                                                                                                     aaaCAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGtt                                                >  1:399460/1‑150 (MQ=255)
                                                                                                      aaCAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTc                                               >  2:90354/1‑150 (MQ=255)
                                                                                                           aaTGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGaa                                          >  1:1389980/1‑150 (MQ=255)
                                                                                                              gCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAAc                                       <  2:1435217/150‑1 (MQ=255)
                                                                                                               cAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACa                                      <  1:222384/150‑1 (MQ=255)
                                                                                                               cAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATATCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGAAAGCATAACCTTTACATAATTTAGTTCCAGAAAACa                                      >  1:1575940/1‑150 (MQ=255)
                                                                                                                 aCAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAAt                                    <  2:443737/150‑1 (MQ=255)
                                                                                                                      tAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTc                               <  1:729698/150‑1 (MQ=255)
                                                                                                                       aaaGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCg                              <  2:1025375/150‑1 (MQ=255)
                                                                                                                       aaaGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGATTTAAATATGAAAGCATAACATTTACATAATTTAGTTCAAGAAAACAATCATTCg                              >  2:954199/1‑150 (MQ=255)
                                                                                                                         aGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGa                            <  1:1408832/150‑1 (MQ=255)
                                                                                                                             gctaaaaaTCAGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGaaaaa                        >  2:193784/8‑150 (MQ=39)
                                                                                                                              ctaaaaaTCAGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAAt                       >  2:1925341/7‑150 (MQ=255)
                                                                                                                               aaaCTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATg                      >  1:172344/1‑150 (MQ=255)
                                                                                                                                aaCTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGGTTTATTTATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGa                     <  2:1051113/150‑1 (MQ=255)
                                                                                                                                  aaaTCAGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGAtt                   >  1:674910/3‑150 (MQ=255)
                                                                                                                                   tATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTc                  >  1:1113520/1‑150 (MQ=255)
                                                                                                                                   tATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTc                  >  1:1120532/1‑150 (MQ=255)
                                                                                                                                     tCAGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAg                >  1:1073380/1‑150 (MQ=255)
                                                                                                                                     tCAGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAg                >  1:376688/1‑150 (MQ=255)
                                                                                                                                     tCAGCCAGGCGATAATGACCTGGCTGATTTTTAGATTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATGGCTTTAAATATGACATCATAACATTTACATAATTTAGTAAAGGAAACAATAATTCGGAAAAattattacgt               >  2:851056/1‑146 (MQ=255)
                                                                                                                                      cAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGt               <  1:474584/150‑1 (MQ=255)
                                                                                                                                      cAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGt               >  1:936641/1‑150 (MQ=255)
                                                                                                                                      cAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGt               >  2:1170186/1‑150 (MQ=255)
                                                                                                                                       aGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTc              <  2:501438/150‑1 (MQ=255)
                                                                                                                                       aGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTc              <  2:601823/150‑1 (MQ=255)
                                                                                                                                       aGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTc              <  1:1350062/150‑1 (MQ=255)
                                                                                                                                        gCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCa             >  2:419505/1‑150 (MQ=255)
                                                                                                                                         ccAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCaa            <  1:954199/150‑1 (MQ=255)
                                                                                                                                           aGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAaca          <  1:1182267/150‑1 (MQ=255)
                                                                                                                                           aGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAaca          <  2:1613545/150‑1 (MQ=255)
                                                                                                                                                  tATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTAttt   >  1:512458/1‑150 (MQ=255)
                                                                                                                                                  tATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTAttt   >  1:552904/1‑150 (MQ=255)
                                                                                                                                                  tATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTAttt   <  1:1709493/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGCGAGAATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTC  >  NC_000913/1870111‑1870407

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