Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,067,747 T→A 7.7% intergenic (‑39/‑334) wrbA ← / → ymdF NAD(P)H:quinone oxidoreductase/conserved protein YmdF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,067,7470TA 7.7% 265.3 / 8.5 88intergenic (‑39/‑334)wrbA/ymdFNAD(P)H:quinone oxidoreductase/conserved protein YmdF
Reads supporting (aligned to +/- strand):  ref base T (45/36);  new base A (2/5);  total (47/41)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

TAACGACAACTTCAGCGCCATCCACTTTGCTTGCACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATCGCGTTTTTTTAATCATAAGAG  >  NC_000913/1067600‑1067894
                                                                                                                                                   |                                                                                                                                                   
tAACGACAACTTCAGCGCCATCCACTTTGCTTGCACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTc                                                                                                                                                   >  2:485631/1‑150 (MQ=255)
  aCGACAACTTCAGCGCCATCCACTTTGCTTGCACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGc                                                                                                                                                 <  2:779481/150‑1 (MQ=255)
  aCGACAACATCATCGCCATCCACTATGCTAGCACCCTCAGCGACTGCGCGTCCCATCCTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGc                                                                                                                                                 >  2:1234708/1‑150 (MQ=255)
   cGACAACTTCAGCGCCATCCACTTTGCTTGCACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCt                                                                                                                                                <  2:1850885/150‑1 (MQ=255)
     aCAACTTCAGCGCCATCCACTTTGCTTGCACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTAc                                                                                                                                              <  2:568487/150‑1 (MQ=255)
        aCTTCAGCGCCATCCACTTTGCTTGCACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCAc                                                                                                                                           <  1:597937/150‑1 (MQ=255)
        aCTTCAGCGCCATCCACTTTGCTTGCACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGAATCGCTAcgga                                                                                                                                           <  1:230896/150‑4 (MQ=255)
         cTTCAGCGCCATCCACTTTGCTTGCACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGACATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCACTACCACt                                                                                                                                          >  1:877419/1‑150 (MQ=255)
          ttCAGCGCCATCCACTTTGCTTGCACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACtt                                                                                                                                         >  1:1925819/1‑150 (MQ=255)
          ttCAGCGCCATCCACTTTGCTTGCACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGAATCGCTAcggaat                                                                                                                                         <  1:667242/150‑6 (MQ=255)
           tCAGCGCCATCCACTTTGCTTGCACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATTGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTAATCCGTAGCGAATCGCTAcggaata                                                                                                                                        <  1:1899824/150‑7 (MQ=255)
           tCAGCGCCATCCACTTTGCTTGCACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGAATCGCTAcggaata                                                                                                                                        <  2:1200845/150‑7 (MQ=255)
             aGCGCCATCCACTTTGCTTGCACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTAtt                                                                                                                                      >  1:1302022/1‑150 (MQ=255)
                gcCATCCACTTTGCTTGCACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTaa                                                                                                                                   <  1:1075992/150‑1 (MQ=255)
                gcCATCCACTTTGCTTGCACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTaa                                                                                                                                   >  2:312766/1‑150 (MQ=255)
                  cATCCACTTTGCTTGCACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAg                                                                                                                                 <  2:785124/150‑1 (MQ=255)
                    tCCACTTTGCTTGCACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACAT‑GAATAATAAAGCACCACAACTTTTGCCCTTTGTTACCCCGGCGAGTGTTAAGGCTATTCCGTAGCGAATCGCGACCACATAATTAAAGATa                                                                                                                              >  1:812361/1‑150 (MQ=255)
                               tgcACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCtt                                                                                                                    >  1:553253/1‑150 (MQ=255)
                                gcACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCttt                                                                                                                   >  2:1951554/1‑150 (MQ=255)
                                gcACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCttt                                                                                                                   <  2:1090473/150‑1 (MQ=255)
                                  aCCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTCTTc                                                                                                                 <  1:528247/150‑1 (MQ=255)
                                   cccTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCa                                                                                                                >  1:1447648/1‑150 (MQ=255)
                                    ccTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCag                                                                                                               <  1:1204476/150‑1 (MQ=255)
                                     cTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCaga                                                                                                              >  1:834110/1‑150 (MQ=255)
                                     cTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCaga                                                                                                              <  1:1271568/150‑1 (MQ=255)
                                      tCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCagag                                                                                                             <  2:8872/150‑1 (MQ=255)
                                        aGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTg                                                                                                           <  1:1530295/150‑1 (MQ=255)
                                            aCTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCTATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCaaa                                                                                                       <  2:296702/150‑1 (MQ=255)
                                              tGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAAtt                                                                                                     >  1:1198824/1‑150 (MQ=255)
                                              tGCGCGTGCCATAGTTTCAATATGTCCGTACATGGAATAATAAATCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAAtt                                                                                                     <  1:833891/150‑1 (MQ=255)
                                               gcgcGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAAttt                                                                                                    >  1:1557513/1‑150 (MQ=255)
                                                cgcgTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTc                                                                                                   <  2:211731/150‑1 (MQ=255)
                                                    tGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACaa                                                                                               >  2:1111238/1‑150 (MQ=255)
                                                     gCCATCGTTTCAATATGTCCTTACAT‑GAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAAcc                                                                                             <  2:812361/150‑1 (MQ=255)
                                                            tttCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTAt                                                                                       >  2:1398293/1‑150 (MQ=255)
                                                              tCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTAttt                                                                                     >  2:818537/1‑150 (MQ=255)
                                                               cAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTg                                                                                    >  1:1702848/1‑150 (MQ=255)
                                                                    tGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTta                                                                               >  1:199356/1‑150 (MQ=255)
                                                                     gTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTtat                                                                              <  2:1720008/150‑1 (MQ=255)
                                                                      tCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTtata                                                                             <  2:916736/150‑1 (MQ=255)
                                                                      tCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATATCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTtata                                                                             <  1:1243911/150‑1 (MQ=255)
                                                                        cGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATaac                                                                           >  2:1867048/1‑150 (MQ=255)
                                                                        cGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATaac                                                                           >  1:1558794/1‑150 (MQ=255)
                                                                          tACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATaacaa                                                                         >  1:596250/1‑150 (MQ=255)
                                                                          tACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATaacaa                                                                         >  1:1621197/1‑150 (MQ=255)
                                                                           aCATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATaacaac                                                                        >  1:1707484/1‑150 (MQ=255)
                                                                              tGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACttt                                                                     >  2:1908035/1‑150 (MQ=255)
                                                                                gAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTca                                                                   >  1:860611/1‑150 (MQ=255)
                                                                                   taataaAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACaa                                                                >  1:1501948/1‑150 (MQ=255)
                                                                                   taataaAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACaa                                                                <  1:1951554/150‑1 (MQ=255)
                                                                                     ataaAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGc                                                              <  1:82984/150‑1 (MQ=255)
                                                                                        aaGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACg                                                           <  2:238466/150‑1 (MQ=255)
                                                                                        aaGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACg                                                           >  1:1823851/1‑150 (MQ=255)
                                                                                        aaGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACg                                                           >  2:1436109/1‑150 (MQ=255)
                                                                                           cACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTaa                                                        >  2:576081/1‑150 (MQ=255)
                                                                                           cACCAGAACTTTAGCCATTTCTAACCACTCCTAGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTaa                                                        >  2:74256/1‑150 (MQ=255)
                                                                                            aCCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAAt                                                       >  1:742695/1‑150 (MQ=255)
                                                                                              cAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAAttt                                                     >  1:260506/1‑150 (MQ=255)
                                                                                                     ttAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCa                                              <  1:487897/150‑1 (MQ=255)
                                                                                                      tAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGGAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCaa                                             >  1:162734/1‑150 (MQ=255)
                                                                                                       aGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCaaa                                            >  2:1500400/1‑150 (MQ=255)
                                                                                                       aGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCaaa                                            <  1:1398293/150‑1 (MQ=255)
                                                                                                        gCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCaaaa                                           >  2:1836607/1‑150 (MQ=255)
                                                                                                        gCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCaaaa                                           >  2:1891420/1‑150 (MQ=255)
                                                                                                                 aaCCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACAtt                                  <  2:860611/150‑1 (MQ=255)
                                                                                                                       tcctcGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTAtc                            <  1:818537/150‑1 (MQ=255)
                                                                                                                        cctcGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTa                              >  2:1550636/1‑147 (MQ=255)
                                                                                                                        cctcGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTa                              <  1:1550636/147‑1 (MQ=255)
                                                                                                                             tgtTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCCTcg                      <  1:1891420/150‑1 (MQ=255)
                                                                                                                             tgtTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATcg                      >  1:147930/1‑150 (MQ=255)
                                                                                                                             tgtTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATcg                      >  2:1677251/1‑150 (MQ=255)
                                                                                                                                   ctctATTCCTTATCGATTCGATACCACTTATTTAAAGATAAGAAGTCATTTTCAGAGTGCAAATTTAACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATCGCGtttt                >  2:1595058/1‑150 (MQ=255)
                                                                                                                                   ctctATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATCGCGtttt                >  2:280930/1‑150 (MQ=255)
                                                                                                                                      tATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATCGCGttttttt             <  1:1111238/150‑1 (MQ=255)
                                                                                                                                      tATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATCGCGttttttt             <  2:1302022/150‑1 (MQ=255)
                                                                                                                                       aTTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATCGCGTTTTTTTa            <  1:74256/150‑1 (MQ=255)
                                                                                                                                       aTTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATCGCGTTTTTTTa            <  1:1500400/150‑1 (MQ=255)
                                                                                                                                         gtggTAGCGAATCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATCGCGTTTTTTTAAt          >  1:774061/4‑150 (MQ=255)
                                                                                                                                          ccGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATCGCGTTTTTTTAATc         <  2:596250/150‑1 (MQ=255)
                                                                                                                                           ggTAGCGAATCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATCGCGTTTTTTTAATCa        <  2:774061/149‑1 (MQ=255)
                                                                                                                                           cGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATCGCGTTTTTTTAATCa        <  2:1116798/150‑1 (MQ=255)
                                                                                                                                             tAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATCGCGTTTTTTTAATCATa      <  1:1836607/150‑1 (MQ=255)
                                                                                                                                               gCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATCGCGTTTTTTTAATCATAag    >  2:1215484/1‑150 (MQ=255)
                                                                                                                                                 gATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATCGCGTTTTTTTAATCATAagag  >  1:1597833/1‑150 (MQ=255)
                                                                                                                                                 gATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATCGCGTTTTTTTAATCATAagag  >  1:1411434/1‑150 (MQ=255)
                                                                                                                                                 gATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATCGCGTTTTTTTAATCATAagag  <  2:1707484/150‑1 (MQ=255)
                                                                                                                                                 gATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATCGCGTTTTTTTAATCATAagag  >  1:661610/1‑150 (MQ=255)
                                                                                                                                                 gATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATCGCGTTTTTTTAATCATAagag  >  1:110135/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TAACGACAACTTCAGCGCCATCCACTTTGCTTGCACCCTCAGCGACTGCGCGTGCCATCGTTTCAATATGTCCGTACATGGAATAATAAAGCACCAGAACTTTAGCCATTTCTAACCACTCCTCGTGTTATCTCTATTCCGTAGCGATTCGCTACCACTTATTTAAAGATAAGACGTCCTTTTCAGAGTGCAAATTTCACAACCACTTATTTGATTTATAACAACTTTCACAAGCACGTAATTTTGTCGCAAAATGACACATTTTTATCTCATCGCGTTTTTTTAATCATAAGAG  >  NC_000913/1067600‑1067894

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