Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,062,371 A→T 15.1% N42I (AAT→ATT)  yihN → MFS transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,062,3710AT15.1% 581.9 / 85.8 231N42I (AAT→ATT) yihNMFS transporter family protein
Reads supporting (aligned to +/- strand):  ref base A (107/89);  new base T (17/18);  total (124/107)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-02

Only 100 of 236 total aligned reads displayed.

TATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACA  >  NC_000913/4062272‑4062469
                                                                                                   |                                                                                                  
 aTTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAAt                                                                                                   <  1:1384366/100‑1 (MQ=255)
   ttAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATgg                                                                                                 <  1:3508335/100‑1 (MQ=255)
    tAGGCTATTAACACTGTGTGGCGGTACAATGTAGAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTCTAACAATGGt                                                                                                >  1:2900711/1‑100 (MQ=255)
         tATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGCTGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACGAATGGTCaagt                                                                                           >  1:2532572/1‑98 (MQ=255)
         tATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAt                                                                                           <  2:1652884/100‑1 (MQ=255)
         tATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAt                                                                                           <  1:36362/100‑1 (MQ=255)
          aTTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAtg                                                                                          <  1:3384636/100‑2 (MQ=255)
          aTTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAtg                                                                                          <  1:2249839/100‑2 (MQ=255)
          aTTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAAtt                                                                                          >  2:1840359/1‑100 (MQ=255)
          aTTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAAtt                                                                                          >  2:3422832/1‑100 (MQ=255)
           ttAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAtgg                                                                                         <  1:3695770/100‑3 (MQ=255)
           ttAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAtgg                                                                                         <  2:3920235/100‑3 (MQ=255)
            tAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAtgga                                                                                        <  2:679930/100‑4 (MQ=255)
            tAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAtgga                                                                                        <  1:379629/100‑4 (MQ=255)
             aaCACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAtggaa                                                                                       <  1:1242335/100‑5 (MQ=255)
             aaCACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGt                                                                                       >  1:2953556/1‑100 (MQ=255)
             aaCACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGt                                                                                       <  2:4040214/100‑1 (MQ=255)
             aaCACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGt                                                                                       <  2:3412652/100‑1 (MQ=255)
             aaCACTGTGTGGCGGTACAATTTATAAATTACCGTCACTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGt                                                                                       >  2:3108065/1‑100 (MQ=255)
               cacTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAtggaaat                                                                                     <  2:2323437/100‑7 (MQ=255)
               cacTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTaa                                                                                     >  2:1613274/1‑100 (MQ=255)
                acTGTGTGGCGGTACAATGTATAAATTACCGTCGCTGAAAGCTGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAAt                                                                                    >  1:1847318/1‑100 (MQ=255)
                 cTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATg                                                                                   >  2:3144700/1‑100 (MQ=255)
                  tgtgtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGc                                                                                  <  1:2732738/100‑1 (MQ=255)
                  tgtgtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGc                                                                                  <  2:3209937/100‑1 (MQ=255)
                   gtgtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCt                                                                                 >  1:2903100/1‑100 (MQ=255)
                   gtgtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCt                                                                                 <  2:3579735/100‑1 (MQ=255)
                   atctgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCt                                                                                 <  1:2081304/97‑1 (MQ=255)
                    tgtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTa                                                                                <  2:4171829/100‑1 (MQ=255)
                     gtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTAt                                                                               >  1:4149831/1‑100 (MQ=255)
                     gtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTAt                                                                               <  2:1941719/100‑1 (MQ=255)
                      tgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATg                                                                              >  1:155166/1‑100 (MQ=255)
                       ggCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTagat                                                                             >  2:2081312/1‑97 (MQ=255)
                       ggCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGt                                                                             >  2:1802640/1‑100 (MQ=255)
                       ggCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGt                                                                             >  1:696967/1‑100 (MQ=255)
                        gCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTc                                                                            <  2:395251/100‑1 (MQ=255)
                         cGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCg                                                                           >  2:522553/1‑100 (MQ=255)
                          ggTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCgg                                                                          >  2:2074887/1‑100 (MQ=255)
                            tACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTa                                                                        <  2:1393044/100‑1 (MQ=255)
                            tACAATGTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGATCAAATTGGTAATGCTATGTCGGTa                                                                        >  1:1712480/1‑100 (MQ=255)
                             aCAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTaa                                                                       <  1:583977/100‑1 (MQ=255)
                              cAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTATTGCTATGTCGGTaaa                                                                      >  2:4328066/1‑100 (MQ=255)
                               aaTTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAAc                                                                     >  2:3823484/1‑100 (MQ=255)
                               aaTTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAAc                                                                     <  2:2792697/100‑1 (MQ=255)
                                     taAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCAttt                                                               <  1:964175/100‑1 (MQ=255)
                                         ttACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTca                                                           <  2:277228/100‑1 (MQ=255)
                                          tACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTccc                                                          <  2:2979693/100‑3 (MQ=255)
                                          tACCGTCGATGAAAGCTGCGTTTTATATAACGATGCAGGACTCTTTCCATTTGACCATTGGTCAAATTGGTAATGCTTTGTCGGTAAACGCATTTGTcac                                                          <  1:2561413/100‑1 (MQ=255)
                                            ccGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTcacca                                                        >  1:2408740/1‑100 (MQ=255)
                                              gTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACcaca                                                      >  2:630883/1‑100 (MQ=255)
                                                cGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTAGTAATGCTAAGTCGGTAAACTCATTTGTCACcaca                                                      >  2:457374/1‑98 (MQ=255)
                                                  cTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTgg                                                  >  2:1444524/1‑100 (MQ=255)
                                                    gAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGc                                                >  1:628017/1‑100 (MQ=255)
                                                        gATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCtttt                                            <  1:4153525/100‑1 (MQ=255)
                                                           gCGTTTTATATCCCGATTCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAACGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTCTTTTc                                         <  1:1452977/100‑1 (MQ=255)
                                                           gCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTc                                         <  1:180318/100‑1 (MQ=255)
                                                            cGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCt                                        <  2:1226680/100‑1 (MQ=255)
                                                            cGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCt                                        >  1:349115/1‑100 (MQ=255)
                                                            cGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACAACAGTGGGCTTTTTTCt                                        >  1:3971944/1‑100 (MQ=255)
                                                             gTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTg                                       <  1:4333282/100‑1 (MQ=255)
                                                             gTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTg                                       <  2:1156578/100‑1 (MQ=255)
                                                                  atatCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCtat                                  <  2:952651/100‑1 (MQ=255)
                                                                   tatCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCtatt                                 >  2:1003882/1‑100 (MQ=255)
                                                                     tCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCtattta                               >  2:1126720/1‑100 (MQ=255)
                                                                     tCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCtattta                               >  2:2598596/1‑100 (MQ=255)
                                                                     tCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCtattta                               >  1:911106/1‑100 (MQ=255)
                                                                       ccGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAATCTCATTTGTCACCACAGTGGGCTTTTTTCTGTCtatttatt                             >  1:3846810/1‑100 (MQ=255)
                                                                       ccGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCtatttatt                             >  2:4243870/1‑100 (MQ=255)
                                                                       ccGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCtatttatt                             >  1:4089216/1‑100 (MQ=255)
                                                                         gATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTAtttt                           <  1:896626/100‑1 (MQ=255)
                                                                            gCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGcc                        <  1:4326017/100‑1 (MQ=255)
                                                                             cAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCg                       >  1:1500389/1‑100 (MQ=255)
                                                                             cAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCg                       >  1:2021308/1‑100 (MQ=255)
                                                                               ggAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGAt                     >  2:3011125/1‑100 (MQ=255)
                                                                               ggAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGAt                     >  1:4229125/1‑100 (MQ=255)
                                                                               ggAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGAt                     >  1:2226534/1‑100 (MQ=255)
                                                                               ggAATATTTCCATTTGACCAATGGTCAAATTGGTAATGATATGTCGGTAAACTCATTTGTCACCACAATGGGCTGTTTTCTGTCTATTTATTTTGCCGAt                     >  1:3924126/1‑100 (MQ=255)
                                                                                 aaTATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATaa                   >  1:2680830/1‑100 (MQ=255)
                                                                                 aaTATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATaa                   <  1:3129642/100‑1 (MQ=255)
                                                                                  cattaccaaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATaaa                  >  2:3846961/9‑100 (MQ=39)
                                                                                    atTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACt                <  1:161214/100‑1 (MQ=255)
                                                                                     tttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTa               >  2:2718125/1‑100 (MQ=255)
                                                                                      ttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTAc              <  2:2285824/100‑1 (MQ=255)
                                                                                      accaaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTAc              >  2:1429948/5‑100 (MQ=255)
                                                                                       tCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTAcc             <  1:2605368/100‑1 (MQ=255)
                                                                                       ccaaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTAcc             >  2:2157798/4‑100 (MQ=255)
                                                                                        ccATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACcg            >  1:3407562/1‑100 (MQ=255)
                                                                                        caaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACcg            >  2:168542/3‑100 (MQ=255)
                                                                                         cATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACcgc           >  2:3659451/1‑100 (MQ=255)
                                                                                         aaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACcgc           >  2:3747297/2‑100 (MQ=255)
                                                                                         aaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACcgc           >  2:2468663/2‑100 (MQ=255)
                                                                                         aaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACcgc           >  2:889095/2‑100 (MQ=255)
                                                                                         aaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACcgc           >  1:773660/2‑100 (MQ=255)
                                                                                            ttGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCa        <  2:4001327/100‑1 (MQ=255)
                                                                                             tGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAg       <  2:773662/100‑1 (MQ=255)
                                                                                               aCCAATGGGCAAATTGGTAATGCAATGACGGTAAACTCATTTGTGACAAAAGAGGGCTTTTTTCTGTCTATAAATTTTGCCGATAAACTACCGCGCAGAt     >  1:3257631/1‑100 (MQ=255)
                                                                                                ccAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGAt     >  2:3324556/1‑99 (MQ=255)
                                                                                                ccAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATa    >  2:243627/1‑100 (MQ=255)
                                                                                                ccAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATa    >  1:1890098/1‑100 (MQ=255)
                                                                                                  aaTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATaca  <  1:3067778/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
TATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACA  >  NC_000913/4062272‑4062469

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