Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,666,466 T→G 23.8% Y384* (TAT→TAG clcB → H(+)/Cl(‑) exchange transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,666,4660TG23.8% 213.3 / 73.7 117Y384* (TAT→TAGclcBH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base T (46/42);  new base G (15/13);  total (62/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

CTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCA  >  NC_000913/1666368‑1666564
                                                                                                  |                                                                                                  
cTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAt                                                                                                    >  2:2714985/1‑99 (MQ=255)
cTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATc                                                                                                   <  1:2094724/100‑1 (MQ=255)
cTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATc                                                                                                   >  2:2342349/1‑100 (MQ=255)
  tttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAg                                                                                                 <  2:82914/100‑1 (MQ=255)
   ttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGc                                                                                                <  2:1618519/100‑1 (MQ=255)
    tACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCt                                                                                               <  2:343610/100‑1 (MQ=255)
      cTCGGCTTGACCGGGATGGCGACACTATCGGCGGCGACCACGCACGCGACGAGTCTGTCGAGGTTGATGAAATGTGAAATGACCGGCGAGGAACAggt                                                                                               >  1:13743/1‑96 (MQ=255)
      cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCCCGCGCCGATTATGTCGACGTTGAGGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                             >  1:1405158/1‑100 (MQ=255)
      cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                             >  2:2310925/1‑100 (MQ=255)
      cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                             >  2:1022743/1‑100 (MQ=255)
      cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                             >  2:147738/1‑100 (MQ=255)
        cGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTc                                                                                           >  1:824546/1‑100 (MQ=255)
           aTTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTcccc                                                                                        <  1:236702/100‑1 (MQ=255)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGAGGTTGATGATATGTGAAATGACCGAGGAGTAGCTGATACTCCCTgg                                                                                      <  1:367655/100‑1 (MQ=255)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCgg                                                                                      <  2:451756/100‑1 (MQ=255)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCgg                                                                                      <  2:2119664/100‑1 (MQ=255)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATGCTCCCCgg                                                                                      <  1:593049/100‑1 (MQ=255)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                      <  2:2436342/100‑1 (MQ=39)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                      <  2:1792048/100‑1 (MQ=39)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                      <  1:1263877/100‑1 (MQ=39)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACACACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCgg                                                                                      <  2:2637407/100‑1 (MQ=255)
              gACCTGGATTGCGACACTGTTGGCCGCAACCACTCACGCTCCGATTATGTCGACGTTGCTTATATTTGCAATGACCGGGGGGTAGCTGATACTCCCCgg                                                                                      <  1:1914377/99‑1 (MQ=255)
              gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGt                                                                                     >  2:1644683/1‑100 (MQ=255)
              gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGt                                                                                     >  2:515344/1‑100 (MQ=255)
              gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGt                                                                                     <  2:2723823/100‑1 (MQ=39)
              gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGAGGACATGTGAAATGACCGGGGAGTATCAGCTACTCCCCgg                                                                                      >  1:938444/1‑99 (MQ=255)
                 cGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTa                                                                                  >  2:1245168/1‑100 (MQ=255)
                  gggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTat                                                                                 >  2:465134/1‑100 (MQ=255)
                       ggCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgatt                                                                            <  2:1640739/100‑1 (MQ=255)
                       ggCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgatt                                                                            >  2:2720506/1‑100 (MQ=255)
                       ggCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGtcatttcac                                                                             <  2:659457/99‑9 (MQ=37)
                        gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGtcatttcacat                                                                           <  2:1847746/100‑11 (MQ=25)
                        ccgatcttGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgattg                                                                           <  2:902675/93‑1 (MQ=255)
                         cGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGc                                                                          >  2:2511086/1‑100 (MQ=255)
                         cGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGc                                                                          >  1:2857409/1‑100 (MQ=255)
                          gACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGcc                                                                         <  2:2568576/100‑1 (MQ=255)
                          gACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTCGCTGATATGCGAAATGACCGGGGAGCAGCCGCTACTCCCCGGTTTATGGATTGc                                                                          >  1:2486350/1‑99 (MQ=255)
                              cTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCg                                                                     >  1:1045253/1‑100 (MQ=255)
                              cTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCg                                                                     <  2:854934/100‑1 (MQ=255)
                               tGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgattgagatcgg                                                                    >  1:902672/1‑93 (MQ=255)
                               tGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGt                                                                    <  1:225793/100‑1 (MQ=255)
                                gTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTa                                                                   <  2:2588607/100‑1 (MQ=255)
                                  tGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAt                                                                 >  2:398685/1‑100 (MQ=255)
                                  tGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAt                                                                 <  1:2214614/100‑1 (MQ=255)
                                   ggcggcGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAtt                                                                >  2:2076137/1‑100 (MQ=255)
                                    gcggcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTg                                                               <  1:1293521/100‑1 (MQ=255)
                                    gcggcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTg                                                               <  2:231423/100‑1 (MQ=255)
                                    gcggcgACCACGCACGCGCCGATTATGTCGACGTCGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTg                                                               <  1:1220388/100‑1 (MQ=255)
                                      ggcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGCAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCg                                                             >  1:82449/1‑100 (MQ=255)
                                      ggcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATTACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCGGCGTAATTGCg                                                             <  1:488555/100‑1 (MQ=255)
                                          aCCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCgg                                                         <  1:2216573/100‑1 (MQ=255)
                                             acgcacgcGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTaa                                                      >  1:2332169/1‑100 (MQ=255)
                                              cgcacgcGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAAt                                                     >  2:303631/1‑100 (MQ=255)
                                              cgcacgcGCCGATTATGTCAAAGTTGATGATGGGTGAAATGGCGGGGGAGTATCAGCTACTCCCAGTTTTATTGATTGCCTGCGTAAGTGCGACGGTAAt                                                     >  1:1763492/1‑100 (MQ=255)
                                                cacgcGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAAttt                                                   >  2:1992651/1‑100 (MQ=255)
                                                  cgcgCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcg                                                 <  2:2180270/100‑1 (MQ=255)
                                                  cgcgCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcg                                                 >  2:1137258/1‑100 (MQ=255)
                                                    cgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTGATTTcgcg                                               >  1:2057842/1‑100 (MQ=255)
                                                    cgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                               <  2:2727804/100‑1 (MQ=255)
                                                     gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCgg                                              <  2:2450346/100‑1 (MQ=255)
                                                       cGATTCTGTCGCCATTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGa                                             >  1:1804261/1‑99 (MQ=255)
                                                       cGATTATGTCGACGTTGATGATATGTGAAATTACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                               <  2:2321063/97‑1 (MQ=255)
                                                        gATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACg                                           >  2:1647395/1‑100 (MQ=255)
                                                         aTTATGTCGACGTTGATGATATGTGAAATGACCGGGGCGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCAGACGt                                          <  1:1967595/100‑1 (MQ=255)
                                                          ttATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGtt                                         >  1:933433/1‑100 (MQ=255)
                                                              gTCGACGTTGATGATATGTGAAATGACTGGGGAGTATCAGCTACTCCCCGGTTTATTTATTGCCTGCGTAATTGCGTCGGTACTTTCGCGGACGTTac                                       >  1:1897579/1‑98 (MQ=255)
                                                              gTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTacac                                     >  2:322591/1‑100 (MQ=255)
                                                              gTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTacac                                     <  1:2610226/100‑1 (MQ=255)
                                                                       gatgatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGActct                            >  2:1901643/1‑100 (MQ=255)
                                                                        atgatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtcta                           <  2:2733641/100‑1 (MQ=255)
                                                                          gatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatc                         <  2:635951/100‑1 (MQ=255)
                                                                          gatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatc                         <  2:2158169/100‑1 (MQ=255)
                                                                          gatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatc                         <  1:1881310/100‑1 (MQ=255)
                                                                           agatGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatct                        >  2:1978384/3‑100 (MQ=255)
                                                                             atGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAc                      >  1:1691415/1‑100 (MQ=255)
                                                                              tgtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAcc                     >  1:1163698/1‑100 (MQ=255)
                                                                               gtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCg                    >  2:2172795/1‑100 (MQ=255)
                                                                               atcaataaaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCg                    >  1:49703/9‑100 (MQ=37)
                                                                                  aaaTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcn                  <  1:2379296/99‑2 (MQ=255)
                                                                                     aaaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGca              >  2:1340662/3‑100 (MQ=38)
                                                                                       aCCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGcaca            >  1:134825/1‑100 (MQ=255)
                                                                                       aCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGcaca            >  2:1716554/1‑100 (MQ=39)
                                                                                        ccGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac           >  1:1446807/1‑100 (MQ=255)
                                                                                        ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAgatcgg           >  2:1278328/1‑95 (MQ=38)
                                                                                        ccGGGGAGTAGCTGATACTCCCCGGTTTAGTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac           >  1:175025/1‑100 (MQ=39)
                                                                                         cGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACt          >  2:405258/1‑100 (MQ=255)
                                                                                         cGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACt          >  2:442896/1‑100 (MQ=255)
                                                                                         cGGGGAGTAGCTGATACTCAAAGGATTATTGATTGCCTGCGTAAGTGCGGCGGTAATTTCGCGGACGGTCAACCGGGACACTATGTACCGccgccccact          >  1:1277787/1‑92 (MQ=255)
                                                                                          ggggAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACt          >  1:1741733/1‑99 (MQ=255)
                                                                                          ggggAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg         >  1:2840953/1‑100 (MQ=255)
                                                                                          ggggAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg         >  1:1437485/1‑100 (MQ=255)
                                                                                          ggggAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg         >  2:1182521/1‑100 (MQ=255)
                                                                                          ggggAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg         >  2:161778/1‑100 (MQ=255)
                                                                                          ggggAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg         >  1:1100350/1‑100 (MQ=255)
                                                                                           gggAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTgc        <  2:2320324/100‑1 (MQ=255)
                                                                                            ggAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTTCGTCGGTAATTTCGCGGACGTTCCACCGTGACTCTGTCTACCGCCAGCacacgg         >  1:548195/1‑96 (MQ=255)
                                                                                              aGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgca     <  2:834340/100‑1 (MQ=255)
                                                                                               gTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag    >  2:866105/1‑100 (MQ=255)
                                                                                               gTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag    <  1:1933108/100‑1 (MQ=255)
                                                                                                 aTCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcagca  <  2:2226272/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
CTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCA  >  NC_000913/1666368‑1666564

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