Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,178,848 A→T 5.6% K86* (AAA→TAA) ‡ nagK → N‑acetyl‑D‑glucosamine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,178,8480AT 5.6% 253.3 / 5.1 71K86* (AAA→TAA) ‡nagKN‑acetyl‑D‑glucosamine kinase
Reads supporting (aligned to +/- strand):  ref base A (34/33);  new base T (2/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

ACGCATTTTTAGATGCAGTGTGTGAGCTGGTAGCTGAAGCTGACCAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGGATGATGAATTTACTCAATATCCACTGGTGATGGGGTTGATTCTCGGCACCG  >  NC_000913/1178702‑1178989
                                                                                                                                                  |                                                                                                                                             
aCGCATTTTTAGATGCAGTGTGTGAGCTGGTAGCTGAAGCTGACCAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTaaa                                                                                                                                             <  1:314157/149‑1 (MQ=255)
  gCATTTTTAGATGCAGTGTGTGAGCTGGTAGCTGAAGCTGACCAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAAcc                                                                                                                                           >  2:141125/1‑149 (MQ=255)
  gCATTTTTAGATGCAGTGTGTGAGCTGGTAGCTGAAGCTGACCAACGTTTTGGCTGTAAAGGTACTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAAcc                                                                                                                                           >  3:356141/1‑149 (MQ=255)
    aTTTTTAGATGCAGTGTGTGAGTTGGTAGCTGAAGCTGACCAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGc                                                                                                                                         <  3:504094/149‑1 (MQ=255)
     tttttAGATGCAGTGTGTGAGCTGGTAGCTGAAGCTGACCAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCt                                                                                                                                        >  3:388126/1‑149 (MQ=255)
       tttAGATGCAGTGTGTGAGCTGGTAGCTGAAGCTGACCAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGCTTTACCGCTgg                                                                                                                                      <  3:536484/149‑2 (MQ=255)
        ttAGATGCAGTGTGTGAGCTGGTAGCTGAAGCTGACCAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCg                                                                                                                                     >  2:559295/1‑149 (MQ=255)
               cAGTGTGTGAGCTGGTAGCTGAAGCTGACCAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGAc                                                                                                                              >  4:372448/1‑149 (MQ=255)
                     gtgAGCTGGTAGCTGAAGCTGACCAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAgc                                                                                                                        <  4:314019/149‑1 (MQ=255)
                       gAGCTGGTAGCTGAAGCTGACCAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAgcgc                                                                                                                      >  1:329109/1‑149 (MQ=255)
                       gAGCTGGTAGCTGAAGCTGACCAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAgcgc                                                                                                                      >  4:235251/1‑149 (MQ=255)
                        aGCTGGTAGCTGAAGCTGACCAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCa                                                                                                                     >  4:304240/1‑149 (MQ=255)
                           tGGTAGCTGAAGCTGACCAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGt                                                                                                                  >  4:502656/1‑149 (MQ=255)
                                  tgaagctgaCCAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATc                                                                                                           >  2:434132/1‑149 (MQ=255)
                                  tgaagctgaCCAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATc                                                                                                           <  4:599814/149‑1 (MQ=255)
                                       ctgaCCAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGAt                                                                                                      <  2:188202/149‑1 (MQ=255)
                                         gaCCAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGt                                                                                                    <  1:341445/149‑1 (MQ=255)
                                            cAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTTTACCGCTgg                                                                                                                                      >  2:433260/1‑111 (MQ=255)
                                            cAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTTTACCGCTgg                                                                                                                                      <  1:433260/112‑2 (MQ=255)
                                            cAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGc                                                                                                                                         >  1:105734/1‑109 (MQ=255)
                                            cAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGc                                                                                                                                         <  2:105734/109‑1 (MQ=255)
                                            cAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACg                                                                                                 >  1:523687/1‑149 (MQ=255)
                                                 ttttGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTg                                                                                            >  4:79107/1‑149 (MQ=255)
                                                  tttGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGa                                                                                           <  2:232254/149‑1 (MQ=255)
                                                     ggCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATaa                                                                                        <  2:130882/149‑1 (MQ=255)
                                                      gCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGCTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAAc                                                                                       <  1:447695/149‑1 (MQ=255)
                                                       cTGTAAAGGTTCTGTCGACATCGGTATTCCGGGAATGCCGGAAACAGAACATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGTGCGCACGTCTTGATCGCGATGTACGCCTTGATAACg                                                                                      <  4:523521/149‑1 (MQ=255)
                                                        tGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGa                                                                                     <  2:236336/149‑1 (MQ=255)
                                                          tAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGa                                                                                     <  1:212223/147‑1 (MQ=255)
                                                          tAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGa                                                                                     >  2:212223/1‑147 (MQ=255)
                                                              ggTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGTCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCaa                                                                               >  1:95799/1‑149 (MQ=255)
                                                                    gTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCTAATGTCCCGGCTGCCAGCGGAAAACGGATGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGttt                                                                         >  2:86667/1‑149 (MQ=255)
                                                                    gTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGa                                                                                                             <  2:73558/113‑1 (MQ=255)
                                                                    gTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGa                                                                                                             >  1:73558/1‑113 (MQ=255)
                                                                    gTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGttt                                                                         >  4:557001/1‑149 (MQ=255)
                                                                     tCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGtttt                                                                        <  1:122705/149‑1 (MQ=255)
                                                                     tCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGtttt                                                                        <  1:396223/149‑1 (MQ=255)
                                                                     tCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGtttt                                                                        >  1:63457/1‑149 (MQ=255)
                                                                        gCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGcc                                                                     >  4:341494/1‑149 (MQ=255)
                                                                        gCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGcc                                                                     >  2:234518/1‑149 (MQ=255)
                                                                        gCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGcc                                                                     <  3:372448/149‑1 (MQ=255)
                                                                        gCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGcc                                                                     >  2:134219/1‑149 (MQ=255)
                                                                         cATCGGTATTCCGGGAATGCCGGAAACAGATGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGccc                                                                    <  1:404466/149‑1 (MQ=255)
                                                                            cGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCttt                                                                 <  4:388126/149‑1 (MQ=255)
                                                                            cGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCttt                                                                 >  3:241309/1‑149 (MQ=255)
                                                                               tATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAg                                                              <  4:241309/149‑1 (MQ=255)
                                                                                      ggAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAGCGATGCCAACTGTTTTGCCCTGTCAGAAGCCTg                                                       <  2:523687/149‑1 (MQ=255)
                                                                                        aaTGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTggg                                                     <  3:235251/149‑1 (MQ=255)
                                                                                            ccGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGgatga                                                 <  1:141125/149‑1 (MQ=255)
                                                                                             cGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGgatgat                                                >  4:283792/1‑149 (MQ=255)
                                                                                             cGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGgatgat                                                >  1:374216/1‑149 (MQ=255)
                                                                                               gAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGGAtgatga                                              >  1:366051/1‑149 (MQ=255)
                                                                                                 aaCAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGGATGATGAAt                                            <  2:374216/149‑1 (MQ=255)
                                                                                                     gaagaTGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGtttt                                                                        >  3:172911/1‑117 (MQ=255)
                                                                                                     gaagaTGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGtttt                                                                        <  4:172911/117‑1 (MQ=255)
                                                                                                     gaagaTGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGGATGATGAATTTAc                                        >  1:17536/1‑149 (MQ=255)
                                                                                                        gaTGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGGATGATGAATTTACTCa                                     >  1:610851/1‑149 (MQ=255)
                                                                                                                  tGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGGATGATGAATTTACTCAATATCCACTg                           <  3:304240/149‑1 (MQ=255)
                                                                                                                        ccgccAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGGATGATGAATTTACTCAATATCCACTGGTGATg                     >  2:577232/1‑149 (MQ=255)
                                                                                                                          gccAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGGATGATGAATTTACTCAATATCCACTGGTGATggg                   <  1:234518/149‑1 (MQ=255)
                                                                                                                              aTGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGGATGATGAATTTACTCAATATCCACTGGTGATGGGGTTg               <  1:57456/149‑1 (MQ=255)
                                                                                                                                   cccGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGGATGATGAATTTACTCAATATCCACTGGTGATGGGGTTGATtct          <  1:434132/149‑1 (MQ=255)
                                                                                                                                    ccGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAg                                                           <  1:356872/99‑1 (MQ=255)
                                                                                                                                    ccGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAg                                                           >  2:356872/1‑99 (MQ=255)
                                                                                                                                    ccGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGGATGATGAATTTACTCAATATCCACTGGTGATGGGGTTGATtctc         <  2:366051/149‑1 (MQ=255)
                                                                                                                                     cGCTGCCAGCGGTAAACCGCTGCGTGTCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGGATGATGAATTTACTCAATATCCACTGGTGATGGGGTTGATTCTCg        <  2:95799/149‑1 (MQ=255)
                                                                                                                                          gcAGCGGTTTACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGGATGATGAATTTACTCAATATCCACTGGTGATGGGGTTGATTCTCGGCAcc   >  3:32983/2‑149 (MQ=255)
                                                                                                                                          ccAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGGATGATGAATTTACTCAATATCCACTGGTGATGGGGTTGATTCTCGGCAcc   >  2:340602/1‑149 (MQ=255)
                                                                                                                                          ccAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGGATGATGAATTTACTCAATATCCACTGGTGATGGGGTTGATTCTCGGCAcc   >  2:272664/1‑149 (MQ=255)
                                                                                                                                           cAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGGATGATGAATTTACTCAATATCCACTGGTGATGGGGTTGATTCTCGGCACCg  <  1:513475/149‑1 (MQ=255)
                                                                                                                                                gTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGGATGATGAATTTa                                         <  1:32874/105‑1 (MQ=255)
                                                                                                                                                gTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGGATGATGAATTTa                                         >  2:32874/1‑105 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
ACGCATTTTTAGATGCAGTGTGTGAGCTGGTAGCTGAAGCTGACCAACGTTTTGGCTGTAAAGGTTCTGTCGGCATCGGTATTCCGGGAATGCCGGAAACAGAAGATGGTACGCTGTATGCCGCCAATGTCCCCGCTGCCAGCGGTAAACCGCTGCGTGCCGACCTGAGCGCACGTCTTGATCGCGATGTACGCCTTGATAACGATGCCAACTGTTTTGCCCTTTCAGAAGCCTGGGATGATGAATTTACTCAATATCCACTGGTGATGGGGTTGATTCTCGGCACCG  >  NC_000913/1178702‑1178989

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