Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,409,738 T→C A130A (GCA→GCG yfbT ← putative phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,409,7380TC96.8% 290.8 / ‑1.6 95A130A (GCA→GCGyfbTputative phosphatase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (47/45);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 101 total aligned reads displayed.

CGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCTGCTTTATTC  >  NC_000913/2409648‑2409825
                                                                                          |                                                                                       
cgcAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGgcgc                                                                                       >  2:3366941/1‑93 (MQ=255)
 gcAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGcgcg                                                                                      >  1:1015184/1‑93 (MQ=255)
 gcAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGcgcg                                                                                      >  1:90779/1‑93 (MQ=255)
   aaGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCggg                                                                                    <  2:1085914/93‑1 (MQ=255)
    aGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGa                                                                                   >  2:1084718/1‑93 (MQ=255)
     gCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGaa                                                                                  <  2:2045396/93‑1 (MQ=255)
     gCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGaa                                                                                  <  2:1286617/93‑1 (MQ=255)
      ccccAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAg                                                                                 >  1:1968412/1‑93 (MQ=255)
        ccAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGcc                                                                               <  2:2701290/93‑1 (MQ=255)
        ccAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGcc                                                                               <  1:1835540/93‑1 (MQ=255)
         cagcagCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAATccc                                                                              <  2:3190768/93‑1 (MQ=255)
         cagcagCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGccc                                                                              <  1:2519974/93‑1 (MQ=255)
         cagcagCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGccc                                                                              <  1:517704/93‑1 (MQ=255)
         cagcagCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGccc                                                                              <  1:2058131/93‑1 (MQ=255)
              gcTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTa                                                                         <  2:1009880/93‑1 (MQ=255)
              gcTGCGCGCCTAACAGATACGCATCAGGTTCTGGGTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCtt                                                                         >  1:3484326/1‑92 (MQ=255)
                 gcgcgcCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTAttt                                                                      <  1:3507257/93‑1 (MQ=255)
                   gcgcCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTa                                                                    <  2:1249653/93‑1 (MQ=255)
                    cgcCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTAt                                                                   <  1:1042331/93‑1 (MQ=255)
                    cgcCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTAt                                                                   <  1:2991314/93‑1 (MQ=255)
                      ccTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATgg                                                                 >  1:2756086/1‑93 (MQ=255)
                       cTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGgc                                                                <  1:1303240/93‑1 (MQ=255)
                         aaCAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGgcgc                                                              >  1:4410908/1‑93 (MQ=255)
                          aCAGATACGCATCAGGTTCTGGTTTTCTGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGgcgcg                                                             >  2:4365936/1‑93 (MQ=255)
                          aCAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAGCCCAGCTATTTTATGgcgcg                                                             <  1:745400/93‑1 (MQ=255)
                          aCAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGgcgcg                                                             >  2:3055631/1‑93 (MQ=255)
                          aCAGATACGCATCAGGTTCTGGTTTTCCGCACTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGgcgcg                                                             >  2:3157459/1‑93 (MQ=255)
                           cAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGgcgcgc                                                            <  2:4139338/93‑1 (MQ=255)
                           cAGATACGCATCAGGTTCTGGTTTTCCCCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCTCAGATATTTTATGgcgcgg                                                            >  2:165540/1‑92 (MQ=255)
                            aGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGcgcgct                                                           >  2:655630/1‑93 (MQ=255)
                            aGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGcgcgct                                                           >  1:2755608/1‑93 (MQ=255)
                             gATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGcgcgctc                                                          <  2:670590/93‑1 (MQ=255)
                               tACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGcgcgctcgc                                                        <  2:693398/93‑1 (MQ=255)
                               tACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGcgcgctcgc                                                        >  1:3576432/1‑93 (MQ=255)
                                aCGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGcgcgctcgcg                                                       >  1:1723457/1‑93 (MQ=255)
                                 cGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGcgcgctcgcgc                                                      <  1:4013287/93‑1 (MQ=255)
                                    aTCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTAc                                                   <  2:705112/93‑1 (MQ=255)
                                     tCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTAcc                                                  >  2:1066568/1‑93 (MQ=255)
                                     tCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTAcc                                                  >  2:788072/1‑93 (MQ=255)
                                       aGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCgg                                                >  2:2264283/1‑93 (MQ=255)
                                        ggTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGc                                               >  1:2331152/1‑93 (MQ=255)
                                        ggTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGc                                               >  2:4486074/1‑93 (MQ=255)
                                          ttCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCAt                                             <  1:1043145/93‑1 (MQ=255)
                                          ttCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCAt                                             >  1:3179473/1‑93 (MQ=255)
                                            cTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATgg                                           <  2:3337303/93‑1 (MQ=255)
                                            cTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATgg                                           <  1:3027873/93‑1 (MQ=255)
                                             tGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGa                                          <  2:227302/93‑1 (MQ=255)
                                              ggTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAg                                         <  1:3899693/93‑1 (MQ=255)
                                               gTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGc                                        >  1:1158366/1‑93 (MQ=255)
                                                ttttCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGcc                                       <  1:92399/93‑1 (MQ=255)
                                                ttttCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGcc                                       >  1:3898173/1‑93 (MQ=255)
                                                ttttCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGcc                                       <  2:159473/93‑1 (MQ=255)
                                                   tCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGa                                    <  1:2537973/93‑1 (MQ=255)
                                                    ccGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGaa                                   <  1:3053042/93‑1 (MQ=255)
                                                     cgcgCTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAg                                  >  2:2471969/1‑93 (MQ=255)
                                                      gcgcTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGt                                 >  1:1531998/1‑93 (MQ=255)
                                                      ccgcTTCGCTCGCTCAGCGGTTACAAGCACCTCTGGTGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGt                                 <  1:3897955/92‑1 (MQ=255)
                                                       cgcTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTc                                <  2:2070412/93‑1 (MQ=255)
                                                       cgcTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTc                                <  1:3578629/93‑1 (MQ=255)
                                                        gcTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTca                               >  2:3204376/1‑93 (MQ=255)
                                                        gcTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTca                               >  1:3984314/1‑93 (MQ=255)
                                                        gcTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTca                               >  1:3937758/1‑93 (MQ=255)
                                                         cTTCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTcac                              >  1:2821004/1‑93 (MQ=255)
                                                          ttCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTcaca                             >  2:2271607/1‑93 (MQ=255)
                                                           tCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACaa                            <  1:4677524/93‑1 (MQ=255)
                                                           tCACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACaa                            <  1:4133800/93‑1 (MQ=255)
                                                           tCACTCGCGCAGCGGTTACAAACACCTTTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACaa                            <  1:1126334/93‑1 (MQ=255)
                                                            cACTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAAt                           <  1:3839777/93‑1 (MQ=255)
                                                             aCTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATg                          >  2:655866/1‑93 (MQ=255)
                                                              cTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATgg                         >  1:2596130/1‑93 (MQ=255)
                                                              cTCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATgg                         >  2:1780745/1‑93 (MQ=255)
                                                               tCGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGc                        >  1:3710882/1‑93 (MQ=255)
                                                                cGCTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGcc                       <  2:2194529/93‑1 (MQ=255)
                                                                 gcgcAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGccc                      <  2:4233702/90‑1 (MQ=255)
                                                                  cTCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCa                     <  1:3592139/93‑1 (MQ=255)
                                                                  cTCAGCGGTTACAAACACCTCTGGCGCGGGAGCCCAGCTATTTTATTGCGGGCTCGCGCTACCGGCATGGAGCCAGAAGGCACAATGGCCCAc                    >  2:234192/1‑93 (MQ=255)
                                                                   tCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCAc                    >  1:3725069/1‑93 (MQ=255)
                                                                   tCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCAc                    <  1:1098645/93‑1 (MQ=255)
                                                                   tCAGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAACTCACAATGGCTCAc                    >  2:1785210/1‑93 (MQ=255)
                                                                     aGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGTCATGGAGCCAGAAGTCACAATGGCCCACgg                  >  1:774958/1‑93 (MQ=255)
                                                                     aGCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACgg                  >  1:133657/1‑93 (MQ=255)
                                                                      gCGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGa                 >  2:2479414/1‑93 (MQ=255)
                                                                       cGGTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGaa                >  1:4312142/1‑93 (MQ=255)
                                                                        ggTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAAt               <  1:910843/93‑1 (MQ=255)
                                                                        ggTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAAt               >  1:574911/1‑93 (MQ=255)
                                                                        ggTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAAt               <  1:3532237/93‑1 (MQ=255)
                                                                         gTTACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATa              >  2:4096355/1‑93 (MQ=255)
                                                                          ttACAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATAc             >  1:220842/1‑93 (MQ=255)
                                                                           tACAAACACCCCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATAcc            <  1:3943083/93‑1 (MQ=255)
                                                                            tcAAACTCCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCTCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCt           <  2:1632282/92‑1 (MQ=255)
                                                                            aCAAACACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCt           >  1:2521516/1‑93 (MQ=255)
                                                                               aaCACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCTGCt        >  2:2496861/1‑93 (MQ=255)
                                                                               aaCACCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCTGCt        >  2:4520779/1‑93 (MQ=255)
                                                                                acacCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCTGCtt       >  2:2242441/1‑93 (MQ=255)
                                                                                acacCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCTGCtt       >  2:3476278/1‑93 (MQ=255)
                                                                                acacCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCTGCtt       <  1:2203199/93‑1 (MQ=255)
                                                                                acacCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCTGCtt       <  1:1531191/93‑1 (MQ=255)
                                                                                  acCTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCTGCTTTa     >  1:2639204/1‑93 (MQ=255)
                                                                                   ccTCTGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCTGCTTTAt    <  2:2454426/93‑1 (MQ=255)
                                                                                     tctGGCGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCTGCTTTATTc  >  1:1225031/1‑93 (MQ=255)
                                                                                          |                                                                                       
CGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAGCCCAGCTATTTTATGGCGCGCTCGCGCTACCGGCATGGAGCCAGAAGTCACAATGGCCCACGGAATACCTGCTTTATTC  >  NC_000913/2409648‑2409825

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