Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,185,307 A→T 5.6% intergenic (+34/‑43) rpoB → / → rpoC RNA polymerase subunit beta/RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,185,3070AT 5.6% 579.0 / 15.9 178intergenic (+34/‑43)rpoB/rpoCRNA polymerase subunit beta/RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base A (83/85);  new base T (2/8);  total (85/93)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 185 total aligned reads displayed.

ATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGG  >  NC_000913/4185158‑4185456
                                                                                                                                                     |                                                                                                                                                     
ataAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTaa                                                                                                                                                       <  2:124886/150‑1 (MQ=255)
ataAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTaa                                                                                                                                                       >  1:1163672/1‑150 (MQ=255)
 taAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTaaa                                                                                                                                                      >  1:345628/1‑150 (MQ=255)
      aCATCTTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCg                                                                                                                                                 <  2:119414/150‑1 (MQ=255)
      aCATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCg                                                                                                                                                 <  1:414/150‑1 (MQ=255)
        aTCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGc                                                                                                                                               <  2:294712/150‑1 (MQ=255)
        aTCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGc                                                                                                                                               <  1:1366519/150‑1 (MQ=255)
          cGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGCCGGGTTTTAACCCGACAg                                                                                                                                             >  1:1191564/1‑150 (MQ=255)
           gTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTTTTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGCCGGGTTTTAACCCGACAGc                                                                                                                                            <  1:1324944/150‑1 (MQ=255)
           gTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGc                                                                                                                                            >  2:138469/1‑150 (MQ=255)
            tGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCg                                                                                                                                           <  2:26186/150‑1 (MQ=255)
             ggACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGCCGGGTTTTAACCCGACAGcag                                                                                                                                          <  2:1089592/150‑3 (MQ=255)
              gACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGa                                                                                                                                         <  1:1379669/150‑1 (MQ=255)
               aCGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGCCGGGTTTTAACCCGACAGcagtg                                                                                                                                        <  1:404307/150‑5 (MQ=255)
                cGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGAtt                                                                                                                                       >  1:797660/1‑150 (MQ=255)
                  gCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATtgt                                                                                                                                     >  1:1095250/1‑150 (MQ=255)
                   cAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCTGCAGCGGATtttg                                                                                                                                    <  2:345628/150‑1 (MQ=255)
                   cAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATtgtg                                                                                                                                    <  2:1364487/150‑1 (MQ=255)
                    aaCCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGc                                                                                                                                   <  2:302878/150‑1 (MQ=255)
                       cATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTaa                                                                                                                                >  1:881905/1‑150 (MQ=255)
                        aTCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAAc                                                                                                                               <  2:1095250/150‑1 (MQ=255)
                           aGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTcc                                                                                                                            >  2:1208677/1‑150 (MQ=255)
                               ggAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACg                                                                                                                        <  1:15703/150‑1 (MQ=255)
                                gAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACgg                                                                                                                       >  1:1059153/1‑150 (MQ=255)
                                   ccGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAg                                                                                                                    <  1:904627/150‑1 (MQ=255)
                                   ccGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAg                                                                                                                    <  2:1111964/150‑1 (MQ=255)
                                   ccGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAg                                                                                                                    <  2:1157686/150‑1 (MQ=255)
                                    cgtgCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTGAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGc                                                                                                                   >  2:313443/4‑150 (MQ=255)
                                        cATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAAt                                                                                                               <  2:1303559/150‑1 (MQ=255)
                                        cATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACACGTCACTGCTGTCGGCTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAAt                                                                                                               >  2:899784/1‑150 (MQ=255)
                                           gCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCg                                                                                                            >  2:688527/1‑150 (MQ=255)
                                           gCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCg                                                                                                            >  1:177408/1‑150 (MQ=255)
                                                aaTCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGaaa                                                                                                       <  2:853350/150‑1 (MQ=255)
                                                   ccTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACAGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGAt                                                                                                    >  1:1351231/1‑150 (MQ=255)
                                                      tCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATtta                                                                                                 >  1:1037929/1‑150 (MQ=255)
                                                       cAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATttat                                                                                                >  1:818940/1‑150 (MQ=255)
                                                        aaCGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAATAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGt                                                                                                           <  2:1000836/138‑1 (MQ=255)
                                                        aaCGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAATAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGt                                                                                                           >  1:1000836/1‑138 (MQ=255)
                                                        aaCGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATttatt                                                                                               >  1:509724/1‑150 (MQ=255)
                                                         aCGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAg                                                                                                                                             <  1:442404/103‑1 (MQ=255)
                                                         aCGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATttatta                                                                                              >  1:858037/1‑150 (MQ=255)
                                                           gTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTaaa                                                                                            <  1:705784/150‑1 (MQ=255)
                                                           gTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTaaa                                                                                            <  2:438503/150‑1 (MQ=255)
                                                           gTAGTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTACTTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTaaa                                                                                            <  2:187353/150‑1 (MQ=255)
                                                            tATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAg                                                                                           >  1:216411/1‑150 (MQ=255)
                                                              ttgttgAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAGGATTTATTAAAGtt                                                                                         <  2:521981/150‑1 (MQ=255)
                                                               tgttgAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGttt                                                                                        >  1:391826/1‑150 (MQ=255)
                                                                gttgAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTc                                                                                       >  2:610680/1‑150 (MQ=255)
                                                                gttgAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTc                                                                                       >  1:634421/1‑150 (MQ=255)
                                                                  tgAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTg                                                                                     <  2:235778/150‑1 (MQ=255)
                                                                  tgAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTg                                                                                     <  2:700089/150‑1 (MQ=255)
                                                                   gAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGa                                                                                    >  1:227542/1‑150 (MQ=255)
                                                                    aaaGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGaa                                                                                   <  1:138469/150‑1 (MQ=255)
                                                                           ttcgttcgCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAg                                                                            <  2:124434/150‑1 (MQ=255)
                                                                             cgttcgCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGAc                                                                          <  2:1423327/150‑1 (MQ=255)
                                                                             cgttcgCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGAc                                                                          <  1:1295571/150‑1 (MQ=255)
                                                                                    tGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAAcc                                                                   <  2:1355198/150‑1 (MQ=255)
                                                                                    tGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAAcc                                                                   >  1:924507/1‑150 (MQ=255)
                                                                                    tGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAAcc                                                                   <  1:1208677/150‑1 (MQ=255)
                                                                                       gTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGaa                                                                >  1:753495/1‑150 (MQ=255)
                                                                                         aTCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAAAAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAaga                                                              >  2:1099794/1‑150 (MQ=255)
                                                                                          tCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCTGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAagag                                                             <  1:594196/150‑1 (MQ=255)
                                                                                           cAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAGACCCGGCAGTGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGt                                                            >  1:747977/1‑150 (MQ=255)
                                                                                            aaCATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGtt                                                           >  2:1257573/1‑150 (MQ=255)
                                                                                            aaCATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGtt                                                           >  1:1019321/1‑150 (MQ=255)
                                                                                                  gAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGc                                                     <  1:1247950/150‑1 (MQ=255)
                                                                                                     cTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGAt                                                  <  2:128851/150‑1 (MQ=255)
                                                                                                      tGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATc                                                 <  2:356233/150‑1 (MQ=255)
                                                                                                       cgAATACGAGTAATTCTCGCTCAAACAGGTCACTCCTCTCGGGTTAAAAACAGGCAGCGGATTGTGCTAACTCCGACGGGAGAAAATACGTCAAAGAATTATAAAAGTTTCTGAAAGCGCAGACTAAAACAGAAGAGTTTGATGAGAACa                                                >  2:1293820/2‑150 (MQ=255)
                                                                                                         aaGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCaaa                                              <  2:874836/150‑1 (MQ=255)
                                                                                                             cGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAAt                                            >  2:1179067/1‑148 (MQ=255)
                                                                                                             cGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAAt                                            <  1:1179067/148‑1 (MQ=255)
                                                                                                             cGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTg                                          >  2:200683/1‑150 (MQ=255)
                                                                                                                ttAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGctc                                       <  2:753495/149‑1 (MQ=255)
                                                                                                                gTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGctc                                       >  1:309815/1‑150 (MQ=255)
                                                                                                                  aaTTATAGCTCAAACAGGTCACTGCTGTCGGCTTAAAAACCGTCAGCGGATTGTGATAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTg                                     >  2:7744/1‑150 (MQ=255)
                                                                                                                      ctcGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCGGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCtt                                 >  2:408884/1‑150 (MQ=255)
                                                                                                                        cGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCg                               >  2:952667/1‑150 (MQ=255)
                                                                                                                        cGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCg                               <  2:645660/150‑1 (MQ=255)
                                                                                                                         gCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGc                              <  2:924507/150‑1 (MQ=255)
                                                                                                                           tCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCa                            <  1:1313758/150‑1 (MQ=255)
                                                                                                                            cAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAg                           <  2:1058034/150‑1 (MQ=255)
                                                                                                                             aaaCAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGa                          >  1:1282537/1‑150 (MQ=255)
                                                                                                                                cAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGGTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACAt                       <  1:969791/150‑1 (MQ=255)
                                                                                                                                 aGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATg                      <  2:1323225/150‑1 (MQ=255)
                                                                                                                                   gTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGAt                    >  1:572325/1‑150 (MQ=255)
                                                                                                                                      aCTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCg                 <  1:302910/150‑1 (MQ=255)
                                                                                                                                       ctgctgTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGt                >  1:1321192/1‑150 (MQ=255)
                                                                                                                                         gctgTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAGAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGCCATGATCCGTTc              >  2:1217947/1‑150 (MQ=255)
                                                                                                                                         gctgTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTc              >  1:122675/1‑150 (MQ=255)
                                                                                                                                         gctgTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTc              >  2:1070419/1‑150 (MQ=255)
                                                                                                                                         gctgCCGGGTTTTAACCCGACAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTc              >  2:1033411/1‑150 (MQ=255)
                                                                                                                                          ctgTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCa             <  1:663724/150‑1 (MQ=255)
                                                                                                                                               gggTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTc        <  2:487310/150‑1 (MQ=255)
                                                                                                                                                ggTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCg                               >  2:722319/1‑126 (MQ=255)
                                                                                                                                                ggTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCg                               <  1:722319/126‑1 (MQ=255)
                                                                                                                                                 gTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCtt      >  1:533884/1‑150 (MQ=255)
                                                                                                                                                   tAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCTTGTAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTc    <  1:96059/150‑1 (MQ=255)
                                                                                                                                                   tAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTc    <  2:775710/150‑1 (MQ=255)
                                                                                                                                                     aaaCCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCgg  >  1:267222/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGG  >  NC_000913/4185158‑4185456

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: