Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,187,989 A→G 9.4% pseudogene (563/2439 nt) yqiG → pseudogene; fimbrial export usher family;putative membrane; Not classified; putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,187,9890AG 9.4% 308.2 / 20.3 117pseudogene (563/2439 nt)yqiGpseudogene; fimbrial export usher family;putative membrane; Not classified; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base A (55/51);  new base G (6/5);  total (61/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.25e-01

Only 100 of 117 total aligned reads displayed.

GAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTC  >  NC_000913/3187841‑3188136
                                                                                                                                                    |                                                                                                                                                   
gAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCATATGGGCAGGCTGGCGTTAATCTTGCGGCCTGGCGCCTGCGTACTGAt                                                                                                                                                    >  2:325526/1‑150 (MQ=255)
  aCTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTa                                                                                                                                                  >  2:326045/1‑150 (MQ=255)
   cTTTGGAATCATGGTATTGCTGGTGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCATAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTAc                                                                                                                                                 <  2:1235248/150‑1 (MQ=255)
      tGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAg                                                                                                                                              <  2:359621/150‑1 (MQ=255)
        gAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGGTAATCAGTa                                                                                                                                            <  1:181662/150‑1 (MQ=255)
          aTCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACg                                                                                                                                          <  2:927357/150‑1 (MQ=255)
             aTGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATc                                                                                                                                       <  2:900903/150‑1 (MQ=255)
             aTGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATc                                                                                                                                       >  2:1540535/1‑150 (MQ=255)
               ggTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAg                                                                                                                                     <  2:387999/150‑1 (MQ=255)
                gTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTAGGCACCACATCAGCGCGATAATTCGCAAAATATAACTTCCTATGGGCAGGCTGGGGTTATTCATGGGGCCTGGCGCGTGCGAACTGATTACCAGTACCATCact                                                                                                                                    >  1:316318/1‑148 (MQ=255)
                 tATTGCTGGGGCATTTCTTGACTACAATTTATTTGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGTTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTc                                                                                                                                   <  2:1421334/150‑1 (MQ=255)
                  aTTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCa                                                                                                                                  >  2:119659/1‑150 (MQ=255)
                    tGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCAtt                                                                                                                                <  2:1133929/150‑1 (MQ=255)
                    tGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCAtt                                                                                                                                >  1:1700398/1‑150 (MQ=255)
                    tGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCAtt                                                                                                                                <  1:1429952/150‑1 (MQ=255)
                     gCTGGGTCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGTCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCAttt                                                                                                                               <  1:100010/150‑1 (MQ=255)
                     gCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCATATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCAttt                                                                                                                               >  1:144845/1‑150 (MQ=255)
                      cTGGGGCATTTCTTGACTACAATTAATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGCGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGAACTGATTACCAGAACGATCAGTCATTTa                                                                                                                              >  2:1686033/1‑150 (MQ=255)
                       tGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTaa                                                                                                                             <  1:776379/150‑1 (MQ=255)
                       tGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTAAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTaa                                                                                                                             <  2:1021361/150‑1 (MQ=255)
                        ggggCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAAc                                                                                                                            <  1:1710924/150‑1 (MQ=255)
                         gggCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACa                                                                                                                           >  1:535672/1‑150 (MQ=255)
                           gCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAAt                                                                                                                         <  2:408708/150‑1 (MQ=255)
                             aTTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCTCAAATTATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATgg                                                                                                                       <  1:1184156/150‑1 (MQ=255)
                             aTTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATgg                                                                                                                       <  2:1685412/150‑1 (MQ=255)
                              tttCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGc                                                                                                                      >  2:575236/1‑150 (MQ=255)
                              tttCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGc                                                                                                                      <  2:781232/150‑1 (MQ=255)
                              tttCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGc                                                                                                                      <  2:1594737/150‑1 (MQ=255)
                                   tGACTACAATATATATGCCTCTCATTATGCACCACATCAGGGCGAAAATTCGCAAAATATAAGGTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAg                                                                                                                 <  2:1101564/150‑1 (MQ=255)
                                       tACAATTTATATGCCTCTCATTATGCACCACATCATGGCGATAAATCGCAAAATATATGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAg                                                                                                             <  2:1062300/150‑1 (MQ=255)
                                        aCAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAgg                                                                                                            <  2:1409179/150‑1 (MQ=255)
                                         cAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGAGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGc                                                                                                           >  1:530038/1‑150 (MQ=255)
                                             ttATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGAcc                                                                                                       <  2:916548/150‑1 (MQ=255)
                                              tataTGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGGTAATCAGTACgcagg                                                                                                                                      <  1:78043/118‑5 (MQ=255)
                                              tataTGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGGTAATCAGTACgcagg                                                                                                                                      >  2:82514/1‑114 (MQ=255)
                                              tataTGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGGTAATCAGTACgcagg                                                                                                                                      >  2:78043/1‑114 (MQ=255)
                                              tataTGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGGTAATCAGTACgcagg                                                                                                                                      <  1:82514/118‑5 (MQ=255)
                                               atatGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCGGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCaa                                                                                                     <  2:1251714/150‑1 (MQ=255)
                                                tatGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCCACCAAc                                                                                                    >  1:1537965/1‑150 (MQ=255)
                                                 atGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAAcc                                                                                                   <  1:597247/150‑1 (MQ=255)
                                                  tGCCTCTCATTCTGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGTCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTCATTAACAGTAATATAAGTCATTTAACAATGGCAAAAGCCAGGCTCCGAACCt                                                                                                  >  1:273933/1‑150 (MQ=255)
                                                  tGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCt                                                                                                  >  1:1341003/1‑150 (MQ=255)
                                                    ccTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACATgg                                                                                                <  2:170259/150‑1 (MQ=255)
                                                    ccTCTCATTATGCACCACATCACGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAACCCAGGCGACCAACCTgg                                                                                                <  1:1671575/150‑1 (MQ=255)
                                                      tctcATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGAt                                                                                              <  2:583134/150‑1 (MQ=255)
                                                      tctcATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTAGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATATAACAATGGCAAAAGCCAGGCGACCAACCTGGAt                                                                                              >  1:356342/1‑150 (MQ=255)
                                                        tcATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTTGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGAttt                                                                                            >  2:842709/1‑150 (MQ=255)
                                                          aTTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTc                                                                                          >  1:1277136/1‑150 (MQ=255)
                                                            tATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCcg                                                                                        <  1:326045/150‑1 (MQ=255)
                                                                caccacATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGta                                                                                    >  2:747618/1‑150 (MQ=255)
                                                                 accacATCAGGGCGATAATTCGCAAAAAATAAGTTCCTATGGGCAGGCTGGGGTTGATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCAATTAACAATGGCAAAAGCCAGGCGACAAACCTGGATATTCCGCGtat                                                                                   >  1:650311/1‑150 (MQ=255)
                                                                  ccacATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtatt                                                                                  >  1:1298413/1‑150 (MQ=255)
                                                                     caTCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtatttat                                                                               >  1:1215041/1‑150 (MQ=255)
                                                                      aTCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtatttatt                                                                              >  1:1224689/1‑150 (MQ=255)
                                                                             cGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCg                                                                       >  1:1263022/1‑150 (MQ=255)
                                                                               aTAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGcc                                                                     <  1:1027012/150‑1 (MQ=255)
                                                                               aTAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGcc                                                                     <  1:1336081/150‑1 (MQ=255)
                                                                                tAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGccc                                                                    <  1:1050438/150‑1 (MQ=255)
                                                                                      gCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATccc                                                              >  2:1364180/1‑150 (MQ=255)
                                                                                      gCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATccc                                                              >  2:1612849/1‑150 (MQ=255)
                                                                                       cAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGAGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCa                                                             >  2:1646695/1‑150 (MQ=255)
                                                                                       cAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTAGGGGCCTCGCGCCTGCGTACTGATTACCACTACGAACAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCa                                                             >  1:1663217/1‑150 (MQ=255)
                                                                                          aaTATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCa                                                          <  2:528438/150‑1 (MQ=255)
                                                                                          aaTATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCa                                                          >  1:1530876/1‑150 (MQ=255)
                                                                                               aaGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaatgaa                                                     <  1:695876/150‑1 (MQ=255)
                                                                                                    ccTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCaa                                                >  2:937485/1‑150 (MQ=255)
                                                                                                       aTGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTGCGCCCAATCCCAGCAATGAATGCAAaac                                             <  1:240648/150‑1 (MQ=255)
                                                                                                          ggCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaactaa                                          >  1:676306/1‑150 (MQ=255)
                                                                                                            cAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGGCCAGCCTGGATTTTCCGCGTATTTATTTGGTTCGCCCAATCCCAGCATTGAATGCAAgactaact                                        >  2:421077/1‑150 (MQ=255)
                                                                                                            cAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaactatct                                        >  1:180633/1‑150 (MQ=255)
                                                                                                              ggCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACtat                                      >  1:978971/1‑150 (MQ=255)
                                                                                                              ggCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACtat                                      >  2:1719966/1‑150 (MQ=255)
                                                                                                              ggCTGGGGTTAATCTTGGGGCCTGGAGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCTACCAACCTTGATTTTACGCGTATTTATTTGTTTAGCACAATCCCAGCAATGAATGCAAAAATAACtat                                      >  2:1713459/1‑150 (MQ=255)
                                                                                                               gCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACtata                                     >  1:1381258/1‑150 (MQ=255)
                                                                                                               gCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACtata                                     <  1:460716/150‑1 (MQ=255)
                                                                                                                  ggggTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGt                                  <  1:1364180/150‑1 (MQ=255)
                                                                                                                   gggTTAATCTTTGGGCCTGGCGGCTGCGTACTGATTACCAGTACGATCAGTCATTTTACAATGGCAAAATCCAGGAGACCAAGCTGGATTTTCCGCGTATTTATTTGTTTCTCCCAATCCCTGCTATGTATGCAAAACTAACTATAGGTc                                 <  1:1089317/150‑1 (MQ=255)
                                                                                                                   gggTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTc                                 >  2:1124422/1‑150 (MQ=255)
                                                                                                                   gggTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTc                                 >  2:212209/1‑150 (MQ=255)
                                                                                                                   gggTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTc                                 <  2:704917/150‑1 (MQ=255)
                                                                                                                   gggTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTc                                 <  2:144845/150‑1 (MQ=255)
                                                                                                                    ggTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCa                                >  2:1712373/1‑150 (MQ=255)
                                                                                                                           cTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGAt                         >  1:1465363/1‑150 (MQ=255)
                                                                                                                               gggCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAATCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTg                     <  2:530038/150‑1 (MQ=255)
                                                                                                                                    tGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATcc                >  1:1033291/1‑150 (MQ=255)
                                                                                                                                      gcgcCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCctc              <  2:1046033/150‑1 (MQ=255)
                                                                                                                                       cgcCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGACACTGAATCctct             >  2:188369/1‑150 (MQ=255)
                                                                                                                                       cgcCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCACAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCctct             >  1:1519179/1‑150 (MQ=255)
                                                                                                                                         tgatcGTACTGGTAATCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTAt           >  1:1316305/5‑150 (MQ=255)
                                                                                                                                          gatcGTACTGGTAATCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaacta                                           <  1:1331064/114‑1 (MQ=255)
                                                                                                                                          gatcGTACTGGTAATCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaacta                                           <  1:849175/114‑1 (MQ=255)
                                                                                                                                          gatcGTACTGGTAATCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaacta                                           >  2:1331064/4‑117 (MQ=255)
                                                                                                                                          gatcGTACTGGTAATCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTAtt          >  2:69253/4‑150 (MQ=255)
                                                                                                                                            gCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGAATTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTAtttt        <  1:325526/150‑1 (MQ=255)
                                                                                                                                             cGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTc       <  2:533864/150‑1 (MQ=255)
                                                                                                                                              gTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCg      <  2:1298413/150‑1 (MQ=255)
                                                                                                                                               tACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGa     >  1:1116821/1‑150 (MQ=255)
                                                                                                                                               tACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGa     >  1:981762/1‑150 (MQ=255)
                                                                                                                                                aCTGATTACCAGTACGATCTGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGAc    >  2:967062/1‑150 (MQ=255)
                                                                                                                                                  tGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGActc  >  1:805721/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTC  >  NC_000913/3187841‑3188136

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: