Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,150,057 T→C 5.6% S319G (AGC→GGC)  gapA ← glyceraldehyde‑3‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,150,0570TC 5.6% 207.5 / 2.2 71S319G (AGC→GGC) gapAglyceraldehyde‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (31/36);  new base C (2/2);  total (33/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 = 8.98e-01

CGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAGCGGGTTGTGGTTGAAATCGCAGGAAACCAATG  >  NC_002947/1149909‑1150176
                                                                                                                                                    |                                                                                                                       
cgATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCt                                                                                                                         <  4:10905/149‑1 (MQ=255)
   tCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGaa                                                                                                                      <  7:294137/149‑1 (MQ=255)
     ccTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAAc                                                                                                                     <  5:125496/148‑1 (MQ=255)
      cTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAAccc                                                                                                                   >  3:231627/1‑149 (MQ=255)
       tGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAAcccc                                                                                                                  >  2:253695/1‑149 (MQ=255)
           ggAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACt                                                                                                              <  1:215608/149‑1 (MQ=255)
               cGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGtt                                                                                                          <  1:253695/149‑1 (MQ=255)
               cGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGtt                                                                                                          >  6:177757/1‑149 (MQ=255)
                gggCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTg                                                                                                         >  5:196939/1‑149 (MQ=255)
                gggCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGGCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGCACCCCCACTCGTTg                                                                                                         >  5:192003/1‑149 (MQ=255)
                  gCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGtt                                                                                                          <  6:331531/146‑1 (MQ=255)
                  gCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGtt                                                                                                          >  5:331531/1‑146 (MQ=255)
                    ttGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCa                                                                                                      >  1:150400/1‑148 (MQ=255)
                    aagCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCTGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCa                                                                                                      >  6:171750/3‑148 (MQ=255)
                         cgcgAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAg                                                                                                >  5:953/1‑149 (MQ=255)
                           cgAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGc                                                                                              <  5:296519/149‑1 (MQ=255)
                              acacGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCa                                                                                            <  6:152701/148‑1 (MQ=255)
                              acacGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCACACCAGGCCAg                                                                                           >  7:201871/1‑149 (MQ=255)
                                  gggCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCCGAACCCCCACTCGTTGTCATACCAGGCCAGCAcc                                                                                       >  1:300796/1‑149 (MQ=255)
                                  gggCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCCGAACCCCCACTCGTTGTCATACCAGGCCAGCAcc                                                                                       >  1:295938/1‑149 (MQ=255)
                                     cGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAg                                                                                           >  4:97739/1‑142 (MQ=255)
                                     cGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAg                                                                                           <  3:97739/142‑1 (MQ=255)
                                      gAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTGGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCGACTCGTTGTCATACCAGGCCAGCACCTTGa                                                                                   <  2:181313/149‑1 (MQ=255)
                                      gAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGa                                                                                   <  2:150400/149‑1 (MQ=255)
                                                ccATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGccc                                                                         <  4:231627/149‑1 (MQ=255)
                                                ccATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGccc                                                                         >  1:163421/1‑149 (MQ=255)
                                                  aTCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGt                                                                       >  5:8370/1‑149 (MQ=255)
                                                     caccacGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTgg                                                                    >  6:64074/1‑149 (MQ=255)
                                                       ccacGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTg                                                                     <  2:175115/146‑1 (MQ=255)
                                                       ccacGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTg                                                                     >  1:175115/1‑146 (MQ=255)
                                                          tggcacgggCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCgg                                                               <  5:347803/142‑1 (MQ=255)
                                                             tgtTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTa                                                            >  1:243927/1‑149 (MQ=255)
                                                               ttGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTg                                                                                    >  4:253267/1‑123 (MQ=255)
                                                               ttGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTg                                                                                    <  3:253267/123‑1 (MQ=255)
                                                               ttGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATc                                                                              >  8:199312/1‑129 (MQ=255)
                                                               ttGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATc                                                                              <  7:199312/129‑1 (MQ=255)
                                                                cgtgttCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATgg                                                         <  2:160553/145‑1 (MQ=255)
                                                                  gCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGtt                                                       <  3:93705/149‑1 (MQ=255)
                                                                    ttCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGtt                                                       <  1:333446/147‑1 (MQ=255)
                                                                    ttCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCCTTGGCCCGGGTATGGtt                                                       >  2:333446/1‑147 (MQ=255)
                                                                           gCAAGGCCGCTCCCACACGGGTCGCGTCAGNCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCg                                                >  1:33088/1‑147 (MQ=255)
                                                                             aaGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGTAGCGCCAGGCAGTTGTCCAACATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGa                                                                                   <  4:10932/110‑1 (MQ=255)
                                                                                  cGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATcgac                                       >  1:279821/1‑149 (MQ=255)
                                                                                  cGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATcgac                                       <  6:953/149‑1 (MQ=255)
                                                                                  cGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCCGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATcgac                                       <  2:300796/149‑1 (MQ=255)
                                                                                  cGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCCGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATcgac                                       <  2:295938/149‑1 (MQ=255)
                                                                                    cTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATcgacga                                     >  2:134413/1‑149 (MQ=255)
                                                                                       ccACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAg                                  >  4:189898/1‑149 (MQ=255)
                                                                                        cacaCGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAGc                                 >  6:7869/1‑149 (MQ=255)
                                                                                         acacGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAGCg                                >  3:276359/1‑149 (MQ=255)
                                                                                         acacGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAGCg                                <  5:177757/149‑1 (MQ=255)
                                                                                          cacGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAGCgg                               <  6:8370/149‑1 (MQ=255)
                                                                                          cacGGGTCGCGGCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAGCgg                               <  6:192003/149‑1 (MQ=255)
                                                                                          cacGGGTCGCGGCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAGCgg                               <  6:196939/149‑1 (MQ=255)
                                                                                                tCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAGCGGGTtgtg                         <  2:279821/149‑1 (MQ=255)
                                                                                                 cgcgTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAGCGGGTTGTgg                        <  3:189898/149‑1 (MQ=255)
                                                                                                  gcgTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAGCGGGTTGTgg                        <  1:134413/148‑1 (MQ=255)
                                                                                                    gTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCACACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAGCGGGTTGTGGTTg                     <  8:201871/149‑1 (MQ=255)
                                                                                                     tCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAGCGGGTTGTGGTTGa                    <  5:64074/149‑1 (MQ=255)
                                                                                                       agccagccCGGTGCAGCGCCTGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCATTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAGCGGGTTGTGGTTGaa                   <  5:171750/148‑1 (MQ=255)
                                                                                                        gccagccCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCg                                                <  4:305275/118‑1 (MQ=255)
                                                                                                        gccagccCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCg                                                >  3:305275/1‑118 (MQ=255)
                                                                                                         ccagccCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATcgac                                       >  3:247308/1‑126 (MQ=255)
                                                                                                         ccagccCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATcgac                                       <  4:247308/126‑1 (MQ=255)
                                                                                                          cagccCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACa                                   <  7:75327/129‑1 (MQ=255)
                                                                                                          cagccCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACa                                   >  8:75327/1‑129 (MQ=255)
                                                                                                            gccCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAGCGGGTTGTGGTTGAAATCGCa             <  2:33088/149‑1 (MQ=255)
                                                                                                             cccGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAGCGGGTTGTGGTTGAAATCGCAg            >  5:343890/1‑149 (MQ=255)
                                                                                                                 gTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAGCGGGTTGTGGTTGAAATCGCAGGaaa        <  4:276359/149‑1 (MQ=255)
                                                                                                                    cAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAGCGGGTTGTgg                        <  4:159701/130‑1 (MQ=255)
                                                                                                                    cAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAGCGGGTTGTgg                        >  3:159701/1‑130 (MQ=255)
                                                                                                                       cgcCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAGCGGGTTGTGGTTCAAATCGCAGGAAACCAATg  >  8:240376/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                       
CGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCATACCAGGCCAGCACCTTGAGCATCCGCCCGTTGGCCCGGGTATGGTTGGCATCGAAGATCGACGACAGCGGGTTGTGGTTGAAATCGCAGGAAACCAATG  >  NC_002947/1149909‑1150176

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