Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,933,843 G→T 5.9% R254R (CGC→CGA)  edd ← 6‑phosphogluconate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,933,8430GT 5.9% 221.4 / 2.2 68R254R (CGC→CGA) edd6‑phosphogluconate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (40/24);  new base T (2/2);  total (42/26)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

GTGGTTAGTGGAACCACCGGTCGCCAGCAGTGCAACGATACCGTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCATGGTAAGACGCGGCTTCTGACTCCAGTAAGG  >  NC_000913/1933696‑1933988
                                                                                                                                                   |                                                                                                                                                 
gtgGTTAGTGGAACCACCGGTCGCCAGCAGTGAAACTATACCGTTCACAACCACTTTCTCATCGATCATCTTACCGATCGGCAACCATTGATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATcgc                                                                                                                                                  >  1:82325/1‑149 (MQ=255)
   gTTAGTGGAACCACCGGTCGCCAGCAGTGCAACTATACCGTTCACCACCACTTTCTCATCGCTCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAACGCAgcgcccc                                                                                                                                               >  2:243257/1‑146 (MQ=255)
   gTTAGTGGAACCACCGGTCGCCAGCAGTGCAACGATACCGTACACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCa                                                                                                                                               >  2:1080539/1‑149 (MQ=255)
       gTGGAACCACCGGTCGCCAGCAGTGCAACGATACCGTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCgg                                                                                                                                           >  2:260612/1‑149 (MQ=255)
        tGGAACCACCGGTCGCCAGCAGTGCAACGATACCGTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGga                                                                                                                                          <  2:1467527/149‑1 (MQ=255)
        tGGAACCACCGGTCGCCAGCAGTGCAACGATACCGTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGga                                                                                                                                          <  1:1535341/149‑1 (MQ=255)
          gAACCACCGGTCGCCAGCAGTGCAACGATACCGTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGgaga                                                                                                                                        >  1:454918/1‑149 (MQ=255)
             ccaccGGTCGCCAGCAGTGCAACGATACCGTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATc                                                                                                                                     <  1:645784/149‑1 (MQ=255)
               accGGTCGCCAGCAGTGCAACGATACCGTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCg                                                                                                                                   >  2:1400436/1‑149 (MQ=255)
               accGGTCGCCAGCAGTGCAACGATACCGTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCg                                                                                                                                   <  2:821460/149‑1 (MQ=255)
                ccGGTCGCCAGCAGTGCAACGATACCGTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCgg                                                                                                                                  >  1:1362207/1‑149 (MQ=255)
                 cGGTCGCCAGCAGTGCAACGATACCGTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGa                                                                                                                                 <  2:894169/149‑1 (MQ=255)
                 cGGTCGCCAGCAGTGCAACGATACCGTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGa                                                                                                                                 <  1:341128/149‑1 (MQ=255)
                  ggTCGCCAGCCGTGCAACGATACCGTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGAt                                                                                                                                <  2:183722/149‑1 (MQ=255)
                        cagcagTGCAACGATACCGTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGCCGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACat                                                                                                                          >  1:90549/1‑148 (MQ=255)
                             gTGCAACGATACCGTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAaagaa                                                                                                                     <  1:120518/149‑1 (MQ=255)
                                cAACGATACCGTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCTCGCAGCGGAGAATCCGGATGAACAAAAGAagag                                                                                                                  >  2:621348/1‑149 (MQ=255)
                                cAACGATACCGTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCTCGCAGCGGAGAATCCGGATGAACAAAAGAagag                                                                                                                  >  2:622278/1‑149 (MQ=255)
                                         cGTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAAc                                                                                                         >  2:882/1‑149 (MQ=255)
                                          gTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACt                                                                                                        <  2:36048/149‑1 (MQ=255)
                                           ttCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTg                                                                                                       >  2:1193385/1‑149 (MQ=255)
                                              accaccacTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCAt                                                                                                    <  1:1080539/149‑1 (MQ=255)
                                                 accacTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCTGGATGAACAAAAGAAGAGCCTGGCAACTGCATccc                                                                                                 >  1:1397320/1‑149 (MQ=255)
                                                 accacTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCTGGATGAACAAAAGAAGAGCCTGGCAACTGCATccc                                                                                                 >  1:1391217/1‑149 (MQ=255)
                                                   cacTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCa                                                                                               >  2:582716/1‑149 (MQ=255)
                                                     atttCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATa                                                                                             <  1:635531/148‑1 (MQ=255)
                                                           cATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTcc                                                                                       >  1:406004/1‑149 (MQ=255)
                                                           cATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTcc                                                                                       >  2:1452967/1‑149 (MQ=255)
                                                            aTCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTcca                                                                                      >  1:550382/1‑149 (MQ=255)
                                                            aTCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTcca                                                                                      >  1:1547467/1‑149 (MQ=255)
                                                             tCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCAt                                                                                                    >  2:112422/1‑134 (MQ=255)
                                                             tCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCAt                                                                                                    <  1:112422/134‑1 (MQ=255)
                                                                atcatcTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTccaccac                                                                                  >  1:478115/1‑149 (MQ=255)
                                                                atcatcTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTccaccac                                                                                  >  2:1362971/1‑149 (MQ=255)
                                                                    tcTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGACTGGCAACTGCATCCCCATAAACTCCACCACCATc                                                                              >  2:431434/1‑149 (MQ=255)
                                                                     cTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCt                                                                             >  2:487712/1‑149 (MQ=255)
                                                                     cTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCt                                                                             <  2:505897/149‑1 (MQ=255)
                                                                     cTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCATGGCAACTGCATCCACATAAACTCCACCACCATCt                                                                             >  2:610611/1‑149 (MQ=255)
                                                                         ccGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCtggtt                                                                         <  1:362368/149‑1 (MQ=255)
                                                                          cGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCtggttg                                                                        >  2:1563759/1‑149 (MQ=255)
                                                                            aTCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGACCAAAAGAAGAGCCTGGCAACTGCATCCCCATACACTCCACCACCATCtggttggt                                                                      >  1:806462/1‑149 (MQ=255)
                                                                              cGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGgtgt                                                                    <  2:696242/149‑1 (MQ=255)
                                                                                          attaccattaccGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGt                                                        >  2:1291036/1‑149 (MQ=255)
                                                                                          attaccattaccGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGt                                                        >  1:1176145/1‑149 (MQ=255)
                                                                                           ttaccattaccGGTCATGCGTGTTTCCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTa                                                       <  2:926852/149‑1 (MQ=255)
                                                                                               cattaccGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGaaa                                                   >  1:718003/1‑149 (MQ=255)
                                                                                               cataaccGGTCATGCGTGAAACCTCACGCGCAGCTGCGGCGGTCAAAGCATAGCGCAGCGGAGAATCCGGATGAAAAAAAGAAGAGACACGCAACTGCATCCCCATAAACTCCAACACCCTCTGCTTGGTGTAGGCAGTACCGTAGaaa                                                   >  2:116074/1‑149 (MQ=255)
                                                                                                attaccGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCTGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAg                                                  <  2:1397320/149‑1 (MQ=255)
                                                                                                attaccGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCTGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAg                                                  <  2:1391217/149‑1 (MQ=255)
                                                                                                 ttaccGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCCCCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGt                                                 >  1:822110/1‑149 (MQ=255)
                                                                                                 ttaccGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCAGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTTGTTGGTGTTGGCAGTACCGTAGAAAGt                                                 >  2:309446/1‑149 (MQ=255)
                                                                                                   accGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTAc                                               <  1:260612/149‑1 (MQ=255)
                                                                                                      ggTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATg                                            >  1:918672/1‑149 (MQ=255)
                                                                                                            gCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCg                                      >  2:31351/1‑149 (MQ=255)
                                                                                                              gtgtAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGgc                                    >  2:1496542/1‑149 (MQ=255)
                                                                                                               tgtAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGgcg                                   <  1:1534596/149‑1 (MQ=255)
                                                                                                                  aaCCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCAt                                <  1:1452967/149‑1 (MQ=255)
                                                                                                                     cTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCATGGt                             >  1:952834/1‑149 (MQ=255)
                                                                                                                     cTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCAATAAACTCCACCACCATCTGGTTGGTGTAGGCAGTACCGTAGAAAGTACATGTTCCCGGAGCATGGt                             >  1:956014/1‑149 (MQ=255)
                                                                                                                            gcAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGACAACTGCATCCCTATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCATGGTAAGAcgc                      >  1:1224772/1‑149 (MQ=255)
                                                                                                                                cTGCGGCGGTCAAAGCATCGCGCAGCGGAGACTCCGGATGAACAAAAGCAGCGCCTGGCACCTGCATCCCCATAATCTCCACCACCATCTGGGAGGTGTTGGCAGTAGCGTAGAAAGTGCATGGTCCCGGCGCATGGTAAGCCGGGGcg                  >  1:270874/1‑148 (MQ=255)
                                                                                                                                   cggcggTCAAAGCATCGGGAAGCGGAGAATCGGGATGAACAAAAGAAGAGACTTGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACGTTATAAAGTACATGTTCCCGTCGCATGTTAATACGCGGCTTCt               <  2:82325/149‑1 (MQ=255)
                                                                                                                                   cggcggTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCATGGTAAGACGCGGCTTCt               <  2:1362207/149‑1 (MQ=255)
                                                                                                                                           aaaGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGATCCTGGCAACTGCATCCCCATAAACTCCACAAACATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCATGGTAAGACGCGGCTTCTGCCTCCAg       >  1:913198/1‑149 (MQ=255)
                                                                                                                                           aaaGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCATGGTAAGACGCGGCTTCTGACTCCAg       >  1:1112979/1‑149 (MQ=255)
                                                                                                                                              gCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCATGGTAAGACGCGGCTTCTGACTCCAGTaa    <  2:1547467/149‑1 (MQ=255)
                                                                                                                                               cATCTCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCATGGTAAGACGCGGCTTCTGACTCCAGTAAg   <  1:622278/149‑1 (MQ=255)
                                                                                                                                               cATCTCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCATGGTAAGACGCGGCTTCTGACTCCAGTAAg   <  1:621348/149‑1 (MQ=255)
                                                                                                                                                aTCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCATGGTAAGACGCGGCTTCTGACTCCAGTAAgg  >  1:1365629/1‑149 (MQ=255)
                                                                                                                                                aTCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCATGGTAAGACGCGGCTTCTGACTCCAGTAAgc  >  1:604268/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTGGTTAGTGGAACCACCGGTCGCCAGCAGTGCAACGATACCGTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCATGGTAAGACGCGGCTTCTGACTCCAGTAAGG  >  NC_000913/1933696‑1933988

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