Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,423,390 C→T 6.2% intergenic (+196/+33) yhdZ → / ← rrfF putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,3900CT 6.2% 998.2 / 35.3 306intergenic (+196/+33)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base C (139/148);  new base T (10/9);  total (149/157)
Fisher's exact test for biased strand distribution p-value = 8.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 318 total aligned reads displayed.

TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAG  >  NC_000913/3423241‑3423539
                                                                                                                                                     |                                                                                                                                                     
   cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAACTTAAAATACCCTCCCCGCGGGTTTTTACTAAACGGGAGAGTGTTAACCCACCAAAAACAGATAAAACAAAATTCCCac                                                                                                                                                    >  1:3146584/1‑149 (MQ=14)
           ttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGa                                                                                                                                            <  1:3778860/150‑1 (MQ=39)
             tttGTCTGGCTGTTTATGGGGGGCGTCCTGCCCGCCACCCTCCGGGCGGGTGCTTCGAAACGTTAAAATCCGCTCCCGGCGTATTTGTCCTACTCGGGAGTGTGTTCACCGACAAACAACATATAAAACAAAAGGCCCAGTCTTCCGACt                                                                                                                                          <  1:3300886/148‑1 (MQ=14)
               tGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGa                                                                                                                                        <  2:864179/150‑1 (MQ=39)
               tGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTAAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGa                                                                                                                                        <  1:353149/150‑1 (MQ=255)
                  cTGGCTGTTTCGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGTTTTGTCCTACTCGGGTGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGcc                                                                                                                                     <  1:2983444/150‑1 (MQ=11)
                    ggCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCTCAGTCGGAAGACTGGGCCtt                                                                                                                                   <  1:1672461/150‑1 (MQ=11)
                     gCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCttt                                                                                                                                  >  2:6079889/1‑150 (MQ=37)
                      cAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCTCAGTCGGAAGACTGGGCCtttt                                                                                                                                 <  2:3162384/150‑1 (MQ=11)
                      cAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt                                                                                                                                 <  2:302984/150‑1 (MQ=38)
                      cAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt                                                                                                                                 <  2:1723542/150‑1 (MQ=38)
                      cAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt                                                                                                                                 <  2:456174/150‑1 (MQ=38)
                       aGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTGTTg                                                                                                                                >  2:4026248/1‑150 (MQ=38)
                         tttATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtt                                                                                                                              <  1:4106109/150‑1 (MQ=39)
                          ttATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGTGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGttt                                                                                                                             <  1:2409298/150‑1 (MQ=255)
                           tATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTAAAATACGCGCCCGGCGGGTGTGTTCGTCTCGGGGGAGGGGTGACCGACAAACCACAGATAAAACAAAAGGCCCCGTCTTGCGAGGGGGCCGTGTGtttt                                                                                                                            >  1:4725504/1‑150 (MQ=2)
                             tGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAt                                                                                                                          <  1:6060948/150‑1 (MQ=39)
                              ggcgggcgCCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCGACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTAtc                                                                                                                         <  1:2554918/150‑2 (MQ=14)
                               gcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTAtct                                                                                                                        <  2:4066007/150‑3 (MQ=14)
                               gcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAttt                                                                                                                        <  1:3935109/150‑1 (MQ=39)
                                cgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTg                                                                                                                       <  1:3411698/150‑1 (MQ=39)
                                    cgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGt                                                                                                                   <  1:6216251/150‑1 (MQ=40)
                                       ccTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTGTATTTGCTGTCTg                                                                                                                >  2:1533982/1‑150 (MQ=255)
                                        cTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTGTTTTATTTGATGTCTGGc                                                                                                              <  2:3746643/150‑1 (MQ=255)
                                         tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc                                                                                                              >  2:514560/1‑150 (MQ=255)
                                         tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAGTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc                                                                                                              >  2:247844/1‑150 (MQ=255)
                                            ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGt                                                                                                           >  2:2233197/1‑150 (MQ=255)
                                                cACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTccc                                                                                                       <  1:1144944/150‑1 (MQ=255)
                                                cACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTccc                                                                                                       <  2:3137770/150‑1 (MQ=255)
                                                   ccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTAc                                                                                                    >  2:1949774/1‑150 (MQ=255)
                                                     tCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActc                                                                                                  <  2:2090670/150‑1 (MQ=255)
                                                     tCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActc                                                                                                  >  2:1387127/1‑150 (MQ=255)
                                                       cGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACtctc                                                                                                <  1:1582346/150‑1 (MQ=255)
                                                         gaacGTTGCGACGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACCACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCGTGTGTTTTAGTTGATGTCTGGCAGTTCCCTCCTCTCGc                                                                                              >  1:4596541/4‑150 (MQ=14)
                                                             gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAACCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATgg                                                                                          <  2:4789536/150‑1 (MQ=255)
                                                                gcgtCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGga                                                                                       <  1:4947925/147‑1 (MQ=255)
                                                                   tCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGAc                                                                                    <  1:3503855/150‑1 (MQ=255)
                                                                   tCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGAc                                                                                    >  2:3879758/1‑150 (MQ=255)
                                                                   tCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGAGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGAc                                                                                    <  1:6204264/150‑1 (MQ=255)
                                                                    cGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGAcc                                                                                   >  1:5525210/1‑150 (MQ=255)
                                                                     gCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTGTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGAccc                                                                                  >  1:5749139/1‑150 (MQ=255)
                                                                       aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCca                                                                                <  2:5355048/150‑1 (MQ=255)
                                                                       aaCGTTCAAATCCGCTCCCGGCGGAGGTGTCCTACTCGGGAGAGTGGTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCGTTTGTTTTATTTGATGTCTGGCTGTTCCCTACTCTCGCATGGGGAGACCCca                                                                                >  2:708650/1‑150 (MQ=14)
                                                                         cGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCcaca                                                                              <  2:6443755/150‑1 (MQ=255)
                                                                              aattactcTCCCGGAGTTTTTGTCCTTCTCGGGATAGTGTTCACCGACATCCAACTGATAAAAGAAAAGGCCCCGTCGTCGGACTTCGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTTcc                                                                         <  2:1180205/143‑1 (MQ=14)
                                                                               aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCa                                                                        <  2:5698527/150‑1 (MQ=255)
                                                                                    gCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGgc                                                                   <  2:5225859/150‑1 (MQ=255)
                                                                                      tCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGgcgc                                                                 >  1:740271/1‑150 (MQ=255)
                                                                                         cggcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAc                                                              <  2:5384997/150‑1 (MQ=255)
                                                                                           gcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACgg                                                            >  1:4204040/1‑150 (MQ=255)
                                                                                               aTTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGtt                                                        >  2:5729932/1‑150 (MQ=255)
                                                                                               aTTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGtt                                                        <  2:2948863/150‑1 (MQ=255)
                                                                                                  tGTCCTCCTCGGGAGAGTGTTCACCGACAACCAACAGATAAAACACAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCa                                                     <  2:3802057/150‑1 (MQ=255)
                                                                                                   gTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCAc                                                    >  1:4346757/1‑150 (MQ=255)
                                                                                                   gTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCAc                                                    >  1:629089/1‑150 (MQ=255)
                                                                                                       tACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGttt                                                       <  1:272432/143‑1 (MQ=255)
                                                                                                       tACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGttt                                                       >  2:272432/1‑143 (MQ=255)
                                                                                                       tACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                                <  2:3418855/150‑1 (MQ=255)
                                                                                                         cTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                                              <  1:5885181/150‑1 (MQ=255)
                                                                                                          tCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAg                                             <  1:1387127/150‑1 (MQ=255)
                                                                                                             ggAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTc                                          <  1:1494287/150‑1 (MQ=255)
                                                                                                               agagTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCgg                                        <  2:2822629/150‑1 (MQ=255)
                                                                                                                   tgtTCACCGACAAACACCAGCGACACCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCGACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCGACGGCGTTTCACTTCTGAGTTCGGCATg                                    <  1:4507327/150‑1 (MQ=14)
                                                                                                                    gtTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgg                                   >  1:6429404/1‑150 (MQ=39)
                                                                                                                     ttCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTGTCACTGCTGAGGTCGGCATggg                                  >  1:173466/1‑150 (MQ=255)
                                                                                                                        aCCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTc                               >  2:2746884/1‑150 (MQ=39)
                                                                                                                        aCCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTc                               <  2:5509309/150‑1 (MQ=39)
                                                                                                                        aCCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTc                               >  2:4095457/1‑150 (MQ=39)
                                                                                                                         ccGACAAACACCAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACGACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCa                              <  2:4218848/150‑1 (MQ=255)
                                                                                                                          cGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAg                             >  1:84429/1‑150 (MQ=39)
                                                                                                                          cGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAg                             <  2:3541160/150‑1 (MQ=39)
                                                                                                                            acaaacaaCAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt                           >  1:6677009/1‑150 (MQ=39)
                                                                                                                              aaacaaCAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTgg                         <  1:217781/150‑1 (MQ=39)
                                                                                                                              aaacaaCAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTgg                         <  2:3464961/150‑1 (MQ=39)
                                                                                                                              aaacaaCAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTgg                         <  2:4688925/150‑1 (MQ=39)
                                                                                                                                acaacaGATAAAACAACAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCGACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGa                       <  1:4787017/150‑1 (MQ=255)
                                                                                                                                acaacaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGa                       >  1:5974058/1‑150 (MQ=39)
                                                                                                                                 caacaGTTAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGac                      <  2:5645724/150‑1 (MQ=255)
                                                                                                                                 caacaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATggg                                                                                         >  1:3629164/1‑83 (MQ=255)
                                                                                                                                 caacaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGac                      <  1:2328949/150‑1 (MQ=39)
                                                                                                                                   acaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGacca                    >  2:5917789/1‑150 (MQ=39)
                                                                                                                                      gATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACcg                 >  2:4509167/1‑150 (MQ=39)
                                                                                                                                      aaTAAAACAAAAGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACcg                 >  1:103427/2‑150 (MQ=11)
                                                                                                                                       aTAAAACAAATGGCTCAGTCGGAAGTGTGGGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACcgc                >  2:365104/1‑150 (MQ=11)
                                                                                                                                       aTAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACcgc                <  1:2688789/150‑1 (MQ=39)
                                                                                                                                         aaaaCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCgcgc              <  2:5269562/150‑1 (MQ=39)
                                                                                                                                          aaaCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt                           >  2:798824/1‑136 (MQ=40)
                                                                                                                                           aaCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTa            >  1:4528321/1‑150 (MQ=39)
                                                                                                                                            aCAAAAGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAc           >  2:2916685/1‑150 (MQ=11)
                                                                                                                                            aCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAc           >  1:1203041/1‑150 (MQ=39)
                                                                                                                                             cAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACt          <  1:1964536/150‑1 (MQ=38)
                                                                                                                                              aaaaGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTAcc                                                                         <  1:5368105/86‑1 (MQ=14)
                                                                                                                                              aaaaGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCCCCGCGCTACtc         >  2:5044469/1‑149 (MQ=11)
                                                                                                                                              aaaaGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTg         >  2:4619464/1‑150 (MQ=38)
                                                                                                                                              aaaaGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTg         >  1:2276073/1‑150 (MQ=38)
                                                                                                                                                aaGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTgcc       <  2:6329493/150‑1 (MQ=37)
                                                                                                                                                   gCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTgccgcc    >  2:4950884/1‑150 (MQ=37)
                                                                                                                                                   gCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTgccgcc    <  2:5985593/150‑1 (MQ=37)
                                                                                                                                                   gCCCAGTCTTCCGACTGAGCCTTTTGTTTGATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTgccgcc    >  2:5873988/1‑150 (MQ=37)
                                                                                                                                                     ccAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAg  >  1:5497945/1‑150 (MQ=37)
                                                                                                                                                     |                                                                                                                                                     
TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAG  >  NC_000913/3423241‑3423539

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