Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,933,858 C→T 6.1% P249P (CCG→CCA)  edd ← 6‑phosphogluconate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,933,8580CT 6.1% 205.2 / 2.8 68P249P (CCG→CCA) edd6‑phosphogluconate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (39/23);  new base T (2/2);  total (42/26)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-01

ACCGGTCGCCAGCAGTGCAACGATACCGTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCATGGTAAGACGCGGCTTCTGACTCCAGTAAGGCCATGCGGTCCACTTTAC  >  NC_000913/1933711‑1934006
                                                                                                                                                   |                                                                                                                                                    
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:1391217/1‑149 (MQ=255)
                                  accacTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCTGGATGAACAAAAGAAGAGCCTGGCAACTGCATccc                                                                                                                   >  1:1397320/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:1547467/1‑149 (MQ=255)
                                             aTCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTcca                                                                                                        >  1:550382/1‑149 (MQ=255)
                                              tCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCAt                                                                                                                      <  1:112422/134‑1 (MQ=255)
                                              tCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCAt                                                                                                                      >  2:112422/1‑134 (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:505897/149‑1 (MQ=255)
                                                      cTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCt                                                                                               >  2:487712/1‑149 (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:1391217/149‑1 (MQ=255)
                                                                                 attaccGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCTGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAg                                                                    <  2:1397320/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)
                                                                                                                                     cgcAGCGGAGAATCCGGATGAACAAAAGAAGCGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCATGGTAAGACGCGGCTTCTGACTCCAGTAAGGCCAt                >  2:254451/1‑149 (MQ=255)
                                                                                                                                     cgcAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTTGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCATGGTAAGACGCGGCTTCTGACTCCAGTAAGGCCAt                <  1:1193385/149‑1 (MQ=255)
                                                                                                                                          cGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCATGGTAAGACGCGGCTTCTGACTCCAGTAAGGCCATGCGGt           >  2:1202297/1‑149 (MQ=255)
                                                                                                                                          cGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCATGGTAAGACGCGGCTTCTGACTCCAGTAAGGCCATGCGGt           <  2:947770/149‑1 (MQ=255)
                                                                                                                                             agaATCCGGATGAACAAACGAAGAGCCTGGCACCTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCATGGTAAGACGCGGCTTCTGACTCCAGTAAGGCCATGCGGTCCa        <  2:90549/149‑1 (MQ=255)
                                                                                                                                                 tCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTACCGGAGCATGGTAAGACGCGGCTTCTGACTCCAGTAAGGCCTTGCGGTCCACttt    >  2:279195/1‑149 (MQ=255)
                                                                                                                                                   cGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCATGGTAAGACGCGGCTTCTGACTCCAGTAAGGCCATGCGGTCCACTTTAc  <  2:454918/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACCGGTCGCCAGCAGTGCAACGATACCGTTCACCACCACTTTCTCATCGATCATCTTACCGATCGGCATCCATTCATTACCATTACCGGTCATGCGTGTAACCTGACGCGCAGCTGCGGCGGTCAAAGCATCGCGCAGCGGAGAATCCGGATGAACAAAAGAAGAGCCTGGCAACTGCATCCCCATAAACTCCACCACCATCTGGTTGGTGTTGGCAGTACCGTAGAAAGTACATGTTCCCGGCGCATGGTAAGACGCGGCTTCTGACTCCAGTAAGGCCATGCGGTCCACTTTAC  >  NC_000913/1933711‑1934006

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: ‑