Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 11,127 C→G 10.6% S77T (AGT→ACT) 
L100V (CTA→GTA) 
yaaW ←
mbiA →
putative enzyme‑specific chaperone YaaW
uncharacterized protein MbiA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091311,1270CG10.6% 255.1 / 17.6 85S77T (AGT→ACT) 
L100V (CTA→GTA) 
yaaW
mbiA
putative enzyme‑specific chaperone YaaW
uncharacterized protein MbiA
Reads supporting (aligned to +/- strand):  ref base C (39/37);  new base G (3/6);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

CAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCA  >  NC_000913/10980‑11274
                                                                                                                                                   |                                                                                                                                                   
catgTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCTGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATAGTa                                                                                                                                                   <  2:1454046/147‑1 (MQ=255)
cAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTa                                                                                                                                                   >  1:1072532/1‑150 (MQ=255)
 agctATTACGCAGAAATTCTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGACTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTAt                                                                                                                                                  >  2:1098967/4‑150 (MQ=255)
  ggTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGCAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATc                                                                                                                                                 >  2:1002339/1‑150 (MQ=255)
      ttACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGcca                                                                                                                                             <  2:1186915/150‑1 (MQ=255)
       tACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGccac                                                                                                                                            >  2:408061/1‑150 (MQ=255)
        aCGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATAGTATCGccaac                                                                                                                                           <  2:930699/150‑3 (MQ=255)
         cGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGccacca                                                                                                                                          >  2:681546/1‑150 (MQ=255)
             gAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATAGTATCGCCAACaaact                                                                                                                                      <  1:277730/150‑3 (MQ=255)
             gAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAAt                                                                                                                                      <  2:1304701/150‑1 (MQ=255)
               aaTTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGGTGGCGATAGTATCGCCAACAAACTGc                                                                                                                                    <  1:245960/150‑1 (MQ=255)
                aTTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCATTTTGGTGGCGATAGTATCGCCAACAAACTGcg                                                                                                                                   <  2:1109730/150‑2 (MQ=39)
                  tGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCATTTTGGTGGCGATAGTATCGCCAACAAACTGcgcg                                                                                                                                 <  1:25436/150‑4 (MQ=39)
                    ttCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTaa                                                                                                                               >  1:54659/1‑150 (MQ=255)
                    ttCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTACGAGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTaa                                                                                                                               >  2:968005/1‑150 (MQ=255)
                       cAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATtc                                                                                                                            >  2:1392390/1‑150 (MQ=255)
                       cAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATtc                                                                                                                            >  2:1047253/1‑150 (MQ=255)
                         gTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATtctc                                                                                                                          <  1:747496/150‑1 (MQ=255)
                          tAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTcc                                                                                                                         <  2:483405/150‑1 (MQ=255)
                           aaCTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCg                                                                                                                        >  2:433184/1‑150 (MQ=255)
                           aaCTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCg                                                                                                                        >  2:677306/1‑150 (MQ=255)
                            aCTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCgg                                                                                                                       <  1:400016/150‑1 (MQ=255)
                             cTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGc                                                                                                                      <  1:9251/150‑1 (MQ=255)
                               gTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCaa                                                                                                                    >  2:9392/1‑150 (MQ=255)
                                 tGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATc                                                                                                                  >  2:697557/1‑150 (MQ=255)
                                   cTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAg                                                                                                                >  2:81913/1‑150 (MQ=255)
                                     cAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTACGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATtctc                                                                                                                          <  2:1469533/138‑1 (MQ=255)
                                     cAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTACGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATtctc                                                                                                                          >  1:1469533/1‑138 (MQ=255)
                                           cAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGt                                                                                                        <  2:1579443/150‑1 (MQ=255)
                                             aaCGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTg                                                                                                      >  1:1508649/1‑150 (MQ=255)
                                                   gACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGa                                                                                                >  1:1590817/1‑150 (MQ=255)
                                                     cATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATg                                                                                              <  2:397830/150‑1 (MQ=255)
                                                             tGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCgg                                                                                      >  2:768284/1‑150 (MQ=255)
                                                              gTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCTCAACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCggg                                                                                     >  1:287603/1‑150 (MQ=255)
                                                               tCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGAAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGa                                                                                    >  1:888843/1‑150 (MQ=255)
                                                               tCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGAAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGa                                                                                    >  1:883035/1‑150 (MQ=255)
                                                               tCGGCTTTCAGATTCAATCGCTTTGAAACATCGAGAAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATTCTGTAATTCTCCGGAAATAAGCTGACAGTTGCGGCGATGTTGCTCGGGa                                                                                    >  1:1359838/1‑150 (MQ=255)
                                                                cGGCTTTCAGCTTCAATCGCTTTGAAAGATCGAGCAAAATGGCCCGATACACTTTAGCGTGTCCGCGCTGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGAt                                                                                   <  1:187005/150‑1 (MQ=255)
                                                                cGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGAt                                                                                   >  2:873127/1‑150 (MQ=255)
                                                                  gCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCATTTTGGTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGc                                                                                 <  1:1427853/150‑1 (MQ=255)
                                                                   cTTTCAGCTTCAATAGCTTTGAAACATCGAGCAAAATGGACCTATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGcc                                                                                >  2:1089001/1‑150 (MQ=255)
                                                                    tttCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGccc                                                                               >  1:1169283/1‑150 (MQ=255)
                                                                       cAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTc                                                                            >  1:788549/1‑150 (MQ=255)
                                                                       cAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTc                                                                            >  2:1126385/1‑150 (MQ=255)
                                                                        aGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTcc                                                                           <  1:1558761/150‑1 (MQ=255)
                                                                            tCAATCGCTTTGAAACATCGAGGAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCg                                                                       >  1:1489963/1‑150 (MQ=255)
                                                                             cAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGa                                                                      <  2:1110945/150‑1 (MQ=255)
                                                                              aaTCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGAt                                                                     <  2:890747/150‑1 (MQ=255)
                                                                              aaTCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGAt                                                                     <  1:1392390/150‑1 (MQ=255)
                                                                               aTCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGAtt                                                                    <  2:788549/150‑1 (MQ=255)
                                                                               aTCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGAtt                                                                    <  2:794795/150‑1 (MQ=255)
                                                                                 cGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTa                                                                  <  2:1137605/150‑1 (MQ=255)
                                                                                 cGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTa                                                                  <  1:1433168/150‑1 (MQ=255)
                                                                                 cGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTa                                                                  <  1:677306/150‑1 (MQ=255)
                                                                                    tttGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAAc                                                               <  2:1413047/150‑1 (MQ=255)
                                                                                    tttGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAAc                                                               >  1:799816/1‑150 (MQ=255)
                                                                                     ttGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACa                                                              <  2:1072532/150‑1 (MQ=255)
                                                                                      tGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAg                                                             >  1:1389687/1‑150 (MQ=255)
                                                                                         aaCATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTc                                                          <  2:601505/150‑1 (MQ=255)
                                                                                          aCATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCTAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCg                                                         <  1:9392/150‑1 (MQ=255)
                                                                                          aCATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAACTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCg                                                         >  1:1121301/1‑150 (MQ=255)
                                                                                              cGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTgc                                                     >  2:887487/1‑150 (MQ=255)
                                                                                                   aaaaTGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATc                                                <  2:697600/150‑1 (MQ=255)
                                                                                                    aaaTGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCa                                               >  1:422292/1‑150 (MQ=255)
                                                                                                     aaTGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAg                                              <  1:1527772/150‑1 (MQ=255)
                                                                                                       tGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTa                                            >  2:229879/1‑150 (MQ=255)
                                                                                                          cccGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGc                                         >  2:93029/1‑150 (MQ=255)
                                                                                                             gATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGAAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTgg                                      <  2:883035/150‑1 (MQ=255)
                                                                                                             gATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGAAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTgg                                      <  2:888843/150‑1 (MQ=255)
                                                                                                              aTACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATAAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGga                                     >  2:906317/1‑150 (MQ=255)
                                                                                                                          gtgtCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGtttt                         <  1:768284/150‑1 (MQ=255)
                                                                                                                          gtgtCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGtttt                         >  1:58580/1‑150 (MQ=255)
                                                                                                                           tgtCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTg                        <  2:1169283/150‑1 (MQ=255)
                                                                                                                            gtCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGc                       <  1:1438642/150‑1 (MQ=255)
                                                                                                                               acagcATTTTGGTGGCGATAGTATCGCCAACAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGTCGATGTTGCTAGGGATGCCCTTCCATCGATTTAAAAAGTTCGTTGCGCATGAGTACGATGGAGAGGCGAGATTTGCCtt                    >  1:577589/4‑150 (MQ=255)
                                                                                                                                 cgcAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCtttt                  <  2:1590817/150‑1 (MQ=255)
                                                                                                                                  gcAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCttttt                 <  1:968005/150‑1 (MQ=255)
                                                                                                                                     gTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCAt              >  1:1511479/1‑150 (MQ=255)
                                                                                                                                       ttgttgGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTa            <  1:1304146/150‑1 (MQ=255)
                                                                                                                                       ttggtgGCGATAGTATCGCCAACAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGGGGCGAGTTTTGCCTTTTTCATTa            >  1:1028098/1‑150 (MQ=255)
                                                                                                                                        tggtgGCGATAGTATCGCCAACAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTAt           >  2:912118/4‑150 (MQ=255)
                                                                                                                                         gttgGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATg          <  1:229879/150‑1 (MQ=255)
                                                                                                                                             gCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTg      <  2:623883/150‑1 (MQ=255)
                                                                                                                                                aTACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGc   <  1:681546/150‑1 (MQ=255)
                                                                                                                                                 tACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCa  >  1:1254905/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCA  >  NC_000913/10980‑11274

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