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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,185,3070AT 6.3% 626.8 / 24.5 175intergenic (+34/‑43)rpoB/rpoCRNA polymerase, beta subunit/RNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (75/89);  new base T (4/7);  total (79/96)
Fisher's exact test for biased strand distribution p-value = 7.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.21e-01

Only 100 of 182 total aligned reads displayed.

TAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCG  >  NC_000913/4185159‑4185455
                                                                                                                                                    |                                                                                                                                                    
    aaCATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAAcc                                                                                                                                                  <  4:369668/149‑1 (MQ=255)
    aaCATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAAcc                                                                                                                                                  >  1:1419755/1‑149 (MQ=255)
     aCATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAAccc                                                                                                                                                 >  1:1209349/1‑149 (MQ=255)
     aCATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAAccc                                                                                                                                                 <  2:1419755/149‑1 (MQ=255)
      cATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCg                                                                                                                                                >  4:1725266/1‑149 (MQ=255)
      cATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCg                                                                                                                                                <  4:912965/149‑1 (MQ=255)
       aTCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCgg                                                                                                                                               >  2:1384996/1‑149 (MQ=255)
          gTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGCCGGGTTTTAACCCGACAg                                                                                                                                            <  1:648359/149‑1 (MQ=255)
           tGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGc                                                                                                                                           <  1:1840292/149‑1 (MQ=255)
             gACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCa                                                                                                                                             >  4:1066670/1‑145 (MQ=255)
              aCGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGa                                                                                                                                        <  1:727476/149‑1 (MQ=255)
                  cAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATtgt                                                                                                                                    <  3:498978/149‑1 (MQ=255)
                   aaCCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATtgtg                                                                                                                                   <  1:1159650/149‑1 (MQ=255)
                     ccATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGTTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCt                                                                                                                                 <  4:1079224/149‑1 (MQ=255)
                       aTCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTaa                                                                                                                               <  2:135274/149‑1 (MQ=255)
                       aTCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTaa                                                                                                                               <  4:1514823/149‑1 (MQ=255)
                       aTCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGCCGGGTTTTAACCCGACAGca                                                                                                                                          <  1:104651/138‑2 (MQ=255)
                       aTCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGCCGGGTTTTAACCCGACAGca                                                                                                                                          >  2:104651/1‑137 (MQ=255)
                        tCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAAc                                                                                                                              >  1:1287711/1‑149 (MQ=255)
                                 gCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACggg                                                                                                                     >  4:714231/1‑149 (MQ=255)
                                  ccGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGa                                                                                                                    <  3:1298464/149‑1 (MQ=255)
                                  ccGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGa                                                                                                                    <  3:791162/149‑1 (MQ=255)
                                  ccGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGa                                                                                                                    <  1:635862/149‑1 (MQ=255)
                                           ccAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGa                                                                                                                    >  1:789154/1‑140 (MQ=255)
                                            cAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAAccc                                                                                                                                                 >  3:1545349/1‑110 (MQ=255)
                                            cAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGt                                                                                                          >  3:1118156/1‑149 (MQ=255)
                                            cAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGt                                                                                                          <  1:919063/149‑1 (MQ=255)
                                            cAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGt                                                                                                          <  2:1287711/149‑1 (MQ=255)
                                            cAGAATCATTCAACGTGTTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAGccc                                                                                                                                                 <  4:1545349/110‑1 (MQ=255)
                                                aTCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTcc                                                                                                                           <  1:1287791/128‑1 (MQ=255)
                                                aTCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGa                                                                                                                    >  2:1626778/1‑135 (MQ=255)
                                                aTCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGa                                                                                                                    <  1:1626778/135‑1 (MQ=255)
                                                 tCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCt                                                                                                                                 >  3:772337/1‑121 (MQ=255)
                                                  ccTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGa                                                                                                    <  3:714231/149‑1 (MQ=255)
                                                    ttCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAAccc                                                                                                                                                 >  2:346123/1‑102 (MQ=255)
                                                    ttCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAAccc                                                                                                                                                 <  1:346123/102‑1 (MQ=255)
                                                    ttCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAAcc                                                                                                                                                  <  4:1386416/101‑1 (MQ=255)
                                                    ttCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGa                                                                                                                                        >  3:760064/1‑111 (MQ=255)
                                                    ttCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGa                                                                                                                    >  1:305694/1‑131 (MQ=255)
                                                    ttCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGa                                                                                                    <  1:1575293/147‑1 (MQ=255)
                                                    ttCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGa                                                                                                    <  4:773802/147‑1 (MQ=255)
                                                    ttCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGa                                                                                                    >  2:1575293/1‑147 (MQ=255)
                                                    ttCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGAtt                                                                                                  <  1:1106299/149‑1 (MQ=255)
                                                    ttCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTACCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAAcc                                                                                                                                                  >  3:1386416/1‑101 (MQ=255)
                                                       aaCGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGaaa                                                                                                      >  4:1143757/1‑142 (MQ=255)
                                                       aaCGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGaaa                                                                                                      <  3:1143757/142‑1 (MQ=255)
                                                        aCGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGt                                                                                                          >  1:151917/1‑137 (MQ=255)
                                                        aCGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGt                                                                                                          <  2:151917/137‑1 (MQ=255)
                                                           tATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGa                                                                                                                    >  1:1696252/1‑124 (MQ=255)
                                                           tATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTaaa                                                                                           <  1:1384996/149‑1 (MQ=255)
                                                           tATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTaaa                                                                                           <  2:64024/149‑1 (MQ=255)
                                                                 tgAAGGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCt                                                                                     <  2:1676684/149‑1 (MQ=255)
                                                                 tgAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCt                                                                                     <  3:1725266/149‑1 (MQ=255)
                                                                  gAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCa                                                                                                                                             <  3:573415/92‑1 (MQ=255)
                                                                   aaaGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAg                                                                                          <  4:42423/142‑1 (MQ=255)
                                                                   aaaGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGa                                                                                   >  3:237857/1‑149 (MQ=255)
                                                                   aaaGAGATTCGGTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGa                                                                                   >  3:1220934/1‑149 (MQ=255)
                                                                     agagATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGa                                                                                                                    <  2:631459/114‑1 (MQ=255)
                                                                     agagATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGaaa                                                                                 <  2:475168/149‑1 (MQ=255)
                                                                       agaATCGTTCGCTGGGTATCAACATCGTACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAgc                                                                               <  4:1220934/149‑1 (MQ=255)
                                                                        gaTTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAgcg                                                                              >  3:1052532/1‑149 (MQ=255)
                                                                         aTTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACg                                                                                                                       <  3:1536540/107‑1 (MQ=255)
                                                                          ttcgttcgCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTg                                                                                                         <  2:1841483/120‑1 (MQ=255)
                                                                           tcgttcgCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAg                                                                           >  1:529602/1‑149 (MQ=255)
                                                                            cgttcgCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGa                                                                          <  4:237857/149‑1 (MQ=255)
                                                                              ttcgCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACt                                                                        >  1:838140/1‑149 (MQ=255)
                                                                                      gTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGCCGGGTTTTAACCCGACAGca                                                                                                                                          <  2:64080/75‑2 (MQ=37)
                                                                                       tATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGaa                                                               <  1:902361/149‑1 (MQ=255)
                                                                                         tCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGt                                                                                                          <  2:74996/104‑1 (MQ=255)
                                                                                          cAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGt                                                                                                          <  1:184583/103‑1 (MQ=255)
                                                                                           aaCATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTc                                                                                      <  4:761607/122‑1 (MQ=255)
                                                                                           aaCATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCt                                                                                     <  4:543040/123‑1 (MQ=255)
                                                                                              aTCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGt                                                                                                          <  3:330752/99‑1 (MQ=255)
                                                                                                 gAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGa                                                       >  1:406050/1‑147 (MQ=255)
                                                                                                 gAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGa                                                       <  2:406050/147‑1 (MQ=255)
                                                                                                 gAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATg                                                     >  2:1534510/1‑149 (MQ=255)
                                                                                                    cTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGa                                                  <  1:650144/149‑1 (MQ=255)
                                                                                                       gtaGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAAc                                                                   >  4:1297278/3‑129 (MQ=255)
                                                                                                          gACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCaaaa                                            <  1:1534510/149‑1 (MQ=255)
                                                                                                          gACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCaaaa                                            <  2:1737077/149‑1 (MQ=255)
                                                                                                              aGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAGTCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGc                                        >  3:1286576/1‑149 (MQ=255)
                                                                                                               gTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGct                                       >  2:368589/1‑149 (MQ=255)
                                                                                                                tAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGctc                                      <  4:1286576/149‑1 (MQ=255)
                                                                                                                  aTTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTg                                    <  2:838140/149‑1 (MQ=255)
                                                                                                                   ttCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTgg                                   <  2:1209349/149‑1 (MQ=255)
                                                                                                                      tcGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGc                                  <  4:1370443/147‑1 (MQ=255)
                                                                                                                             aaCAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGa                         >  3:1580637/1‑149 (MQ=255)
                                                                                                                             aaCAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGa                         <  4:1580637/149‑1 (MQ=255)
                                                                                                                              aCAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGAtt                                                                                                  <  4:617666/75‑1 (MQ=255)
                                                                                                                                   tCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGAtt                                                                                                  <  3:446772/70‑1 (MQ=255)
                                                                                                                                     aCTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGc                             <  4:482427/137‑1 (MQ=255)
                                                                                                                                       tgctgTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGt               >  2:1123688/1‑149 (MQ=255)
                                                                                                                                        gctgTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGtt              <  1:442220/149‑1 (MQ=255)
                                                                                                                                           gTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCAt           >  2:441142/1‑149 (MQ=255)
                                                                                                                                              gggTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGt        <  1:1123688/149‑1 (MQ=255)
                                                                                                                                                  tAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCttt    <  2:277529/149‑1 (MQ=255)
                                                                                                                                                   aaaaCCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTc                               >  4:621820/1‑121 (MQ=255)
                                                                                                                                                   aaaaCCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTc                               >  3:991276/1‑121 (MQ=255)
                                                                                                                                                   aaaaCCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTc                               <  3:621820/121‑1 (MQ=255)
                                                                                                                                                    aaaCCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCg  >  3:1889701/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCG  >  NC_000913/4185159‑4185455

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