Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,676,195 A→G 5.2% intergenic (+87/+95) yhjE → / ← yhjG putative MFS transporter; membrane protein/putative inner membrane‑anchored periplasmic AsmA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,676,1950AG 5.2% 460.3 / 9.7 153intergenic (+87/+95)yhjE/yhjGputative MFS transporter; membrane protein/putative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base A (93/51);  new base G (4/4);  total (97/56)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

Only 100 of 155 total aligned reads displayed.

TGGCTGGCTTGACGTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAA  >  NC_000913/3676047‑3676344
                                                                                                                                                    |                                                                                                                                                     
tggctggcTTGACGTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCAATACAATTTATTGAATTAGCAAGGATTTGTAGGCCGGATaa                                                                                                                                                      >  2:464075/1‑150 (MQ=255)
 ggctggctTGACGTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAg                                                                                                                                                     >  2:781562/1‑150 (MQ=255)
 ggctggctTGACGTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAg                                                                                                                                                     >  2:433577/1‑150 (MQ=255)
  gctggctTGACGTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAgg                                                                                                                                                    >  2:1265168/1‑150 (MQ=255)
  gctggctTGACGTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAgg                                                                                                                                                    >  1:1483861/1‑150 (MQ=255)
  gctggctTGACGTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAgg                                                                                                                                                    >  1:1477184/1‑150 (MQ=255)
   ctggctTGACTTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGc                                                                                                                                                   <  1:149093/150‑1 (MQ=255)
     ggctTGACGTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGt                                                                                                                                                 <  2:1477184/150‑1 (MQ=255)
       gttGACGTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGGTc                                                                                                                                               <  2:208678/149‑1 (MQ=255)
        ttGACGGTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCa                                                                                                                                              <  2:85418/150‑1 (MQ=255)
          gACGTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCAcg                                                                                                                                            >  1:1295152/1‑150 (MQ=255)
          gACGTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCAcg                                                                                                                                            >  2:479685/1‑150 (MQ=255)
           aCGTTAATCGCCCTGCTGCTGAACCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCAcgc                                                                                                                                           <  1:1388955/150‑1 (MQ=255)
              ttAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAATGCGTTCAcgccgc                                                                                                                                        <  1:927343/150‑1 (MQ=255)
               tAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCa                                                                                                                                       >  1:910502/1‑150 (MQ=255)
                aaTCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                                                                      >  2:891558/1‑150 (MQ=255)
                  tCGCACTGCTGCTTACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                    <  2:242023/150‑1 (MQ=38)
                   cGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATAAATAGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTACACGCCTCATACg                                                                                                                                   >  2:298010/1‑150 (MQ=38)
                    gCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCCTACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                  >  1:1718762/1‑150 (MQ=255)
                      ccTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                                <  2:1329343/150‑1 (MQ=255)
                            gctgACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACa                                                                                                                          <  1:615528/150‑1 (MQ=255)
                            gctgACCCATGAGACGAGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACa                                                                                                                          >  1:986758/1‑150 (MQ=255)
                               gACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACaaac                                                                                                                       >  2:221242/1‑149 (MQ=255)
                                aCCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGc                                                                                                                      <  1:654225/150‑1 (MQ=255)
                                   cATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGATTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGcaca                                                                                                                   >  2:939721/1‑150 (MQ=255)
                                    aTGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGAATCCGGCATGAAGAAAGCacac                                                                                                                  >  2:68834/1‑150 (MQ=255)
                                        gaCGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCATGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTg                                                                                                              >  2:1566915/1‑150 (MQ=255)
                                             gACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCcaac                                                                                                         >  1:1697822/1‑150 (MQ=255)
                                               cATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCcaacaa                                                                                                       >  2:534517/1‑150 (MQ=255)
                                               cATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCcaacaa                                                                                                       >  1:482777/1‑150 (MQ=255)
                                               cATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCcaacaa                                                                                                       >  2:15925/1‑150 (MQ=255)
                                                aTCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAAt                                                                                                      <  2:1469245/150‑1 (MQ=255)
                                                   aGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTa                                                                                                   <  1:829039/150‑1 (MQ=255)
                                                   aGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTa                                                                                                   >  1:1325713/1‑150 (MQ=255)
                                                   aGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTa                                                                                                   >  1:606831/1‑150 (MQ=255)
                                                       gTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATc                                                                                               >  2:1500246/1‑150 (MQ=255)
                                                         tGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCtt                                                                                             >  2:372323/1‑149 (MQ=255)
                                                          gTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGt                                                                                            >  1:1375911/1‑150 (MQ=255)
                                                           taataaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGTACACTTTGCCAACAATCTAAAGCTGtt                                                                                           >  1:1450296/1‑150 (MQ=255)
                                                           taataaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGtt                                                                                           >  2:1144630/1‑150 (MQ=255)
                                                           taataaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGtt                                                                                           >  1:360458/1‑150 (MQ=255)
                                                           taataaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGtt                                                                                           <  2:884626/150‑1 (MQ=255)
                                                             ataaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTa                                                                                         >  1:848878/1‑150 (MQ=255)
                                                             ataaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTa                                                                                         <  1:15925/150‑1 (MQ=255)
                                                             aaaaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGAAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCAAACTTTGCCAACAATCTAATACTGTTTa                                                                                         >  1:873078/3‑150 (MQ=255)
                                                               aaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAAt                                                                                       >  2:1408281/1‑150 (MQ=255)
                                                               aaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAAt                                                                                       >  2:1081001/1‑150 (MQ=255)
                                                                agagTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAAtt                                                                                      >  2:1699124/1‑150 (MQ=255)
                                                                   gTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTCAAACTGTTTAATTgac                                                                                   >  2:292093/1‑148 (MQ=255)
                                                                           ccTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                    <  2:1525016/93‑1 (MQ=25)
                                                                           ccTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                    >  1:1525016/1‑93 (MQ=25)
                                                                             tGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGa                                                                         <  1:781562/150‑1 (MQ=255)
                                                                              gATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGAc                                                                        >  1:1261289/1‑150 (MQ=255)
                                                                               aTGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACa                                                                       >  1:1668039/1‑150 (MQ=255)
                                                                                   gcTACGCTTATCAGGCGTTAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATTAGGCGTGCACGCCGGATCCGGCATGAACTAAGCAGACTTTGCCAACAATCTAATTCTGTTTAATTGCCGGATGTCAGACATCCg                                                                   <  1:56382/150‑1 (MQ=255)
                                                                                   gcTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTAAAAAAACACACATTTCCAACAATAAAAATATTTTTAATTTACCGATGTATGACAtacc                                                                   >  1:1578400/1‑147 (MQ=18)
                                                                                    cTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCgg                                                                  <  1:1518511/150‑1 (MQ=255)
                                                                                     tACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAATTATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCTGACATCCGGc                                                                 <  1:493084/150‑1 (MQ=255)
                                                                                       cGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCaa                                                               >  1:1121085/1‑150 (MQ=255)
                                                                                        gCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCTGACATCCGGCAAt                                                              <  1:496255/150‑1 (MQ=255)
                                                                                        gCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCTGACATCCGGCAAt                                                              <  1:1245589/150‑1 (MQ=255)
                                                                                         cTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAAtt                                                             >  2:982842/1‑150 (MQ=255)
                                                                                             tCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCTGACATCCGGCAATTAaat                                                         >  1:534520‑M2/1‑147 (MQ=255)
                                                                                             tCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCTGACATCCGGCAATTAaac                                                         >  1:534517‑M2/1‑147 (MQ=255)
                                                                                             tCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTatca                                                         <  1:1534740/150‑1 (MQ=255)
                                                                                             tCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTatca                                                         >  1:919837/1‑150 (MQ=255)
                                                                                               aGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTatcatc                                                       <  1:1488578/150‑1 (MQ=255)
                                                                                                ggCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTAtcatca                                                      >  2:580698/1‑150 (MQ=255)
                                                                                                 gCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTatcatccc                                                     >  2:938014/1‑148 (MQ=255)
                                                                                                 gCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCAc                                                     >  2:373803/1‑150 (MQ=255)
                                                                                                            cATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTg                                          >  1:782993/1‑150 (MQ=255)
                                                                                                            cATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTg                                          <  1:482298/150‑1 (MQ=255)
                                                                                                             aTACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGc                                         >  2:1163032/1‑150 (MQ=255)
                                                                                                               aCAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGa                                       >  2:329037/1‑150 (MQ=255)
                                                                                                                     tatTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACaaaaa                                 >  2:843251/1‑149 (MQ=255)
                                                                                                                     tatTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAAt                                 >  1:167959/1‑150 (MQ=255)
                                                                                                                     tatTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAAt                                 >  1:185101/1‑150 (MQ=255)
                                                                                                                       ttGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGTATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCg                               <  2:1350402/150‑1 (MQ=255)
                                                                                                                       ttGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCg                               <  2:484200/150‑1 (MQ=255)
                                                                                                                           aTTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCg                           <  1:1265168/150‑1 (MQ=255)
                                                                                                                             tCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGac                         <  2:1529356/150‑1 (MQ=255)
                                                                                                                             tCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCAGGCAATTATCATCACTTCTTCATCTGCGACTAAATCGTCCGac                         >  2:344900/1‑150 (MQ=255)
                                                                                                                              cGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGaca                        >  2:1620131/1‑150 (MQ=255)
                                                                                                                               gCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGacac                       >  1:1138797/1‑150 (MQ=255)
                                                                                                                                 aaGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTg                     >  2:1496099/1‑150 (MQ=255)
                                                                                                                                    gATTTGTAGGACGGATAAGGCGTTCACGCCGCATACGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGAtt                  <  2:1286942/150‑1 (MQ=255)
                                                                                                                                        ttcatgCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGcc              >  2:556700/6‑150 (MQ=38)
                                                                                                                                        tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGcc              <  1:888061/150‑1 (MQ=255)
                                                                                                                                        tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGcc              <  1:1052486/150‑1 (MQ=255)
                                                                                                                                         gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCt             >  1:1225013/1‑150 (MQ=255)
                                                                                                                                          tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAATGCACACTTTGACAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTc            <  1:298010/150‑1 (MQ=255)
                                                                                                                                           aGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGGCAAAATCGTCCGACACTGATTAGCCTcc           <  1:1144630/150‑1 (MQ=255)
                                                                                                                                           aGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTcc           <  2:167959/150‑1 (MQ=255)
                                                                                                                                           aGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATCTCAGACATCCGGCATTTATCCTCAGTTCTTAATCTTCGACAAAATCGTCCGACACTGAGTAGACTcc           >  2:334481/1‑150 (MQ=37)
                                                                                                                                            tgCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTccc          >  1:751483/2‑150 (MQ=39)
                                                                                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTccc          >  2:425291/1‑150 (MQ=255)
                                                                                                                                              ccGGATAAGGCGTTAACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACTCTGATTAGCCTcccca        >  1:455776/1‑149 (MQ=255)
                                                                                                                                                 gATATGGCGTTCACTCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCg     <  1:221242/150‑1 (MQ=255)
                                                                                                                                                 gATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCg     <  2:848878/150‑1 (MQ=255)
                                                                                                                                                    aaGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGaa  >  2:1571893/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGGCTGGCTTGACGTTAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAA  >  NC_000913/3676047‑3676344

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