Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,612,716 C→T 5.1% G179D (GGC→GAC)  glsB ← glutaminase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,612,7160CT 5.1% 259.2 / 2.4 78G179D (GGC→GAC) glsBglutaminase 2
Reads supporting (aligned to +/- strand):  ref base C (38/36);  new base T (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01

TCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGTGGTGCGCTTA  >  NC_000913/1612573‑1612860
                                                                                                                                               |                                                                                                                                                
tCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCaaa                                                                                                                                              <  4:82176/148‑1 (MQ=255)
 caccacTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAAc                                                                                                                                             >  8:183460/1‑148 (MQ=255)
 caccacTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAAc                                                                                                                                             >  4:303209/1‑148 (MQ=255)
  accacTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACg                                                                                                                                            <  3:392750/148‑1 (MQ=255)
   ccacTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGa                                                                                                                                           >  1:312796/1‑148 (MQ=255)
    cacTGGTTCATCAATATGGATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGAc                                                                                                                                          <  2:234582/148‑1 (MQ=255)
     acTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACt                                                                                                                                         >  2:224863/1‑148 (MQ=255)
      cTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACtt                                                                                                                                        >  8:335304/1‑148 (MQ=255)
        ggTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTtca                                                                                                                                      >  4:121738/1‑148 (MQ=255)
               caATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGc                                                                                                                                 >  7:260801/1‑146 (MQ=255)
                 atatGAAAAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCACTACACAGCACATTTTAAGAGCGCAGTAATGAAAGAAGTATAAAAGAAAGGTTGTCACGTCATGATGGAAATTGCAAAACCACTTAATCAGCCAgt                                                                                                                             >  5:399889/1‑147 (MQ=255)
                  tatGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGc                                                                                                                            >  1:43728/1‑148 (MQ=255)
                          gCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCg                                                                                                                    <  1:224863/148‑1 (MQ=255)
                             ttCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCAt                                                                                                                 >  1:42281/1‑148 (MQ=255)
                                    gATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTcgcgc                                                                                                          >  4:140376/1‑148 (MQ=255)
                                     aTTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTcgcgc                                                                                                          <  2:42281/147‑1 (MQ=255)
                                      ttAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCgg                                                                                                        <  3:303209/148‑1 (MQ=255)
                                        aGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGa                                                                                                       <  2:312796/147‑1 (MQ=255)
                                          ccAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATg                                                                                                    >  5:151340/1‑148 (MQ=255)
                                            aGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGtt                                                                                                  >  4:391615/1‑148 (MQ=255)
                                             gAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGcc                                                                                                                                <  7:67468/117‑1 (MQ=255)
                                             gAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGcc                                                                                                                                >  8:67468/1‑117 (MQ=255)
                                               aaGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCa                                                                                                >  3:221546/1‑147 (MQ=255)
                                               aaGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCa                                                                                                >  8:188588/1‑147 (MQ=255)
                                                   cAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTc                                                                                           >  4:334984/1‑148 (MQ=255)
                                                      aCGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGa                                                                                        <  5:421506/148‑1 (MQ=255)
                                                             ggCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCt                                                                                 <  3:106607/148‑1 (MQ=255)
                                                              gCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTa                                                                                <  1:220077/148‑1 (MQ=255)
                                                               ccAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTac                                                                               <  1:327825/148‑1 (MQ=255)
                                                               ccAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTac                                                                               >  6:404809/1‑148 (MQ=255)
                                                                  gCTCTCCACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGGTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTaccac                                                                            <  3:301104/148‑1 (MQ=255)
                                                                   ctctACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTCTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTaccacc                                                                           >  7:258478/1‑148 (MQ=255)
                                                                        cacaGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTAt                                                                       <  4:221546/147‑1 (MQ=255)
                                                                          caGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCaa                                                                                                                                               >  5:292922/1‑73 (MQ=255)
                                                                          caGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCaa                                                                                                                                               <  6:292922/73‑1 (MQ=255)
                                                                          caGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCg                                                                     >  4:361745/1‑147 (MQ=255)
                                                                          caGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGt                                                                    <  5:325332/148‑1 (MQ=255)
                                                                              tCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAgg                                                                 >  1:435643/1‑147 (MQ=255)
                                                                                 ttttCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGaaa                                                              >  3:324072/1‑147 (MQ=255)
                                                                                 gtttCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTAGCGCGGAATGTTCAAATTCGGAAAGCGCTACCACCGTATCGTAGTaaa                                                              >  3:55577/2‑147 (MQ=255)
                                                                                     cAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATc                                                          <  4:324072/147‑1 (MQ=255)
                                                                                      agagCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGTCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAatat                                                           >  3:231382/1‑145 (MQ=255)
                                                                                      agagCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGTCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAatat                                                           >  7:276916/1‑145 (MQ=255)
                                                                                      agagCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAg                                                        <  2:13978/148‑1 (MQ=255)
                                                                                       gagCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGa                                                       <  7:335304/148‑1 (MQ=255)
                                                                                           gcAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTtca                                                                                                                                      <  2:316090/65‑1 (MQ=255)
                                                                                           gcAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTtca                                                                                                                                      >  1:316090/1‑65 (MQ=255)
                                                                                               tAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGTCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGc                                               <  4:231382/148‑1 (MQ=255)
                                                                                               tAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGTCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGc                                               <  8:276916/148‑1 (MQ=255)
                                                                                                aaTGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAAcgc                                                                                    >  2:444640/1‑110 (MQ=255)
                                                                                                aaTGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAAcgc                                                                                    <  1:444640/110‑1 (MQ=255)
                                                                                                aaTGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATc                                                                      >  7:440094/1‑124 (MQ=255)
                                                                                                aaTGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATc                                                                      <  8:440094/124‑1 (MQ=255)
                                                                                                aaTGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTa                                                                   <  4:308188/127‑1 (MQ=255)
                                                                                                aaTGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTa                                                                   >  3:308188/1‑127 (MQ=255)
                                                                                                 aTGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCtt                                             <  7:183460/148‑1 (MQ=255)
                                                                                                   gAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCtt                                             >  1:1179/1‑146 (MQ=255)
                                                                                                    aaaGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAg                                          >  2:389007/1‑148 (MQ=255)
                                                                                                     aaGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGc                                         >  7:255712/1‑148 (MQ=255)
                                                                                                       gtagtTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGc                                                                                                                   >  8:420181/1‑72 (MQ=255)
                                                                                                       gtagtTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGc                                                                                                                   <  7:420181/72‑1 (MQ=255)
                                                                                                       gtagtTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCcg                                       <  3:121738/148‑1 (MQ=255)
                                                                                                           ttttGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCa                                   <  5:404809/148‑1 (MQ=255)
                                                                                                            tttGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCAc                                  >  1:32491/1‑148 (MQ=255)
                                                                                                                  gAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGaa                                                               <  1:155299/113‑1 (MQ=255)
                                                                                                                  gAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGaa                                                               >  2:155299/1‑113 (MQ=255)
                                                                                                                   aaCGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTcc                           <  7:188588/148‑1 (MQ=255)
                                                                                                                   aaCGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTcc                           <  2:43728/148‑1 (MQ=255)
                                                                                                                    aCGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTc                            >  7:73835/1‑146 (MQ=255)
                                                                                                                      ggTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTCCAGc                        <  2:1179/148‑1 (MQ=255)
                                                                                                                           tCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTCCAGCATAc                    <  3:334984/147‑1 (MQ=255)
                                                                                                                                tCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCt                                                                                 >  4:91492/1‑81 (MQ=255)
                                                                                                                                tCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCt                                                                                 <  3:91492/81‑1 (MQ=255)
                                                                                                                                    gatgGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCg                                                                                          >  7:26292/1‑68 (MQ=255)
                                                                                                                                    gatgGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCg                                                                                          <  8:26292/68‑1 (MQ=255)
                                                                                                                                      tgGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAat                                                             >  7:8227/1‑95 (MQ=255)
                                                                                                                                      tgGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAat                                                             <  8:8227/95‑1 (MQ=255)
                                                                                                                                            ttGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGTGGTGCGCTTa  <  3:391615/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
TCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGAACGCGCTACCACCGTATCGTAGGAAATATCAGACACACCGCTTAAGCCGCGCACGACTTCCAGCATACGTTGCCGTGGTGCGCTTA  >  NC_000913/1612573‑1612860

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