Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 89,066 C→A 5.6% intergenic (+34/‑568) cra → / → mraZ transcriptional repressor‑activator for carbon metabolism/RsmH methytransferase inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091389,0660CA 5.6% 325.5 / 5.3 109intergenic (+34/‑568)cra/mraZtranscriptional repressor‑activator for carbon metabolism/RsmH methytransferase inhibitor
Reads supporting (aligned to +/- strand):  ref base C (51/51);  new base A (2/5);  total (53/56)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 113 total aligned reads displayed.

CAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGT  >  NC_000913/88917‑89215
                                                                                                                                                      |                                                                                                                                                     
cAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTcc                                                                                                                                                       <  2:530131/150‑1 (MQ=255)
 agagCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTccc                                                                                                                                                      <  2:342813/150‑1 (MQ=255)
       tgtgCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGgc                                                                                                                                                >  1:782636/1‑150 (MQ=255)
       tgtgCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGgc                                                                                                                                                <  1:1659108/150‑1 (MQ=255)
        gtgCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGgcg                                                                                                                                               <  2:1413010/150‑1 (MQ=255)
          gCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCGAAGCCGCGATCAAAAATACGCGCCAGGTGAATTTCCCTCTGGcgcg                                                                                                                                             >  1:853035/1‑150 (MQ=255)
           cTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGt                                                                                                                                            >  2:1261004/1‑150 (MQ=255)
            tGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTa                                                                                                                                           >  1:966781/1‑150 (MQ=255)
             ggAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGAGGGAAATTCACCTGGCGCGTat                                                                                                                                          <  2:1170973/150‑2 (MQ=255)
             ggAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTag                                                                                                                                          <  2:1359887/150‑1 (MQ=255)
                gttcgTCCAGGCTTGCCAGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGt                                                                                                                                       >  1:596561/5‑150 (MQ=255)
                gttcgTCCAGGCTTGCCAGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGt                                                                                                                                       >  1:1105149/5‑150 (MQ=255)
                gaTTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGAGGGAAATTCACCTGGCGCGTatttt                                                                                                                                       <  2:1368984/150‑5 (MQ=38)
                gaTTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGAGGGAAATTCACCTGGCGCGTatttt                                                                                                                                       <  2:848856/150‑5 (MQ=38)
                gaTTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGt                                                                                                                                       <  2:502468/150‑1 (MQ=255)
                   tGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACg                                                                                                                                    <  1:815273/150‑1 (MQ=255)
                    gTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACgg                                                                                                                                   <  1:1624944/150‑1 (MQ=255)
                    gTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACgg                                                                                                                                   <  2:1571803/150‑1 (MQ=255)
                    gTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACgg                                                                                                                                   <  1:848397/150‑1 (MQ=255)
                      ccTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGa                                                                                                                                 >  1:1703272/1‑150 (MQ=255)
                      ccTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGa                                                                                                                                 >  1:356230/1‑150 (MQ=255)
                      ccTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGa                                                                                                                                 >  1:509958/1‑150 (MQ=255)
                       cTGGCCAGCCTGGACGAACCGCGGAAGCCAAAACCTGGTTTAACGCACATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGCAGCTAAGCCGCGAACAAAAATACGCGCCAGGGGAATTTCCCTCTGGCGCGTAGAGTACGGGAc                                                                                                                                <  1:1130029/150‑1 (MQ=255)
                        tGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACt                                                                                                                               >  1:1513207/1‑150 (MQ=255)
                          gCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAAGCTCTATCGCCGCGGCGGGCTCCGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTgg                                                                                                                             <  1:1424933/150‑1 (MQ=255)
                               ccTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATc                                                                                                                        >  2:1157487/1‑150 (MQ=255)
                                cTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCa                                                                                                                       >  2:114019/1‑150 (MQ=255)
                                cTGGACGAACCGCGTAAGCCAAAACCTCGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCAGTACCTAAGCCGCGAAAAAAAATACGAGACAGGTGAATTTCCCTCTGGCGAGAAGACTACGGGCCTGGACATCa                                                                                                                       >  2:329865/1‑150 (MQ=255)
                                  ggccGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAat                                                                                                                     >  1:1731332/4‑150 (MQ=255)
                                   gACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAata                                                                                                                    <  2:663879/150‑1 (MQ=255)
                                      gAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGc                                                                                                                 >  1:1028636/1‑150 (MQ=255)
                                       aaCCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCt                                                                                                                <  1:811925/150‑1 (MQ=255)
                                        aCCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCtt                                                                                                               >  1:987863/1‑150 (MQ=255)
                                         ccGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGGGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTa                                                                                                              <  1:734699/150‑1 (MQ=255)
                                            cgTAAGCCAAAACCTGGTATAACGCGCATAAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCTAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAg                                                                                                           <  1:1543140/150‑1 (MQ=255)
                                              tAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACTGGACAGAACATCAATATGCTTAAAGTa                                                                                                         >  1:322255/1‑150 (MQ=255)
                                               aaGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTaa                                                                                                        <  1:1647305/150‑1 (MQ=255)
                                                aGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTAAAAGTaaa                                                                                                       >  2:1167692/1‑150 (MQ=255)
                                                 gCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGAGGGAAATTCACCTGGCGCGTa                                                                                                                                           <  2:898181/113‑1 (MQ=38)
                                                 gCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGAGGGAAATTCACCTGGCGCGTa                                                                                                                                           >  1:898181/1‑113 (MQ=38)
                                                 gCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAAt                                                                                                      >  2:1269806/1‑150 (MQ=255)
                                                 gCCAAAACCTGGATTAACGCGCATTAAACGTAATCTCTATCGCCGCCGCGTGCTAAGCCGAAGCTAAGCCGCGATCAAAAATACGCGCCAAGAGAATTTCCCTCTGGCGCGAAGAGTACGGGAAAGGACATCAATATGCTTAAAGTaaaa                                                                                                      >  1:828725/1‑149 (MQ=255)
                                                 cccAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAAt                                                                                                      <  1:455011/149‑1 (MQ=255)
                                                  ccAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATa                                                                                                     >  2:1643944/1‑150 (MQ=255)
                                                  ccAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATa                                                                                                     <  2:689642/150‑1 (MQ=255)
                                                  ccAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATa                                                                                                     >  2:798352/1‑150 (MQ=255)
                                                    aaaaCCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATCCGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAg                                                                                                   <  1:1496064/150‑1 (MQ=255)
                                                    aaaaCCTGGTTTAACGCGCATTAAACGTAATCTCTAGCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAg                                                                                                   <  1:1067866/150‑1 (MQ=255)
                                                      aaCCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCTAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGAc                                                                                                 <  1:452367/150‑1 (MQ=255)
                                                       aCCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTACGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCGGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACt                                                                                                >  2:823241/1‑150 (MQ=255)
                                                         cTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGCCTGGACATCAATATGCTTAAAGTAAATAAGACTAt                                                                                              >  1:692418/1‑150 (MQ=255)
                                                         cTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGAGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTAt                                                                                              >  2:886249/1‑150 (MQ=255)
                                                           ggTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTc                                                                                            >  2:160490/1‑150 (MQ=255)
                                                            gTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTcc                                                                                           >  2:30972/1‑150 (MQ=255)
                                                                  cgcgCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACta                                                                                     >  2:466311/1‑150 (MQ=255)
                                                                    cgcATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACtatt                                                                                   >  1:1410328/1‑150 (MQ=255)
                                                                      cATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACtattat                                                                                 >  1:183320/1‑150 (MQ=255)
                                                                              gTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAAt                                                                         <  1:754642/150‑1 (MQ=255)
                                                                               tAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATg                                                                        <  1:498748/150‑1 (MQ=255)
                                                                                 aTCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCt                                                                      <  2:596561/150‑1 (MQ=255)
                                                                                 aTCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCt                                                                      <  2:963064/150‑1 (MQ=255)
                                                                                  tctcTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCtt                                                                     >  2:1317381/1‑150 (MQ=255)
                                                                                   ctctATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGCCGATTCCTGACTATTATTGATAAATGCttt                                                                    >  1:479373/1‑150 (MQ=255)
                                                                                   ctctATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATTAATGCttt                                                                    >  1:473165/1‑150 (MQ=255)
                                                                                   ctctATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCttt                                                                    <  1:671127/150‑1 (MQ=255)
                                                                                    tctATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCtttt                                                                   <  2:272081/150‑1 (MQ=255)
                                                                                     ctATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCACTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTa                                                                  >  2:518850/1‑150 (MQ=255)
                                                                                       aTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATTAATGCTTTTaaa                                                                >  1:1028871/1‑150 (MQ=255)
                                                                                       aTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCtcgtcac                                                                >  2:38939/1‑144 (MQ=255)
                                                                                        tCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTTGCGCGTAGAGTACGTGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAAc                                                               <  1:302369/150‑1 (MQ=255)
                                                                                         cgccgcGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAAcc                                                              >  2:1105764/1‑150 (MQ=255)
                                                                                         cgccgcGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAAcc                                                              >  1:110750/1‑150 (MQ=255)
                                                                                          gccgcGGCGTGCTCAGCCGTAGCTAAGCCGCGAAAAAAAATACGCGCCAGTTGAATTTCCATATGGCGCGGAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAAccc                                                             >  2:1356249/1‑150 (MQ=255)
                                                                                              cggcgTGCTCAGCAGTAGCTAAGCCGCGAACAAAAATACGCGCCAGTTGAATTTCCCTCTGGCGCGTAGAGTACGTGAGTGGACATCAATTTGCTTATAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAAcccgccc                                                         <  2:566839/150‑1 (MQ=255)
                                                                                                 cgTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGtt                                                      <  2:356230/150‑1 (MQ=255)
                                                                                                       cAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGttaattaa                                                >  2:1413203/1‑150 (MQ=255)
                                                                                                        aGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAAc                                               <  1:74572/150‑1 (MQ=255)
                                                                                                        aGCAGTAGCTAATCCGCGAACAAAAATACGCTCCAGGGGAATTTCCCTCTGGCGCGTAGAGTAAGTGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAAc                                               <  2:322255/150‑1 (MQ=255)
                                                                                                          ccGTAGCTAAGCCGCTAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTc                                             <  1:1694894/150‑1 (MQ=255)
                                                                                                           cGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACTGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCa                                            <  1:965060/150‑1 (MQ=255)
                                                                                                             tAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCAcc                                          <  1:1269806/150‑1 (MQ=255)
                                                                                                                  aaGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTg                                     <  1:1261004/150‑1 (MQ=255)
                                                                                                                    gCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGaa                                   >  1:1486756/1‑150 (MQ=255)
                                                                                                                         taaCAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGAGTAGAGTACTGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTca                              <  1:63162/149‑1 (MQ=255)
                                                                                                                                aaTACGCTCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACaataat                       <  2:1094799/150‑1 (MQ=255)
                                                                                                                                 aTACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAAtt                      >  1:679669/1‑150 (MQ=255)
                                                                                                                                  tACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTa                     <  1:1105764/150‑1 (MQ=255)
                                                                                                                                   aCGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCt                                                                                          >  1:289092/1‑80 (MQ=255)
                                                                                                                                   aCGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCt                                                                                          <  2:289092/80‑1 (MQ=255)
                                                                                                                                     gcgcCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGTGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGt                  <  2:509958/150‑1 (MQ=255)
                                                                                                                                     gcgcCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGt                  <  2:1105149/150‑1 (MQ=255)
                                                                                                                                      cgcCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTg                 >  1:1595427/1‑150 (MQ=255)
                                                                                                                                         cAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGata              <  2:1576451/150‑1 (MQ=255)
                                                                                                                                             gTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCTATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCg           >  2:1366042/1‑150 (MQ=255)
                                                                                                                                             gTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAAACCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCg           >  2:1350825/1‑150 (MQ=255)
                                                                                                                                               gAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCTAc         >  2:1514432/1‑150 (MQ=255)
                                                                                                                                               gAATTTCCCTCTGGCGAGTAGAGTACGGGAATGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTTACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGAc         <  1:1051516/150‑1 (MQ=255)
                                                                                                                                                    tCCCTCTGGCGCGTAGAGTACGGGTCTGGACTTCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAACCCCGCCCGGTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAgcgc    <  1:864057/150‑1 (MQ=255)
                                                                                                                                                      ccTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGt  >  2:118428/1‑150 (MQ=255)
                                                                                                                                                      ccTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGt  >  1:834109/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGT  >  NC_000913/88917‑89215

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