Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,985,112 C→G 11.9% intergenic (‑43/+27) araG ← / ← araF L‑arabinose ABC transporter ATPase/L‑arabinose ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,985,1120CG11.9% 182.9 / 16.7 59intergenic (‑43/+27)araG/araFL‑arabinose ABC transporter ATPase/L‑arabinose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (27/25);  new base G (5/2);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-01

CCTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTGTAAA  >  NC_000913/1984967‑1985260
                                                                                                                                                 |                                                                                                                                                     
ccTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATg                                                                                                                                                    >  1:340807/1‑149 (MQ=255)
     ccGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCaga‑g                                                                                                                                              <  1:217962/149‑1 (MQ=255)
       ggCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑ggg                                                                                                                                            <  2:103333/149‑1 (MQ=255)
       ggCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑ggg                                                                                                                                            <  2:419139/149‑1 (MQ=255)
          aTAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGaa                                                                                                                                         >  2:51693/1‑149 (MQ=255)
          aTAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGaa                                                                                                                                         <  4:487751/149‑1 (MQ=255)
                 tCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTcc                                                                                                                                  >  3:116312/1‑149 (MQ=255)
                    aaaCTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGc                                                                                                                               <  3:319840/149‑1 (MQ=255)
                          atatCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCaga‑g                                                                                                                                              >  4:352044/1‑128 (MQ=255)
                          atatCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCaga‑g                                                                                                                                              <  3:352044/128‑1 (MQ=255)
                                    gcgcCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCt                                                                                                               >  3:26637/1‑149 (MQ=255)
                                        cTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAAc                                                                                                           >  1:510441/1‑149 (MQ=255)
                                         ttAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAAcc                                                                                                          >  3:538812/1‑149 (MQ=255)
                                          tAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCt                                                                                                         >  1:92112/1‑149 (MQ=255)
                                                 gggAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTtct                                                                                                  >  2:439473/1‑149 (MQ=255)
                                                 gggAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTtct                                                                                                  >  1:28301/1‑149 (MQ=255)
                                                   gAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTcc                                                                                                >  1:263977/1‑149 (MQ=255)
                                                   gAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTcc                                                                                                >  3:373869/1‑149 (MQ=255)
                                                    aaaCGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCa                                                                                               <  4:321418/149‑1 (MQ=255)
                                                    aaaCGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCa                                                                                               <  2:154005/149‑1 (MQ=255)
                                                     aaCGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAg                                                                                              >  2:15281/1‑149 (MQ=255)
                                                                 aTGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTaaa                                                                                  >  1:45814/1‑149 (MQ=255)
                                                                           aTGAGAGATACGGGGTAGACTGTTGCATAATGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAcgc                                                                        >  2:486165/1‑149 (MQ=255)
                                                                            tGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAcgcg                                                                       >  3:427407/1‑149 (MQ=255)
                                                                             gagagaTACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGt                                                                      <  2:92112/149‑1 (MQ=255)
                                                                                     cGGGGTAGACTGTTGCATAATGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTGCTTACCGCCTAAACCTTTTTTCTCCAGTGCTTCTTTAAAGTTGTCACGCGTGATCAGTa                                                              <  1:486165/149‑1 (MQ=255)
                                                                                           aGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGt                                                        >  1:189588/1‑149 (MQ=255)
                                                                                            gACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTc                                                       <  4:538812/149‑1 (MQ=255)
                                                                                                    cATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTaaa                                                                                  >  2:575777/1‑114 (MQ=255)
                                                                                                    cATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTaaa                                                                                  <  1:575777/114‑1 (MQ=255)
                                                                                                    cATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACtt                                               <  3:294752/149‑1 (MQ=255)
                                                                                                    cATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACtt                                               <  1:520941/149‑1 (MQ=255)
                                                                                                      tAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCg                                             >  2:61489/1‑149 (MQ=255)
                                                                                                      tAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCg                                             >  4:525638/1‑149 (MQ=255)
                                                                                                       aaCGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCgg                                            <  1:15281/149‑1 (MQ=255)
                                                                                                           tctcCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTaaa                                        <  3:525638/149‑1 (MQ=255)
                                                                                                             tcCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAAtt                                      >  2:592525/1‑149 (MQ=255)
                                                                                                               cGTGAATCACTGGTCTTTCCCTCCCCCTCTGCATCATGCAGA‑GGGGAATTTTCCCGGCAAATTACTTACCTCCTAAAGCTTTTTTCTCCAGTTCTTCTTTAAAGTTCTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAAtttt                                    <  2:205907/149‑1 (MQ=255)
                                                                                                                gTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTg                                   <  4:422596/149‑1 (MQ=255)
                                                                                                                gTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTg                                   >  1:337154/1‑149 (MQ=255)
                                                                                                                gTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTg                                   <  4:116312/149‑1 (MQ=255)
                                                                                                                 tGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTgg                                  <  2:122227/149‑1 (MQ=255)
                                                                                                                    aTCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCgg                               >  3:371441/1‑149 (MQ=255)
                                                                                                                        cTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCa                           >  1:146417/1‑149 (MQ=255)
                                                                                                                             cGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTc                      <  2:223310/149‑1 (MQ=255)
                                                                                                                               ttCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCtt                    >  1:463237/1‑149 (MQ=255)
                                                                                                                                  acacCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTg                 >  1:455755/1‑149 (MQ=255)
                                                                                                                                    ttccccTCTGCATGATGCAGAGGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCtttt                  >  4:86739/3‑147 (MQ=255)
                                                                                                                                    ttccccTCTGCATGATGCAGAGGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCtttt                  <  3:86739/145‑1 (MQ=255)
                                                                                                                                    acCCCCTCTGCATCATGCAGC‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCt               >  4:366114/1‑149 (MQ=255)
                                                                                                                                    acCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGTCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCt               <  2:263977/149‑1 (MQ=255)
                                                                                                                                    acCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCt               <  4:373869/149‑1 (MQ=255)
                                                                                                                                      ccccTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAc             >  1:285346/1‑149 (MQ=255)
                                                                                                                                      ccccTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAc             <  2:288375/149‑1 (MQ=255)
                                                                                                                                       cccTCTGCATGATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAcc            >  3:161501/1‑149 (MQ=255)
                                                                                                                                       cccTCTGCATGATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAcc            >  1:260181/1‑149 (MQ=255)
                                                                                                                                        ccTCTGCATGATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAccc           >  2:33364/1‑149 (MQ=255)
                                                                                                                                            tGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGtt       <  1:61489/149‑1 (MQ=255)
                                                                                                                                            tGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGtt       <  2:146417/149‑1 (MQ=255)
                                                                                                                                              catcatGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTGt     <  1:51693/149‑1 (MQ=255)
                                                                                                                                               atcatGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTGTa    >  3:302563/1‑149 (MQ=255)
                                                                                                                                               atcatGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTGTa    >  2:25654/1‑149 (MQ=255)
                                                                                                                                                 catGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTGTaaa  <  4:26637/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CCTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTGTAAA  >  NC_000913/1984967‑1985260

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑