Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,062,371 A→T 18.4% 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,3710AT18.4% 295.2 / 59.2 130N42I (AAT→ATT) yihNMFS transporter family protein
Reads supporting (aligned to +/- strand):  ref base A (48/57);  new base T (12/12);  total (60/70)
Fisher's exact test for biased strand distribution p-value = 8.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

Only 100 of 130 total aligned reads displayed.

ATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATAC  >  NC_000913/4062273‑4062468
                                                                                                  |                                                                                                 
aTTTAGTCTATTAACACTGTGTGGCGTTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGAAGGAATATTTCCATTTGACCAAt                                                                                                  >  1:1297304/1‑100 (MQ=255)
 tttAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATg                                                                                                 >  1:1861170/1‑100 (MQ=255)
   tAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGt                                                                                               >  2:1702111/1‑100 (MQ=255)
    aGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTc                                                                                              <  2:1886136/100‑1 (MQ=255)
    aGTCTATTAACACTGTGTAGCGGTACAAATTCTAAATTACCGTCGCTGAAAAATGCGTCTTATATCCCGATGCAGAAATATTTCCATTTGACCACTGGt                                                                                               <  1:2225793/99‑1 (MQ=255)
     gTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCa                                                                                             >  2:789294/1‑100 (MQ=255)
       cTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCaaa                                                                                           >  2:1114926/1‑100 (MQ=255)
       cTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCaaa                                                                                           >  2:1832263/1‑100 (MQ=255)
        tATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAAt                                                                                          >  2:408304/1‑100 (MQ=255)
        tATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAAt                                                                                          >  1:1766278/1‑100 (MQ=255)
         attaggacTGTGTGGCGGTACAATTTGTAGAGTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTGCCATTTGACCAATGGTCAAAt                                                                                          <  1:2675149/93‑1 (MQ=255)
         aTTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAtg                                                                                         <  1:2646103/100‑2 (MQ=255)
         aTTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAtg                                                                                         <  1:1274265/100‑2 (MQ=255)
         aTTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAtg                                                                                         <  1:1470370/100‑2 (MQ=255)
         aTTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAtg                                                                                         <  2:229380/100‑2 (MQ=255)
         aTTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAtg                                                                                         <  1:1621486/100‑2 (MQ=255)
         aTTAACACGGTGTGGCGGTACGATTTAAAAATTGCCGTAGCTGAAAGATGCGTTCTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAtc                                                                                         <  1:852782/100‑2 (MQ=255)
          ttAACACTGTGTGGCTGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAtgg                                                                                        <  1:2430665/100‑3 (MQ=255)
          ttAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAtgg                                                                                        <  2:482587/100‑3 (MQ=255)
          ttAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCATTGGTCAAAtgg                                                                                        <  2:1939872/100‑3 (MQ=255)
            aaCACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGt                                                                                      >  2:2323755/1‑100 (MQ=255)
            aaCACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGt                                                                                      <  1:1844631/100‑1 (MQ=255)
             acacTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTa                                                                                     >  2:1365283/1‑100 (MQ=255)
              cacTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTGCCATTTGACCAATGGTCGAATTGGTaa                                                                                    >  1:2869013/1‑100 (MQ=255)
               acTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAAt                                                                                   <  1:1368434/100‑1 (MQ=255)
                 tgtgtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTTAAATTGGTAAtac                                                                                 <  2:1974479/100‑3 (MQ=255)
                 tgtgtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGc                                                                                 <  1:20025/100‑1 (MQ=255)
                  gtgtgGCGGTACAATTTATAAATTACCTTCGCTGAACGATGCGATTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCt                                                                                <  1:1031274/100‑1 (MQ=255)
                  gtgtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTGTATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCt                                                                                <  1:2131443/100‑1 (MQ=255)
                  gtgtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCt                                                                                >  2:2414906/1‑100 (MQ=255)
                   tgtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTa                                                                               <  1:1364473/100‑1 (MQ=255)
                   tgtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTa                                                                               <  1:2850578/100‑1 (MQ=255)
                    gtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTAt                                                                              >  1:1291711/1‑100 (MQ=255)
                    gtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTAt                                                                              <  2:1313991/100‑1 (MQ=255)
                    gtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTAt                                                                              >  2:12491/1‑100 (MQ=255)
                       gCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTc                                                                           <  2:421155/100‑1 (MQ=255)
                       gCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTCGACCAATGGTCAAATTGGTAATGCTATGTc                                                                           <  2:2643249/100‑1 (MQ=255)
                        cGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCg                                                                          <  2:1291713/100‑1 (MQ=255)
                           tACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTa                                                                       <  1:1575524/100‑1 (MQ=255)
                            aCAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTaa                                                                      >  2:2589286/1‑100 (MQ=255)
                              aaTTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAAc                                                                    >  2:1819236/1‑100 (MQ=255)
                                 ttATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCa                                                                 <  2:2057201/100‑1 (MQ=255)
                                   ataAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCAtt                                                               <  2:932953/100‑1 (MQ=255)
                                   ataAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCAtt                                                               <  1:1187624/100‑1 (MQ=255)
                                    taAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCAttt                                                              >  2:2118991/1‑100 (MQ=255)
                                      aaTTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGt                                                            <  1:10555/100‑1 (MQ=255)
                                          aCCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTcacc                                                        <  2:1216731/100‑1 (MQ=255)
                                             gTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACcaca                                                     >  2:1070751/1‑100 (MQ=255)
                                             gTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACcaca                                                     >  1:1076626/1‑100 (MQ=255)
                                                 cTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTgg                                                 >  1:800542/1‑100 (MQ=255)
                                                  tGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTggg                                                >  1:663734/1‑100 (MQ=255)
                                                   gAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGc                                               >  1:1977191/1‑100 (MQ=255)
                                                    aaaGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCt                                              <  2:2802808/100‑1 (MQ=255)
                                                     aaGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCtt                                             <  2:1431510/100‑1 (MQ=255)
                                                       gATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCtttt                                           >  1:926090/1‑100 (MQ=255)
                                                        aTGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCttttt                                          <  2:1931287/100‑1 (MQ=255)
                                                        aTGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCttttt                                          <  2:1544520/100‑1 (MQ=255)
                                                          gCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACAAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTc                                        <  1:105966/100‑1 (MQ=255)
                                                           cGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCt                                       <  1:408302/100‑1 (MQ=255)
                                                            gTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCGCAGTGGGCTTTTTTCTg                                      >  2:578963/1‑100 (MQ=255)
                                                            gTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTg                                      >  1:2256186/1‑100 (MQ=255)
                                                            gTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCAATCTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTTTg                                      >  1:2038623/1‑100 (MQ=255)
                                                              tttATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTc                                    <  1:2589275/100‑1 (MQ=255)
                                                                tataTCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCGTTTTTCTGTTta                                  >  2:2877562/1‑100 (MQ=255)
                                                                    tCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCtattta                              >  1:428741/1‑100 (MQ=255)
                                                                       cGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTatttattt                           >  2:2184082/1‑100 (MQ=255)
                                                                           gCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGcc                       <  2:90609/100‑1 (MQ=255)
                                                                           gCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGcc                       <  2:1407781/100‑1 (MQ=255)
                                                                            tctgAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCg                      <  1:179284/97‑1 (MQ=255)
                                                                            cAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCg                      >  2:725512/1‑100 (MQ=255)
                                                                             aGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGcc                       <  1:1593461/98‑1 (MQ=255)
                                                                               gAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATa                   <  2:2610052/100‑1 (MQ=255)
                                                                                aaTTTTTCCATTTGACCATTGGTCAAATCGGTAGTGTTTTGTTGGTAAACTCTTTTGTCACCACAGTGGGCTTTTTGTTGTATATTTATTTTGCCGATaa                  <  1:1198033/100‑1 (MQ=255)
                                                                                aaTATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATaa                  >  2:2134938/1‑100 (MQ=255)
                                                                                aaTATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATaa                  >  1:1718241/1‑100 (MQ=255)
                                                                                 cattaccaaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATaaa                 >  2:2769626/9‑100 (MQ=39)
                                                                                   atTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACt               >  1:505059/1‑100 (MQ=255)
                                                                                   atTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACt               >  1:2294667/1‑100 (MQ=255)
                                                                                     accaaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGTCTTTTTTCTGTCTATTTATTTTGCCGATAAACTAc             >  2:2475581/5‑100 (MQ=255)
                                                                                     accaaTTTGACCATTGGTCAAATGGAAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTAcc            >  1:2908716/5‑100 (MQ=38)
                                                                                      tCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTAcc            >  2:1084227/1‑100 (MQ=255)
                                                                                      ccaaTTTGACCATTGGTCAAATTGGTCATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTAcc            >  2:1766793/4‑100 (MQ=255)
                                                                                      ccaaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTACACTCATTTGTCACCACAGTGGGCTTTTTTGTGTCTATTTATTTTGCCGATAAACTAcc            >  2:1471944/4‑100 (MQ=255)
                                                                                      ccaaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTAcc            >  2:2586349/4‑100 (MQ=255)
                                                                                      ccaaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTAcc            >  2:2762059/4‑100 (MQ=255)
                                                                                      ccaaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTAcc            >  1:832665/4‑100 (MQ=255)
                                                                                       caaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACcg           >  2:2554445/3‑100 (MQ=255)
                                                                                       caaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTATCTATTTATTTTGCCGATAAACTACcg           >  2:2336925/3‑100 (MQ=255)
                                                                                        cATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACcgc          <  2:1218000/100‑1 (MQ=255)
                                                                                        aaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATGTTGCCTATAAACTACcg           >  2:246985/2‑99 (MQ=255)
                                                                                           ttGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTGCTGTCTATTTATTTTGCCGATAAACTACCGCGCa       <  1:1432429/100‑1 (MQ=255)
                                                                                            tGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAg      <  1:1377461/100‑1 (MQ=255)
                                                                                            tGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAg      <  1:989972/100‑1 (MQ=255)
                                                                                            tGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTCGCCGATAAACTATCGCGCa       >  1:478269/1‑99 (MQ=255)
                                                                                             gACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCgcgc        >  1:448704/1‑97 (MQ=255)
                                                                                             gACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGa     <  2:3934/100‑1 (MQ=255)
                                                                                               ccAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATa   >  2:2568123/1‑100 (MQ=255)
                                                                                               ccAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATa   >  2:892663/1‑100 (MQ=255)
                                                                                               ccAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATa   <  1:2144155/100‑1 (MQ=255)
                                                                                                cAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATac  <  2:1548581/100‑1 (MQ=255)
                                                                                                  |                                                                                                 
ATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATAC  >  NC_000913/4062273‑4062468

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: