Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,601,505 A→C 5.3% R194R (CGT→CGG ftsX ← cell division protein FtsX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,601,5050AC 5.3% 236.1 / 5.1 75R194R (CGT→CGGftsXcell division protein FtsX
Reads supporting (aligned to +/- strand):  ref base A (36/35);  new base C (2/2);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

ACACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACG  >  NC_000913/3601356‑3601654
                                                                                                                                                     |                                                                                                                                                     
acacTGTTACCTATGACGAGGAACTCGGCCTCCACCATCTACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATa                                                                                                                                                       <  2:504507/150‑1 (MQ=255)
acacTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATa                                                                                                                                                       <  2:1150390/150‑1 (MQ=255)
     gTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATc                                                                                                                                                  <  2:167893/150‑1 (MQ=255)
         ccGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATCGTATCACGc                                                                                                                                              <  2:51596/150‑1 (MQ=255)
          cGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCTCTTCGTCAATGCGGTTAATCTTCGTGATCTTATCACGCa                                                                                                                                             <  2:493061/150‑1 (MQ=255)
             tGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCg                                                                                                                                          <  1:1113848/150‑1 (MQ=255)
             tGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCg                                                                                                                                          <  1:1611980/150‑1 (MQ=255)
                 gAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTAtt                                                                                                                                      >  2:483367/1‑150 (MQ=255)
                  aGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTc                                                                                                                                     <  2:341136/150‑1 (MQ=255)
                     aaCACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGTTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCAGTTCGTCAATGCTGTTAATCTGCGTGATAGGATCACGGAGCGTATTCAGt                                                                                                                                  <  2:568492/150‑1 (MQ=255)
                      acacGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTg                                                                                                                                 <  1:1023481/150‑1 (MQ=255)
                      acacGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTg                                                                                                                                 <  1:966167/150‑1 (MQ=255)
                          gtcctccaacaTCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTc                                                                                                                             <  2:592417/145‑1 (MQ=255)
                           gccgccACCATCAACACGCCGATCATCGCCGAACCGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTcc                                                                                                                            >  2:225651/1‑150 (MQ=255)
                             cgccATCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGt                                                                                                                          >  1:586494/1‑150 (MQ=255)
                               ccaccaTCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTcc                                                                                                                        >  2:1646985/1‑150 (MQ=255)
                               ccaccaTCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAAGCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTcc                                                                                                                        <  2:95100/150‑1 (MQ=255)
                                  ccaTCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCt                                                                                                                     <  1:501545/150‑1 (MQ=255)
                                    aTCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTgg                                                                                                                   <  2:586494/150‑1 (MQ=255)
                                     tcgacacGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGa                                                                                                                  <  1:50688/147‑1 (MQ=255)
                                     tCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGa                                                                                                                  <  1:1226599/150‑1 (MQ=255)
                                           cGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGa                                                                                                            <  1:967432/150‑1 (MQ=255)
                                                   atcGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCggg                                                                                                    >  2:1042668/1‑150 (MQ=255)
                                                   atcGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCggg                                                                                                    >  1:1037214/1‑150 (MQ=255)
                                                    tcGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGa                                                                                                   <  1:942364/150‑1 (MQ=255)
                                                     ggCCGAAACGCGCCGGAGCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTGCGTCCCCTGGAAATCGAGTTTCGGGAt                                                                                                  <  2:529654/149‑1 (MQ=255)
                                                       ccGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCTGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATca                                                                                                <  2:957570/150‑1 (MQ=255)
                                                       ccGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATca                                                                                                >  2:134481/1‑150 (MQ=255)
                                                             cgcgCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCg                                                                                          >  2:1618226/1‑150 (MQ=255)
                                                               cgcCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGcc                                                                                        >  1:1459627/1‑150 (MQ=255)
                                                                gcCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAACCCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCCCCACCGCCa                                                                                       >  1:391576/1‑150 (MQ=255)
                                                                gcCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCggg                                                                                                    >  1:424823/1‑137 (MQ=255)
                                                                gcCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCggg                                                                                                    <  2:424823/137‑1 (MQ=255)
                                                                   cGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTg                                                                                    <  1:483367/150‑1 (MQ=255)
                                                                      ccAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCg                                                                                 <  2:1061711/150‑1 (MQ=255)
                                                                       cAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCgg                                                                                >  2:306919/1‑150 (MQ=255)
                                                                        aGCCCGGTCAACGCCGCCAGACGGGCAAGCCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAAACTGCGTGATACGATCACGCAGCg                                                                                                                                          <  1:860416/91‑1 (MQ=255)
                                                                        aGCCCGGTCAACGCCGCCAGACGGGCAAGCCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAAACTGCGTGATACGATCACGCAGCg                                                                                                                                          >  2:860416/1‑91 (MQ=255)
                                                                        aGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAACTCGAGTTTCGGGATCCCCACCGCCACTGCCGGa                                                                               >  2:638329/1‑150 (MQ=255)
                                                                         gCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGaa                                                                              >  2:686674/1‑150 (MQ=255)
                                                                              gTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCggg                                                                         <  1:79520/150‑1 (MQ=255)
                                                                              gTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCggg                                                                         <  1:1646985/150‑1 (MQ=255)
                                                                               tCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGt                                                                        >  2:267528/1‑150 (MQ=255)
                                                                                 aaCGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGttt                                                                      <  1:978982/150‑1 (MQ=255)
                                                                                 aaCGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGACGGAAGCTGGttt                                                                      >  2:704643/1‑150 (MQ=255)
                                                                                  aCGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGtttt                                                                     <  1:686674/150‑1 (MQ=255)
                                                                                   cgccgcCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCg                                                                                          <  1:828691/128‑1 (MQ=255)
                                                                                   cgccgcCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTttc                                                                    >  1:1531010/1‑150 (MQ=255)
                                                                                   cgccgcCAGACGGGCAAACCAGCTGTCATCCATCAGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCg                                                                                          >  2:828691/1‑128 (MQ=255)
                                                                                     ccgccAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGt                                                                                                          <  1:515710/110‑1 (MQ=255)
                                                                                     ccgccAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGt                                                                                                          >  2:515710/1‑110 (MQ=255)
                                                                                           gACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCa                                                            >  1:1080472/1‑150 (MQ=255)
                                                                                           caCGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCa                                                            <  1:1431275/149‑1 (MQ=255)
                                                                                             cGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCata                                                          >  2:294501/1‑150 (MQ=255)
                                                                                              gggCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCatat                                                         >  1:1633164/1‑150 (MQ=255)
                                                                                                gCAAACCATCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATcc                                                       <  2:1411966/150‑1 (MQ=255)
                                                                                                         cTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGcacca                                              >  1:1224837/1‑150 (MQ=255)
                                                                                                           gTCATCCATCCGCACTTCGTCAATGACGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGTGCaccacc                                            >  2:1279726/1‑150 (MQ=255)
                                                                                                             catccatcCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCaa                                          >  2:598467/1‑150 (MQ=255)
                                                                                                               tccatccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCaaaa                                        >  1:215034/1‑150 (MQ=255)
                                                                                                                 catccGCACTTCGTCAATGCCGTTAATCTGCGTGGTACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAacc                                      >  1:263862/1‑150 (MQ=255)
                                                                                                                   tccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAaccag                                    <  1:306919/150‑1 (MQ=255)
                                                                                                                   tccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGGGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCACCATATCCAGCGCACCACCAAAaccag                                    <  1:1353940/150‑1 (MQ=255)
                                                                                                                      gCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAaccagacc                                 >  2:401553/1‑150 (MQ=255)
                                                                                                                        aCTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAaccagaccag                               <  2:1337183/150‑1 (MQ=255)
                                                                                                                           tCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTa                            >  1:864208/1‑150 (MQ=255)
                                                                                                                              tCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACgg                         >  2:1196576/1‑150 (MQ=255)
                                                                                                                                   gCCGTTAATCTGCGTGATACGATCACGCAGCGTATTAAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTc                    >  2:1383373/1‑150 (MQ=255)
                                                                                                                                        tAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCa               >  2:1291297/1‑150 (MQ=255)
                                                                                                                                          cgcTGCGTGATCGTATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCCCCACCAAAACCAGACCAGTTACGGAACTCACCCAGt             >  1:1151621/3‑150 (MQ=255)
                                                                                                                                           tCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTg            <  2:215034/150‑1 (MQ=255)
                                                                                                                                           tCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTg            <  1:386293/150‑1 (MQ=255)
                                                                                                                                            cTGCGTGATCGTATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGc           >  2:559944/1‑150 (MQ=255)
                                                                                                                                            cTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGc           >  1:838107/1‑150 (MQ=255)
                                                                                                                                             tGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCg          <  1:339924/150‑1 (MQ=255)
                                                                                                                                                 tGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCtt      >  1:243864/1‑150 (MQ=255)
                                                                                                                                                  gATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTc     >  1:1301276/1‑150 (MQ=255)
                                                                                                                                                     aCGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACg  >  2:108477/1‑150 (MQ=255)
                                                                                                                                                     aCGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACg  <  2:1037214/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACG  >  NC_000913/3601356‑3601654

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: