Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,185,306 A→T 6.3% intergenic (+33/‑44) 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,3060AT 6.3% 270.5 / 5.6 76intergenic (+33/‑44)rpoB/rpoCRNA polymerase, beta subunit/RNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (35/36);  new base T (2/3);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

CGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCT  >  NC_000913/4185168‑4185451
                                                                                                                                          |                                                                                                                                                 
cGTGGACGGCAACCAACAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTAGTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGCCGGGTTTTAACCCgac                                                                                                                                          <  4:288719/148‑3 (MQ=255)
       ggCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGAt                                                                                                                                   <  4:925634/148‑1 (MQ=255)
        gCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGAtt                                                                                                                                  >  4:581774/1‑148 (MQ=255)
            ccATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTa                                                                                                                                                   >  2:94865/1‑127 (MQ=255)
            ccATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTa                                                                                                                                                   <  1:94865/127‑1 (MQ=255)
               tCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTaa                                                                                                                           <  4:642148/148‑1 (MQ=255)
                cAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAAc                                                                                                                          >  2:848704/1‑148 (MQ=255)
                  gATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTc                                                                                                                        >  4:90671/1‑148 (MQ=255)
                         ccGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACggg                                                                                                                 <  3:90671/148‑1 (MQ=255)
                          cGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGa                                                                                                                <  3:581774/148‑1 (MQ=255)
                            ggCCTGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAAc                                                                                                                                               >  4:917936/1‑115 (MQ=255)
                            ggCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAAc                                                                                                                                               <  3:917936/115‑1 (MQ=255)
                                tGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAAcc                                                                                                                                              >  1:675536/1‑112 (MQ=255)
                                tGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAAcc                                                                                                                                              <  2:675536/112‑1 (MQ=255)
                                   cAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCg                                                                                                       <  1:848704/148‑1 (MQ=255)
                                    aGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGt                                                                                                      >  4:281204/1‑148 (MQ=255)
                                     gAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATtgtg                                                                                                                               <  1:874830/122‑1 (MQ=255)
                                     gAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATtgtg                                                                                                                               >  2:874830/1‑122 (MQ=255)
                                         ccTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTaaa                                                                                                                                                 >  1:177116/1‑100 (MQ=255)
                                         ccTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTaaa                                                                                                                                                 <  2:177116/100‑1 (MQ=255)
                                             cAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGAttt                                                                                             <  3:281204/148‑1 (MQ=255)
                                                cGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATttatt                                                                                          >  3:7774/1‑148 (MQ=255)
                                                    ttgttgAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAg                                                                                      >  2:65590/1‑148 (MQ=255)
                                                        tgAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTc                                                                                  <  4:7774/148‑1 (MQ=255)
                                                           aaGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAg                                                                                                                                        <  4:223429/91‑1 (MQ=255)
                                                           aaGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAg                                                                                                                                        >  3:223429/1‑91 (MQ=255)
                                                           aaGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTaa                                                                                        >  4:638271/1‑139 (MQ=255)
                                                           aaGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTaa                                                                                        <  3:638271/139‑1 (MQ=255)
                                                            agagATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGCCGGGTTTTAACCCgaca                                                                                                                                         >  4:624248/1‑86 (MQ=255)
                                                            agagATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGCCGGGTTTTAACCCGACa                                                                                                                                         <  3:624248/89‑1 (MQ=255)
                                                                aTTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAgcg                                                                          <  1:65590/148‑1 (MQ=255)
                                                                           gggTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATtta                                                                                            <  2:89725/119‑1 (MQ=255)
                                                                           gggTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATtta                                                                                            >  1:89725/1‑119 (MQ=255)
                                                                                 cAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTaa                                                                                        <  3:45360/117‑1 (MQ=255)
                                                                                 cAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTaa                                                                                        >  4:45360/1‑117 (MQ=255)
                                                                                   aCATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTTTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGttt                                                                                   <  3:347108/120‑1 (MQ=255)
                                                                                   aCATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAg                                                                                                               >  3:841792/1‑92 (MQ=255)
                                                                                   aCATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAg                                                                                                               <  4:841792/92‑1 (MQ=255)
                                                                                   aCATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGttt                                                                                   >  4:347108/1‑120 (MQ=255)
                                                                                    cATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGAttt                                                                                             <  3:94577/109‑1 (MQ=255)
                                                                                    cATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGAttt                                                                                             >  4:94577/1‑109 (MQ=255)
                                                                                       cGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGaa                                                                                                   >  1:20131/1‑100 (MQ=255)
                                                                                       cGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGaa                                                                                                   <  2:20131/100‑1 (MQ=255)
                                                                                         aaCTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTg                                                    >  1:933902/1‑145 (MQ=255)
                                                                                         aaCTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTg                                                    <  2:933902/145‑1 (MQ=255)
                                                                                          aCTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGttt                                                                                   >  1:601283/1‑113 (MQ=255)
                                                                                          aCTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGttt                                                                                   <  2:601283/113‑1 (MQ=255)
                                                                                                     aGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGttt                                                                                   >  2:944611/1‑102 (MQ=255)
                                                                                                     aGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGttt                                                                                   <  3:91546/102‑1 (MQ=255)
                                                                                                     aGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGttt                                                                                   >  4:91546/1‑102 (MQ=255)
                                                                                                     aGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGttt                                                                                   <  1:944611/102‑1 (MQ=255)
                                                                                                      gTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCgg                                                                                                                                     <  3:1034592/51‑1 (MQ=255)
                                                                                                      gTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCgg                                                                                                                                     >  4:1034592/1‑51 (MQ=255)
                                                                                                      gTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATtta                                                                                            >  4:1129687/1‑92 (MQ=255)
                                                                                                      gTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATtta                                                                                            <  3:1129687/92‑1 (MQ=255)
                                                                                                      gTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGc                                    <  1:1001476/148‑1 (MQ=255)
                                                                                                      gTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGc                                    >  2:1001476/1‑148 (MQ=255)
                                                                                                       tAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATtta                                                                                            <  3:582104/91‑1 (MQ=255)
                                                                                                       tAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATtta                                                                                            >  4:582104/1‑91 (MQ=255)
                                                                                                                  cAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGAt                                                                                                                                   <  1:312105/41‑1 (MQ=255)
                                                                                                                  cAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGAt                                                                                                                                   >  2:312105/1‑41 (MQ=255)
                                                                                                                           cACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAagag                                                        <  4:512142/107‑1 (MQ=255)
                                                                                                                           cACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAagag                                                        >  3:512142/1‑107 (MQ=255)
                                                                                                                             ctgctgTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGttt                                                                                   <  4:119033/78‑1 (MQ=255)
                                                                                                                             ctgctgTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGttt                                                                                   >  3:119033/1‑78 (MQ=255)
                                                                                                                             ctgctgTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGAt                                                  >  3:862435/1‑111 (MQ=255)
                                                                                                                             ctgctgTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGAt                                                  <  4:862435/111‑1 (MQ=255)
                                                                                                                                 tgccGGGTTTTAACCCGACAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGa                                                            <  4:630017/94‑1 (MQ=255)
                                                                                                                                 tgccGGGTTTTAACCCGACAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGa                                                            >  3:630017/4‑97 (MQ=255)
                                                                                                                                   tCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGctct                                 <  3:569277/122‑1 (MQ=255)
                                                                                                                                   tCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGctct                                 >  4:569277/1‑122 (MQ=255)
                                                                                                                                      ggTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGa                                              <  4:302769/106‑1 (MQ=255)
                                                                                                                                      ggTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGa                                              <  4:479820/106‑1 (MQ=255)
                                                                                                                                      ggTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGa                                              >  3:479820/1‑106 (MQ=255)
                                                                                                                                      ggTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGa                                              >  3:302769/1‑106 (MQ=255)
                                                                                                                                        ttAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGa                >  4:459146/1‑134 (MQ=255)
                                                                                                                                        ttAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGa                <  3:459146/134‑1 (MQ=255)
                                                                                                                                        ttAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCt  >  3:916682/1‑148 (MQ=255)
                                                                                                                                         tAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGa                                                   <  4:693028/98‑1 (MQ=255)
                                                                                                                                         tAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGa                                                   >  3:693028/1‑98 (MQ=255)
                                                                                                                                          aaaaCCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAGAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATg                 >  3:611476/1‑131 (MQ=255)
                                                                                                                                          aaaaCCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAGAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATg                 <  4:611476/131‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
CGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCT  >  NC_000913/4185168‑4185451

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