Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,673,475 A→G 50.4% G9G (GGT→GGC mdtJ ← multidrug efflux system transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,673,4750AG50.4% 1.3 / 147.0 111G9G (GGT→GGCmdtJmultidrug efflux system transporter
Reads supporting (aligned to +/- strand):  ref base A (30/25);  new base G (31/25);  total (61/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

Only 100 of 112 total aligned reads displayed.

AAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCC  >  NC_000913/1673378‑1673573
                                                                                                 |                                                                                                  
aaaTCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGcc                                                                                                  >  1:1504486/1‑100 (MQ=255)
  atcatcACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCta                                                                                                >  2:2768275/1‑100 (MQ=255)
   tcatcaCCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGGCGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCta                                                                                                <  2:2899377/99‑1 (MQ=255)
    catcaCCCGCATTAAAATAAAGCCGCCATTCCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTAGCGGTAATTTCTGTAGCAATAGCCAGACCtaat                                                                                              <  1:3925946/100‑1 (MQ=255)
    catcaCCAGCATTGAAATAAAGCCGCCATTTCCCTCACAGACGCTCGCCCATTTCATTGACAGCGTACCGGTACTTTCTGTAGCAATAGCCAGGCCtaat                                                                                              >  2:2559169/1‑100 (MQ=255)
    catcaCCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCtaat                                                                                              >  1:345274/1‑100 (MQ=255)
     atcaCCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCtaata                                                                                             >  1:58773/1‑100 (MQ=255)
     atcaCCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCtaata                                                                                             >  2:1210103/1‑100 (MQ=255)
      tcaCCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCtaataa                                                                                            <  2:228427/100‑1 (MQ=255)
       caCCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATaaa                                                                                           >  1:2416862/1‑100 (MQ=255)
         ccAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAAt                                                                                         <  1:3513556/100‑1 (MQ=255)
         ccAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAAt                                                                                         <  2:3256024/100‑1 (MQ=255)
          cAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATc                                                                                        <  2:1386607/100‑1 (MQ=255)
          cAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACACTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATc                                                                                        >  1:2906419/1‑100 (MQ=255)
           aGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATcc                                                                                       <  1:901929/100‑1 (MQ=255)
             cATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCaa                                                                                     >  2:1148635/1‑100 (MQ=255)
              aTTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCaat                                                                                    <  2:2755465/100‑1 (MQ=255)
              aTTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCaat                                                                                    <  2:3682604/100‑1 (MQ=255)
              aTTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCaat                                                                                    <  1:2895326/100‑1 (MQ=255)
               ttAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCaata                                                                                   <  2:130963/100‑1 (MQ=255)
               ttAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCaata                                                                                   <  1:2383287/100‑1 (MQ=255)
                tAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCaataa                                                                                  >  1:704802/1‑100 (MQ=255)
                 aaaaTAAAGCCGCCATTTCCCTCACTTACGCTCGCCCCTTTCATTGACAGCGTACCGGTAATTTCTGTACCAATAAACAGGCCTAATAACATCCaattaa                                                                                 >  2:3606547/1‑100 (MQ=255)
                 aaaaTAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCaataaa                                                                                 >  1:3695295/1‑100 (MQ=255)
                 aaaaTAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCaataaa                                                                                 >  2:3903367/1‑100 (MQ=255)
                 aaaaTAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCaataaa                                                                                 <  1:1210111/100‑1 (MQ=255)
                  aaaTAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAACTCCaataaat                                                                                >  2:2617615/1‑100 (MQ=255)
                  aaaTAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCaataaat                                                                                >  2:612123/1‑100 (MQ=255)
                   aaTAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCaataaata                                                                               <  2:1442271/100‑1 (MQ=255)
                         gccgccATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACAtt                                                                         >  1:1092517/1‑100 (MQ=255)
                         gccgccATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACAtt                                                                         >  1:2374925/1‑100 (MQ=255)
                         gccgccATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACAtt                                                                         >  1:3817828/1‑100 (MQ=255)
                          ccgccATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTTTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTg                                                                        <  2:109998/100‑1 (MQ=255)
                          ccgccATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTg                                                                        >  1:1577349/1‑100 (MQ=255)
                          cagccATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATATAATCCAATAAATATACATTg                                                                        <  2:2515395/98‑1 (MQ=255)
                           cgccATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGt                                                                       <  1:183639/100‑1 (MQ=255)
                             ccATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTcc                                                                     <  2:40740/100‑1 (MQ=255)
                              cATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCt                                                                    <  1:3193822/100‑1 (MQ=255)
                              cATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCt                                                                    <  2:378062/100‑1 (MQ=255)
                                  tCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTtctc                                                                <  2:1004322/100‑1 (MQ=255)
                                  tCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTtctc                                                                >  1:2144367/1‑100 (MQ=255)
                                   cccTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTcc                                                               >  2:258121/1‑100 (MQ=255)
                                    ccTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCt                                                              <  2:2494299/100‑1 (MQ=255)
                                         cTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAag                                                         <  2:578106/100‑1 (MQ=255)
                                           gACGCTCGCCCATTTCATTGACAGCGTACCAGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAagag                                                       <  2:548682/100‑1 (MQ=255)
                                             cGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGaa                                                     >  1:477476/1‑100 (MQ=255)
                                             cGCTCCCCCATTTCATTGACAAAGGACCGGAACTGTCTATAGCAAAAACCACGCATAATAACATCCAATAAATATACATTGTCCTTccccgcctaga                                                        >  2:602136/1‑90 (MQ=255)
                                                tCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAatta                                                  <  2:2226242/100‑1 (MQ=255)
                                                 cGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCTTGCAGGAGaaggac                                                 <  1:3026675/100‑5 (MQ=255)
                                                 cGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAattat                                                 <  2:345271/100‑1 (MQ=255)
                                                  gCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAattatt                                                >  1:563637/1‑100 (MQ=255)
                                                  gCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAattatt                                                >  1:3659431/1‑100 (MQ=255)
                                                  gCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGACCTTCTCCTGCAAGAGAattatt                                                >  2:1620574/1‑100 (MQ=255)
                                                  gCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAattatt                                                >  1:2309241/1‑100 (MQ=255)
                                                  gCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAattat                                                 >  2:2277514/1‑99 (MQ=255)
                                                  gCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAACAGCCAGGCCTAATAAAATCCAATTAATATACATTGTCCTTCTCCTGCAAGATAattctt                                                <  2:1763696/100‑1 (MQ=255)
                                                   cccATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTAttt                                               <  2:3994948/100‑1 (MQ=255)
                                                     cATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTa                                             >  1:2116786/1‑100 (MQ=255)
                                                      aTTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTaa                                            <  2:498651/100‑1 (MQ=255)
                                                      aTTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTaa                                            <  2:3238218/100‑1 (MQ=255)
                                                       tttCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAAt                                           >  2:304645/1‑100 (MQ=255)
                                                       tttCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAAt                                           >  2:4109600/1‑100 (MQ=255)
                                                          cATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAAtttt                                        >  1:73343/1‑100 (MQ=255)
                                                            ttGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTGTCCTGCAAGAGAATTATTTTAATTTTCg                                      <  1:1343057/100‑1 (MQ=255)
                                                            ttGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCg                                      <  1:1551279/100‑1 (MQ=255)
                                                            ttGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCg                                      >  1:2119058/1‑100 (MQ=255)
                                                             tGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGc                                     <  2:794827/100‑1 (MQ=255)
                                                             tGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGc                                     <  2:182274/100‑1 (MQ=255)
                                                              gACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCATTCTCCTGCAAGAGAATTATTTTAAtttt                                        >  2:2915148/1‑96 (MQ=255)
                                                                 aGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTaa                                 <  1:1286990/100‑1 (MQ=255)
                                                                   cGTACCGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAAtt                               >  1:3271447/1‑100 (MQ=255)
                                                                   cGTACCGGTAATTTCTGTAGCAATAGCCAGACCAAAGAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTCTTTTAATTTTCGCTTAAtc                               >  2:3454102/1‑99 (MQ=255)
                                                                      aCCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCAGCAAGAGAATTATATTAATTTTAGCTTAATTCAg                            >  2:3653375/1‑100 (MQ=255)
                                                                       ccGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCa                             >  2:2047879/1‑98 (MQ=255)
                                                                        cGGTAATTTCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCg                          <  1:2092413/100‑1 (MQ=255)
                                                                        cGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCg                          <  1:1144507/100‑1 (MQ=255)
                                                                         ggTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGa                         >  1:379158/1‑100 (MQ=255)
                                                                         ggTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGa                         >  1:254246/1‑100 (MQ=255)
                                                                         ggTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGa                         <  1:3651593/100‑1 (MQ=255)
                                                                         gatAATTTCTGTAGCAATAGCCAGAGCCAAGAAAATCCAAGAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGc                           <  2:3289512/96‑1 (MQ=255)
                                                                             aTTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAAc                      >  2:2832655/1‑99 (MQ=255)
                                                                              tttCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCa                    <  1:3063876/100‑1 (MQ=255)
                                                                              tttCTGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTGAATTAAGCGaaaatt                    <  2:3514262/100‑4 (MQ=255)
                                                                              tttCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCa                    <  1:3380882/100‑1 (MQ=255)
                                                                               ttCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaa                   >  1:1080338/1‑100 (MQ=255)
                                                                                 cTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaag                 >  1:2058255/1‑100 (MQ=255)
                                                                                  tGTAGCAGTAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaagt                >  1:2164228/1‑100 (MQ=255)
                                                                                  tGTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaagt                <  1:2047892/100‑1 (MQ=255)
                                                                                   gTAGCAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaagta               >  1:85944/1‑100 (MQ=255)
                                                                                   gTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaagta               >  1:817981/1‑100 (MQ=255)
                                                                                      gCAATATCCAGACCTAAGAAAATCCAATAAATAGACATTGTCCTTCTCCTGCAAGAGAATTATTTTAAGTTTCTCTTAATTCAGCGAAACCaaattaaa             >  2:3448313/1‑99 (MQ=255)
                                                                                      gCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaagtaaag            >  1:3896192/1‑100 (MQ=255)
                                                                                      gCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaagtaaag            <  2:2091798/100‑1 (MQ=255)
                                                                                       cAATAGCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaagtaaagt           >  1:1410643/1‑100 (MQ=255)
                                                                                       cAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaagtaaagt           <  1:3209660/100‑1 (MQ=255)
                                                                                       cAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCATTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaagtaaagt           >  2:2261921/1‑100 (MQ=255)
                                                                                            gCCAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCa      >  1:1841807/1‑100 (MQ=255)
                                                                                              cAGGCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAg    >  1:236401/1‑100 (MQ=255)
                                                                                              cAGACCTCATAAAATCCAATACATTTACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAg    <  2:1778420/100‑1 (MQ=255)
                                                                                                 aCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTTGTCagccg  >  2:804260/1‑95 (MQ=255)
                                                                                                 |                                                                                                  
AAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCC  >  NC_000913/1673378‑1673573

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: