Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,318,470 G→A 8.5% A440V (GCA→GTA) ‡ trpC ← indole‑3‑glycerolphosphate synthetase and N‑(5‑phosphoribosyl)anthranilate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,318,4700GA 8.5% 258.1 / 11.2 82A440V (GCA→GTA) ‡trpCindole‑3‑glycerolphosphate synthetase and N‑(5‑phosphoribosyl)anthranilate isomerase
Reads supporting (aligned to +/- strand):  ref base G (38/37);  new base A (4/3);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

AAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGT  >  NC_000913/1318322‑1318617
                                                                                                                                                    |                                                                                                                                                   
aaGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTg                                                                                                                                                     >  2:1585451/1‑149 (MQ=255)
    ttcttcCAGCTGGCGCAGAGCAGGCATCAGGATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTc                                                                                                                                                 >  1:801488/1‑149 (MQ=255)
     tcttcCAGCTGGCGCAGAGCAGGCATCAGGATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCt                                                                                                                                                <  1:1585451/149‑1 (MQ=255)
      cttcCAGCTGGCGCAGAGCAGGCATCAGGATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCtt                                                                                                                                               >  2:1857239/1‑149 (MQ=255)
        tcCAGCTGGCGCAGAGCAGGCATCAGGATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGCACGTCtttt                                                                                                                                             <  3:1866664/149‑2 (MQ=255)
         ccAGCTGGCGCAGAGCAGGCATCAGGATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGa                                                                                                                                            >  2:538212/1‑149 (MQ=255)
                      agCAGGCATCAGGATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGa                                                                                                                               <  1:1156843/149‑1 (MQ=255)
                       gCAGGCATCAGGATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGAc                                                                                                                              >  1:685853/1‑149 (MQ=255)
                           gCATCAGGATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTa                                                                                                                          >  4:993739/1‑149 (MQ=255)
                           gCATCAGGATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTa                                                                                                                          <  1:1857239/149‑1 (MQ=255)
                             aTCAGGATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTAca                                                                                                                        <  2:785129/149‑1 (MQ=255)
                             aTCAGGATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTAca                                                                                                                        >  1:785129/1‑149 (MQ=255)
                              tCAGGATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGc                                                                                                                                                    <  2:539339/120‑1 (MQ=255)
                              tCAGGATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGc                                                                                                                                                    >  1:539339/1‑120 (MQ=255)
                                  gATTTGTGGCACGTACATGCCGCCAAACTCATCAAAATAGGGGTTAAGTAATGTTATCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGa                                                                                                                   <  1:538212/149‑1 (MQ=255)
                                  gATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTAca                                                                                                                        <  1:708200/144‑1 (MQ=255)
                                  gATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTAca                                                                                                                        >  2:708200/1‑144 (MQ=255)
                                        tgGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGActct                                                                                                                           <  1:486097/135‑1 (MQ=255)
                                        tgGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGActct                                                                                                                           >  2:486097/1‑135 (MQ=255)
                                        tgGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTaa                                                                                                              >  1:1115614/1‑148 (MQ=255)
                                        tgGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTaa                                                                                                              <  2:1115614/148‑1 (MQ=255)
                                         ggCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTaaaa                                                                                                            <  2:801488/149‑1 (MQ=255)
                                            aCGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGa                                                                                                                               <  4:1691538/127‑1 (MQ=255)
                                            aCGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGa                                                                                                                               >  3:1691538/1‑127 (MQ=255)
                                            aCGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGActc                                                                                                                            >  1:423499/1‑130 (MQ=255)
                                            aCGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGActc                                                                                                                            <  2:423499/130‑1 (MQ=255)
                                             cGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCaa                                                                                                        <  2:138631/149‑1 (MQ=255)
                                                    gccgccAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGgc                                                                                                 <  2:935060/149‑1 (MQ=255)
                                                          aaaCTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGaa                                                                                                                  >  2:1142661/1‑126 (MQ=255)
                                                          aaaCTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGaa                                                                                                                  <  1:1142661/126‑1 (MQ=255)
                                                               cACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGttt                                                                                      <  2:685853/149‑1 (MQ=255)
                                                                aCCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTtg                                                                                     >  4:1825113/1‑149 (MQ=255)
                                                                   aaaaTAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGc                                                                                  <  4:189783/149‑1 (MQ=255)
                                                                    aaaTAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGCACGTCttttg                                                                                                                                            >  2:1699149/1‑88 (MQ=255)
                                                                    aaaTAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGCACGTCttttg                                                                                                                                            <  1:1699149/90‑3 (MQ=255)
                                                                       tAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGAt                                                                                                                                           >  1:1048825/1‑88 (MQ=255)
                                                                       tAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGAt                                                                                                                                           <  2:1048825/88‑1 (MQ=255)
                                                                            gTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGtt                                                                                                                                   >  1:1750585/1‑91 (MQ=255)
                                                                            gTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGtt                                                                                                                                   <  2:1750585/91‑1 (MQ=255)
                                                                             ttAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACg                                                                        <  3:1825113/149‑1 (MQ=255)
                                                                                     tgttgtCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGaa                                                                                                                  <  1:900508/99‑1 (MQ=255)
                                                                                     tgttgtCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGaa                                                                                                                  >  2:900508/1‑99 (MQ=255)
                                                                                        tgtCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTgc                                                             >  2:128889/1‑149 (MQ=255)
                                                                                               gTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAg                                                      >  3:177077/1‑149 (MQ=255)
                                                                                                   ctttccttAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGcccc                                                  >  4:429343/1‑149 (MQ=255)
                                                                                                        cttAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGTCGCAGCCGGtt                                                                                       >  4:1162979/1‑107 (MQ=255)
                                                                                                        cttAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGtt                                                                                       <  3:1162979/107‑1 (MQ=255)
                                                                                                        cttAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGc                                                                               <  3:939580/115‑1 (MQ=255)
                                                                                                        cttAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGc                                                                               >  4:939580/1‑115 (MQ=255)
                                                                                                        cttAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGcc                                             >  4:896298/1‑149 (MQ=255)
                                                                                                         ttAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCca                                            >  1:507780/1‑149 (MQ=255)
                                                                                                             tatGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCcagcag                                        >  1:794477/1‑149 (MQ=255)
                                                                                                                 cgcgcAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACg                                    >  2:877791/1‑149 (MQ=255)
                                                                                                                  gcgcAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGt                                   >  3:959862/1‑149 (MQ=255)
                                                                                                                     cAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGc                                <  3:896298/149‑1 (MQ=255)
                                                                                                                      aGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGcc                               <  2:794477/149‑1 (MQ=255)
                                                                                                                          tCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGc                           <  1:128889/149‑1 (MQ=255)
                                                                                                                            tGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACg                                                                        >  4:488979/1‑102 (MQ=255)
                                                                                                                            tGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACg                                                                        <  3:488979/102‑1 (MQ=255)
                                                                                                                            tGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGCGa                         >  2:253580/1‑149 (MQ=255)
                                                                                                                             ggAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGCGAt                        >  2:983224/1‑149 (MQ=255)
                                                                                                                              gAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGCGAtt                       >  2:1044071/1‑149 (MQ=255)
                                                                                                                              gAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTATGCCGCTTCCACGCAGTTATCt                                                               >  1:65210/1‑109 (MQ=255)
                                                                                                                              gAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTATGCCGCTTCCACGCAGTTATCt                                                               <  2:65210/109‑1 (MQ=255)
                                                                                                                                aaaCCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCAc                                                                         <  3:338362/97‑1 (MQ=255)
                                                                                                                                aaaCCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCAc                                                                         >  4:338362/1‑97 (MQ=255)
                                                                                                                                aaaCCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCAc                                                                         >  2:1555006/1‑97 (MQ=255)
                                                                                                                                aaaCCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCAc                                                                         <  1:1555006/97‑1 (MQ=255)
                                                                                                                                    cGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGCGATTGAcctt                 <  3:993739/149‑3 (MQ=255)
                                                                                                                                    cGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGCGATTGACCAt                 >  1:1724843/1‑149 (MQ=255)
                                                                                                                                      aGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGCGATTGACCAttt               <  2:507780/149‑1 (MQ=255)
                                                                                                                                        gCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGCGATTGACCATTTaa             <  4:959862/149‑1 (MQ=255)
                                                                                                                                         atcaaaGACGCGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAAt            >  1:988360/4‑149 (MQ=255)
                                                                                                                                         atcaaaGACGCACGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCca                                            >  4:1375048/4‑117 (MQ=255)
                                                                                                                                         atcaaaGACGCACGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCca                                            <  3:1375048/114‑1 (MQ=255)
                                                                                                                                         atcaaaGACGCACGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAAt            >  3:824239/4‑149 (MQ=255)
                                                                                                                                           aaaaGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGc                                                                                  >  4:1460114/1‑77 (MQ=255)
                                                                                                                                           aaaaGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGc                                                                                  <  3:1460114/77‑1 (MQ=255)
                                                                                                                                            aaaGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGt         <  4:177077/149‑1 (MQ=255)
                                                                                                                                            aaaGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGt         <  1:983224/149‑1 (MQ=255)
                                                                                                                                            aaaGACGCACGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGt         >  4:920956/1‑149 (MQ=255)
                                                                                                                                               gACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGAc      <  3:429343/149‑1 (MQ=255)
                                                                                                                                                   tGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGt  >  4:1599766/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AAGCTTCTTCCAGCTGGCGCAGAGCAGGCATCAGGATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGT  >  NC_000913/1318322‑1318617

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