Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 248,868 C→G 11.8% pseudogene (1203/1713 nt) lfhA ← flagellar system protein, promoterless fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913248,8680CG11.8% 242.9 / 21.7 85pseudogene (1203/1713 nt)lfhAflagellar system protein, promoterless fragment
Reads supporting (aligned to +/- strand):  ref base C (38/35);  new base G (6/4);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

CACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTGTCGATCACCTG  >  NC_000913/248724‑249017
                                                                                                                                                |                                                                                                                                                     
cacaTCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCt                                                                                                                                                  >  2:877771/1‑150 (MQ=255)
  caTCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCGGCGCTCa                                                                                                                                                >  1:461716/1‑150 (MQ=255)
   aTCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCaa                                                                                                                                               >  1:497412/1‑150 (MQ=255)
     cGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCGGCGCTCaatt                                                                                                                                             <  2:461716/150‑3 (MQ=255)
     cGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCGGCGCTCaatt                                                                                                                                             <  2:718906/150‑3 (MQ=255)
        ccgccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTCCACTTTCAGCAACTGGCTGTAATTGAGCGCGGCGCTCaattaca                                                                                                                                          <  1:977745/150‑6 (MQ=255)
               cagAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCAGCGCTCAAATCTTCCGCCAg                                                                                                                                   >  2:1483881/1‑150 (MQ=255)
                   atatGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGAcgc                                                                                                                               >  1:1078096/1‑150 (MQ=255)
                     atGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCgg                                                                                                                             <  2:686100/150‑1 (MQ=255)
                     atGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACCCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCgg                                                                                                                             >  2:1269945/1‑150 (MQ=255)
                      tGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGt                                                                                                                            >  2:1381146/1‑150 (MQ=255)
                      tGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGt                                                                                                                            >  1:263456/1‑150 (MQ=255)
                      tGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGt                                                                                                                            >  2:314165/1‑150 (MQ=255)
                      tGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGt                                                                                                                            <  1:957305/150‑1 (MQ=255)
                      tGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGt                                                                                                                            >  1:1095711/1‑150 (MQ=255)
                       gATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTg                                                                                                                           >  2:95110/1‑150 (MQ=255)
                        aTCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGc                                                                                                                          <  2:397714/150‑1 (MQ=255)
                        aTCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACGTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGc                                                                                                                          >  2:631083/1‑150 (MQ=255)
                          cTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCg                                                                                                                        <  2:1078096/150‑1 (MQ=255)
                          cTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCg                                                                                                                        <  2:1120879/150‑1 (MQ=255)
                          cTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCg                                                                                                                        <  2:1115492/150‑1 (MQ=255)
                           catCGCCACCGTGCTGGTCGCCAGTAGCACGGTGACGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGt                                                                                                                       >  1:1007595/3‑150 (MQ=38)
                            ttCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTc                                                                                                                      <  1:415420/150‑1 (MQ=255)
                              cGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGACTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcga                                                                                                                    >  2:1373977/1‑150 (MQ=255)
                               gTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgac                                                                                                                   >  1:983156/1‑150 (MQ=255)
                               gTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgac                                                                                                                   <  1:23529/150‑1 (MQ=255)
                                   ccGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCCGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATa                                                                                                               >  2:1151767/1‑150 (MQ=255)
                                    cGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATaa                                                                                                              <  1:435942/150‑1 (MQ=255)
                                        cTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACgg                                                                                                          >  1:568390/1‑150 (MQ=255)
                                         tACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGAAGGGAAACGCCTTCGGTCAGCAGCGCCCGTTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGt                                                                                                         <  1:1354273/150‑1 (MQ=255)
                                          aCTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGtt                                                                                                        <  2:478273/150‑1 (MQ=255)
                                                     gCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCg                                                                                             >  1:127988/1‑150 (MQ=255)
                                                     gCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCAGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCg                                                                                             >  2:1487113/1‑150 (MQ=255)
                                                       aCGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTa                                                                                           <  2:418278/150‑1 (MQ=255)
                                                         ggTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGAAAAACGGTTATGCAACTGCGTAat                                                                                         <  2:972287/150‑1 (MQ=255)
                                                           tGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAatat                                                                                       <  2:788961/150‑1 (MQ=255)
                                                               gATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATatcatc                                                                                   <  1:298464/150‑1 (MQ=255)
                                                                 tGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATcatcat                                                                                 <  1:1483881/150‑1 (MQ=255)
                                                                  ggTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATa                                                                                >  1:848302/1‑150 (MQ=255)
                                                                    tGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGt                                                                              >  2:804956/1‑150 (MQ=255)
                                                                    tGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGt                                                                              >  1:11896/1‑150 (MQ=255)
                                                                      aCGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGGCGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTa                                                                            >  1:1020843/1‑150 (MQ=255)
                                                                       cGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTaa                                                                           <  1:877771/150‑1 (MQ=255)
                                                                       cGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTaa                                                                           <  1:744388/150‑1 (MQ=255)
                                                                         atatCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAAc                                                                         <  2:160412/150‑1 (MQ=255)
                                                                             cgcgCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAAt                                                                     >  2:870679/1‑150 (MQ=255)
                                                                             cgcgCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAAt                                                                     >  2:1418978/1‑150 (MQ=255)
                                                                                    ggAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAat                                                              >  2:510778/1‑150 (MQ=255)
                                                                                    ggAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAat                                                              <  2:378740/150‑1 (MQ=255)
                                                                                     gAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAata                                                             >  2:528449/1‑150 (MQ=255)
                                                                                         cGCCTTCGGTCAGCAGCGCCCGGTACACCTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGc                                                         >  2:1271455/1‑150 (MQ=255)
                                                                                           ccTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTg                                                       >  2:508990/1‑150 (MQ=255)
                                                                                             ttCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTgcg                                                     <  1:1202428/150‑1 (MQ=255)
                                                                                                ggTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGcac                                                  >  1:31171/1‑150 (MQ=255)
                                                                                                  tCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACaa                                                >  2:817879/1‑150 (MQ=255)
                                                                                                     gcagcGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCt                                             <  1:778433/150‑1 (MQ=255)
                                                                                                     gcagcGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCt                                             >  1:76904/1‑150 (MQ=255)
                                                                                                       agcGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTa                                           <  2:1007595/150‑1 (MQ=255)
                                                                                                         cgcCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTAtt                                         <  1:528449/150‑1 (MQ=255)
                                                                                                         cgcCCGCTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGGGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCAAAGTTAAACAAATCAGGAATATAGCTGCGCAGAATCtgagt                                         >  2:680914/1‑146 (MQ=255)
                                                                                                          gcCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGGCGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTAttt                                        <  2:1020843/150‑1 (MQ=255)
                                                                                                          gcCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTAttt                                        <  1:1093579/150‑1 (MQ=255)
                                                                                                           cccGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTa                                       >  2:1215661/1‑150 (MQ=255)
                                                                                                           cccGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTa                                       <  2:127988/150‑1 (MQ=255)
                                                                                                               gTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCTACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATg                                   <  1:1151767/150‑1 (MQ=255)
                                                                                                                  cacTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGt                                >  2:998569/1‑150 (MQ=255)
                                                                                                                  cacTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGt                                <  2:760348/150‑1 (MQ=255)
                                                                                                                    cTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAg                              >  1:420168/1‑150 (MQ=255)
                                                                                                                     tttCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGc                             <  1:1418978/150‑1 (MQ=255)
                                                                                                                      ttCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCg                            <  2:694850/150‑1 (MQ=255)
                                                                                                                       tCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGa                           <  1:207906/150‑1 (MQ=255)
                                                                                                                       tCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGa                           <  1:95110/150‑1 (MQ=255)
                                                                                                                                ggCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTgg                  >  2:525262/1‑150 (MQ=255)
                                                                                                                                ggCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTgg                  >  1:1368883/1‑150 (MQ=255)
                                                                                                                                 gCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGgc                 >  1:760616/1‑150 (MQ=255)
                                                                                                                                     tAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTg             <  2:568390/150‑1 (MQ=255)
                                                                                                                                     gatttGAGCGCGGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTg             >  2:92033/4‑150 (MQ=255)
                                                                                                                                     gatttGAGCGCGGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTg             >  2:1005299/4‑150 (MQ=255)
                                                                                                                                      atttGAGCGCGGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTGt            <  2:1283774/148‑1 (MQ=255)
                                                                                                                                      atttGAGCGCGGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTGt            >  1:807109/3‑150 (MQ=255)
                                                                                                                                      atttGAGCGCGGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTGt            >  2:573768/3‑150 (MQ=255)
                                                                                                                                         tGAGCGCGGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTGTCGa         >  1:1523239/1‑150 (MQ=255)
                                                                                                                                          gAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTGTCGAt        >  1:846866/1‑150 (MQ=255)
                                                                                                                                           aGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTGTCGATc       <  2:1368883/150‑1 (MQ=255)
                                                                                                                                            gcgcCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTGTCGATCa      >  2:1512592/1‑150 (MQ=255)
                                                                                                                                               ccgcGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTGTCGATCACCt   >  2:905946/1‑150 (MQ=255)
                                                                                                                                                cgcgCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTGTCGATCACCTg  >  1:1174488/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTGTCGATCACCTG  >  NC_000913/248724‑249017

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: