Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,096,861 A→C 7.7% T61P (ACC→CCC)  rdgB → dITP/XTP pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,096,8610AC 7.7% 236.6 / 2.2 106T61P (ACC→CCC) rdgBdITP/XTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (55/41);  new base C (2/8);  total (57/49)
Fisher's exact test for biased strand distribution p-value = 4.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

TCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCG  >  NC_000913/3096712‑3097010
                                                                                                                                                     |                                                                                                                                                     
tCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCAGCAAAAGTGa                                                                                                                                                       >  1:1038240/1‑150 (MQ=255)
   gTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCg                                                                                                                                                    <  1:1057543/150‑1 (MQ=255)
     aaaGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCt                                                                                                                                                  <  1:1170404/150‑1 (MQ=255)
     aaaGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCt                                                                                                                                                  >  1:915858/1‑150 (MQ=255)
          gCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCCCGCCATGCGGCAAAAGTGACCGCTTTAcc                                                                                                                                             >  2:1850009/1‑150 (MQ=255)
          gCGTGAGCGGGCGTCGCTGCTTATCGACTTCGGGCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAACCCGGCCTGACCTTTATCGAATACGCGATTCTGAAAGCGCGCCATGCGGGAAAAGTGACCGCTTTAcc                                                                                                                                             <  2:1261857/150‑1 (MQ=255)
            gTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCgg                                                                                                                                           <  2:1780973/150‑1 (MQ=255)
                 cTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAAtt                                                                                                                                      >  1:878294/1‑150 (MQ=255)
                  tGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTg                                                                                                                                     <  2:70000/150‑1 (MQ=255)
                  tGGCGTCGCTGCTTAGCGACGTCGGGCTTGCTATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCGGAAGCAACCGGCCTGACCTGTAACGAAAACGCGATTCTGAAAGCGCGCCAGGCGGCAAAAGTGCCCGCGTTACCGGCAAtag                                                                                                                                     >  2:119862/1‑148 (MQ=255)
                     cGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCcg                                                                                                                                  >  1:1769837/1‑150 (MQ=255)
                     cGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCcg                                                                                                                                  <  1:781646/150‑1 (MQ=255)
                       tctttgctTTGCGACTTCTGTCTTTATATCGTGGCCCAAAGAGACCTCGGCGTTGATTCGGCTGAAGAAACCGGCCTGACCTTTAGCGAAAACGCGATTCTGAAAGCGCGCCTTGCGGCAAAAGTGACCGCTTTACCGGCAATTGCcgac                                                                                                                                <  1:268405/146‑1 (MQ=255)
                       tCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGCCCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCcgac                                                                                                                                >  1:847697/1‑150 (MQ=255)
                        cGCTGCTTAGCGACTTAGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCcgacg                                                                                                                               <  2:1139421/150‑1 (MQ=255)
                             ctTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGActct                                                                                                                          <  2:227767/150‑1 (MQ=255)
                             ctTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACATCGGCGTTGATTCCTCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGCCAATTGCCGACGActct                                                                                                                          >  1:572220/1‑150 (MQ=255)
                               tAGCGACTTCGGTCTTTATATCGGGGCCAAAACAGACCTCGGCGTTGATGCCGCTGAAGAACCCGGCCTGACCTTTATCGAAACCGCGATTCTGAAAGCGCGCCTTGCGGAAAAAGTGACCGCTTTACCGGCATTTGCCGACGACtctgg                                                                                                                        <  1:1859005/150‑1 (MQ=255)
                                 gCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCATTATCGAAAACGCGTTTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACtctggtc                                                                                                                      >  1:247934/1‑150 (MQ=255)
                                     cTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGc                                                                                                                  <  1:49696/150‑1 (MQ=255)
                                      tgcGGGCTTGATTTCGGGGCCCAAACCTACCGCGGCGTTGATTCCGCAGAAGAAACCGGCCTGACCTTAATCGAAAACGCTATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCg                                                                                                                 <  2:336939/148‑1 (MQ=255)
                                       tCGGTCTTGATATCGGGGCCCAAACAGCCCTCGGCGTTGATTCCGCGGATGAAACCTGCCTGACCTTAATCTAAACCGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCgg                                                                                                                <  1:1226001/150‑1 (MQ=255)
                                        cGGTCTTGATATCGTGGCCCAAACAGACATCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACAGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGt                                                                                                               >  2:1303587/1‑150 (MQ=255)
                                         ggTCTTGATATCGGGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGACACCGGCCTGACCTTTATCGAAACCGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTa                                                                                                              <  1:1311407/150‑1 (MQ=255)
                                           tCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGa                                                                                                            >  2:915297/1‑150 (MQ=255)
                                             ttGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGCTTCTGAACGCGCGCCATGCGGCAAACGTGCCCGCGTTACCGGCAATTGCAGACGACTCTGGTCTGGCGGTAGAtg                                                                                                          >  2:308555/1‑150 (MQ=255)
                                             ttGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGAtg                                                                                                          <  2:1174482/150‑1 (MQ=255)
                                              tGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACAGGCCTGACATTTATCGAAAACTCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGAtgt                                                                                                         >  2:470270/1‑150 (MQ=255)
                                                  atCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCtt                                                                                                     >  2:784912/1‑150 (MQ=255)
                                                  atCGTGGCCCAAACAGACCTCGGAGTTGATTCCGCTGAAGAAACCTGACTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGAGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGACGTAGTGGTGCtt                                                                                                     >  1:172390/1‑150 (MQ=255)
                                                   tCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAACGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTg                                                                                                    <  1:837687/150‑1 (MQ=255)
                                                    cggcgCCCAAACCGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAACCGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTgg                                                                                                   <  1:841859/146‑1 (MQ=255)
                                                     gggccccAACCACCCCTCGGCGTTTATTCCGCTAAACAACCCGGCCTCACCTATATCAAAAACGCGATTCTTAACGCGCGCCATGCGCAAAACGGGCCCGCTTAACCGGCAATTCCCGACGACTCTGGTCTGGCGGTAGATGTGCTTggc                                                                                                  <  2:1754061/146‑1 (MQ=255)
                                                          ccAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGgcgc                                                                                             >  1:167343/1‑150 (MQ=255)
                                                            aaaCAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGgcgcgc                                                                                           >  2:1772750/1‑150 (MQ=255)
                                                                    cTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATTGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCAtt                                                                                   >  1:809726/1‑150 (MQ=255)
                                                                     tCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCAttt                                                                                  <  2:220771/150‑1 (MQ=255)
                                                                     ccggcggTGATGCCGCTGACGAAACCGGCGTGACCTTTAGCGAAACCGCGATGCTTAACGGGCGCCTTGCGGCAAATGTGACCGCTTTCCGGGCATTTGCGGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCAttt                                                                                  <  2:1704288/149‑1 (MQ=255)
                                                                      cGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTa                                                                                 >  2:1732725/1‑150 (MQ=255)
                                                                       ggCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAACGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTAc                                                                                <  1:1183938/150‑1 (MQ=255)
                                                                         ggtggATTCCGCTGACGAAACCGCCCGGCCCTTTCACGAACCCGCGATGCTGAATCGGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTc                                                                              <  1:358948/146‑1 (MQ=255)
                                                                            tGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGTTCTGGCGGTAGATGTGCTTAGCGGCGCGCCGGGCATTTAATCcgc                                                                           >  2:1393479/1‑150 (MQ=255)
                                                                              cttCCGCTGAACAAACCGCCCTGCCCTTTATCGAAACCGCGATTCTGAAAGCGCGCCATGCGGCAAATGGGACCGCTTTCCCGGCAATTCCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCcgcgc                                                                         <  2:915858/149‑1 (MQ=255)
                                                                              aTTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCcgcgc                                                                         >  1:38654/1‑150 (MQ=255)
                                                                               ttCCGCTTAACACCCCGCCCTGACCTTTATCGAAACCGCGAGACTGAAAGCGCGCCATGCGGCAAAACTGACCGCTTTACCGCCAATTGACGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCcgcgcg                                                                        <  2:362361/150‑1 (MQ=255)
                                                                               ttCCGCTGAAGAACCCGGCCTGACCTTTATCCAAAACGCGATTCTGAACGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCcgcgcg                                                                        <  1:487883/150‑1 (MQ=255)
                                                                                  cGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGCTCAGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTa                                                                     >  1:1370694/1‑150 (MQ=255)
                                                                                   gCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCGGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTAt                                                                    >  1:833886/1‑150 (MQ=255)
                                                                                      taataacccGGCCTGACCTTACTCTAAAACGCGATTCTCAAAGCGCGCCATGCGGCAAAAGTGACCGCTTACCCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTcc                                                                 <  2:813383/143‑1 (MQ=255)
                                                                                      caagaaACCGCCCTGACCTTTATCGAAACCGCGATTCTCAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTcc                                                                 <  1:900710/149‑1 (MQ=255)
                                                                                       aagaaACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCg                                                                >  1:1132374/1‑150 (MQ=255)
                                                                                        agaaACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGCCTCTGGTCTGGCGGGAGAGGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCgg                                                               >  2:1792425/1‑150 (MQ=255)
                                                                                        agaaACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCgg                                                               >  1:159867/1‑150 (MQ=255)
                                                                                        agaaACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGGCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCgg                                                               >  2:854885/1‑150 (MQ=255)
                                                                                        agaaACCGGCCTGACCTTTATCGAAAACGCGATTATGAAAGAGCGCCATCCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCGGGTCTGTCGTTAGAAGTGCTTGGCGGCGCTCCGGGCATTTACTCCGCGCGTTATTCCgg                                                               >  1:1256109/1‑150 (MQ=255)
                                                                                          aacccGGCCTGACCTTTTTCCAAAACGCGATACTGAAAGCGCGCCAGGCGGAAAACGTGCCCGCTTTACGGCCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTg                                                             <  2:785982/147‑1 (MQ=255)
                                                                                          aaaCCGGCCTGACCTTTATCGAAACCGCGATTCGGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTg                                                             <  2:1108176/150‑1 (MQ=255)
                                                                                          aaaCCGGCCTGACCTTTAACAAAACCGCGATTCTCAACCCGCGCCATGCGGCAAAAGGGCCCGCTTTACCGCCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTg                                                             <  2:451350/150‑1 (MQ=255)
                                                                                          aaaCCGGCCTGACCTTACTCGACCCCGCGATTCTGACCGCGCGCCATGCGGCAACAGGGCCCGCTTTCCCGGCACGAGCCGACGACGCCGGGCTGGCGGGAGCTGTGCGTGGCGGCGCGCCGGGCATTTCCTCCGCGCGTTATTCCGGTg                                                             <  1:119862/150‑1 (MQ=255)
                                                                                             ccGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAg                                                          >  1:1524752/1‑150 (MQ=255)
                                                                                             ccGGCCTGACCATACTCGAAAACGCGATTCTCAAAGCGCGCCATGCGGCAAAAGTGACCGCGTTACCGGCTTTTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGTCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAg                                                          <  1:124661/150‑1 (MQ=255)
                                                                                                gCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGACATGCGGCAAAAGTGACCGCTTTACCGCCAATTGCCGACGACTCTGGACTGGCGGTAGATGTGCTTGGCGGCGAGCCGGGCATTTACTCCGCTCGTTATTACGGTGAAGAcg                                                       >  1:1074245/1‑150 (MQ=255)
                                                                                                 ccTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGTGGCAAAAGTGACCGCTTTACCAGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGAcgc                                                      >  1:1621495/1‑150 (MQ=255)
                                                                                                  cTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGAcgcg                                                     >  1:1862773/1‑150 (MQ=255)
                                                                                                  cTGACCTTAAACGACAACCCCATTCAGAAACCCTGCCATGCGCCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACACTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGGTATTCCGGTGAAGAcgcg                                                     <  2:128949/150‑1 (MQ=255)
                                                                                                   tGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGATAAGACGCGa                                                    >  1:1751573/1‑150 (MQ=255)
                                                                                                    taCCTTCAACTAAACCCCGTTTCTGAATCCGCGCAATGCGGCAAAAGTCCCGGCTTTACGGGCATTGGCCCACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGAc                                                   <  1:1822955/149‑1 (MQ=255)
                                                                                                    gACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGAAAAAGTGACCGCTATACCGGCAATTGCCGACGACTCTGGTATGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCAGTTACTCCGCGCGTTATTCCGGTGAAGACGCGAc                                                   >  1:1673698/1‑150 (MQ=255)
                                                                                                    gACATTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGCCCGAGTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGAc                                                   <  1:96593/150‑1 (MQ=255)
                                                                                                      ccTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCg                                                 >  2:163035/1‑150 (MQ=255)
                                                                                                       cTTTATCGAAAACCCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGa                                                <  1:974040/150‑1 (MQ=255)
                                                                                                         ttATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCAGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATc                                              >  2:365572/1‑150 (MQ=255)
                                                                                                          tATCGAACCCGCGATTCTGAACGCCCGCCATGCGGCAAAAGTGACCGCTATACCGGCAATTGCCGACGACTCTGGTCTTGCTGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCa                                             <  2:1309279/150‑1 (MQ=255)
                                                                                                          tATCGAAACCGCGATTCTGAACGCGCGCCATGCGGCAAAAGTCACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGTCGCGACCGATCa                                             <  2:1990440/150‑1 (MQ=255)
                                                                                                          tATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCa                                             >  2:643758/1‑150 (MQ=255)
                                                                                                          tATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCa                                             >  2:1196034/1‑150 (MQ=255)
                                                                                                             cGAAACCGCGATTGTTAAAGCGCGCCTTGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGTTCGGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCaaaa                                          <  1:1283242/150‑1 (MQ=255)
                                                                                                              gAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAg                                         >  2:1701096/1‑150 (MQ=255)
                                                                                                              caaaaCCCTATTGTCAAAGCGCGCCATGCGGAAAAATTTCCCGTTTAACCGCCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAg                                         <  2:1482175/149‑1 (MQ=255)
                                                                                                                aaaCGCGATTCTGTAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGTCATTTACTCCGCGCGTTATTACGGTGAAGACGCGACCGATCAAAAGaa                                       >  2:1428643/1‑150 (MQ=255)
                                                                                                                     cgATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGCTCAAAAGAATCTGc                                  >  1:1206914/1‑150 (MQ=255)
                                                                                                                     cgATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGc                                  >  1:748091/1‑150 (MQ=255)
                                                                                                                      gATTCTTAACGGGCGCCTTGGGGAAAATGTCACCGTTTTACCGGAATTTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCa                                 <  2:921183/150‑1 (MQ=255)
                                                                                                                       aTTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGAGCGATCAAAAGAATCTGCaa                                >  1:1787701/1‑150 (MQ=255)
                                                                                                                          cTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCaaaaa                             >  2:434333/1‑150 (MQ=255)
                                                                                                                           tGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCAGTTACTCCGCGCGTTATTCCGGTGAAGACGCGCCCGATCAAAAGAATCTGCAAAAAc                            >  1:1980956/1‑150 (MQ=255)
                                                                                                                               aGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAActgct                        >  1:1847284/1‑150 (MQ=255)
                                                                                                                                 cgcgcCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCCGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTgg                      >  1:1954706/1‑150 (MQ=255)
                                                                                                                                  gcgcCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGa                     <  1:1335288/150‑1 (MQ=255)
                                                                                                                                    gcCATGCGGCAAAAGTGACCGCTTTACCGGGAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGaaa                   >  1:701860/1‑150 (MQ=255)
                                                                                                                                     ccAGGCGGCAAAAGCGCCCCCTTAACCGGCAATCCCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAAc                  <  2:734637/150‑1 (MQ=255)
                                                                                                                                      cATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGTCTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACa                 >  2:366108/1‑150 (MQ=255)
                                                                                                                                      cATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACa                 <  1:1962675/150‑1 (MQ=255)
                                                                                                                                            gCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGaaa           <  2:618039/150‑1 (MQ=255)
                                                                                                                                            gCAAAAGTGACCGCTTAACCGGAAATTCCCGACGACTCTGGTCGGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTAACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGaaa           <  2:1756768/150‑1 (MQ=255)
                                                                                                                                             cAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCATAAACTGCGGGAAACAATGAAAg          >  2:1313396/1‑150 (MQ=255)
                                                                                                                                             cAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAg          >  2:1437649/1‑150 (MQ=255)
                                                                                                                                             cAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCC‑‑AGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGAc        >  2:2010383/1‑150 (MQ=255)
                                                                                                                                                aaGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACg       >  2:727719/1‑150 (MQ=255)
                                                                                                                                                     aCCGTTTTTCAGGTATTTGCCGTCGACTCTGGTCTGGTGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCTCTTTTTTCCGGTTAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCg  <  1:1850009/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCG  >  NC_000913/3096712‑3097010

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