Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,903,427 G→A R257C (CGT→TGT)  bglG ← BglG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,903,4270GA86.3% 814.0 / 89.7 277R257C (CGT→TGT) bglGBglG
Reads supporting (aligned to +/- strand):  ref base G (17/21);  new base A (134/105);  total (151/126)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 278 total aligned reads displayed.

ATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACTGCTTGTTGTAATGATTC  >  NC_000913/3903336‑3903517
                                                                                           |                                                                                          
aTCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGt                                                                                           <  2:3246641/93‑1 (MQ=255)
 tCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACAtt                                                                                          <  2:5683382/93‑1 (MQ=255)
 tCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACAtt                                                                                          >  2:2445952/1‑93 (MQ=255)
 tCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACAtt                                                                                          >  1:5970634/1‑93 (MQ=255)
  ccTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTg                                                                                         >  2:655780/1‑93 (MQ=255)
  ccTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTg                                                                                         >  2:5332438/1‑93 (MQ=255)
   cTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGa                                                                                        <  1:4246642/93‑1 (MQ=255)
    tttACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGAt                                                                                       >  2:4650235/1‑93 (MQ=255)
     ttACTCAGTAATAATATTTCAGTGTTATTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATa                                                                                      >  2:3504642/1‑93 (MQ=255)
      tACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATAc                                                                                     <  1:1083595/93‑1 (MQ=255)
        cTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTg                                                                                   <  2:1165671/93‑1 (MQ=255)
        cTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTg                                                                                   >  2:4270253/1‑93 (MQ=255)
         tCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGc                                                                                  <  1:2980917/93‑1 (MQ=255)
         tCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGc                                                                                  <  2:2865430/93‑1 (MQ=255)
         tCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGc                                                                                  >  1:1502702/1‑93 (MQ=255)
            gTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCaaa                                                                               >  2:280418/1‑93 (MQ=255)
            gTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCaaa                                                                               <  1:4839098/93‑1 (MQ=255)
            gTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCaaa                                                                               >  1:4231506/1‑93 (MQ=255)
             taataaTATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAAc                                                                              >  1:1037729/1‑93 (MQ=255)
              aataatATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAAcc                                                                             <  1:5267389/93‑1 (MQ=255)
              aataatATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAAcc                                                                             >  2:3275615/1‑93 (MQ=255)
                  atatTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATa                                                                         >  1:5474488/1‑93 (MQ=255)
                   tatTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATaa                                                                        >  1:237860/1‑93 (MQ=255)
                    atTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATaaa                                                                       >  2:3346584/1‑93 (MQ=255)
                    atTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATaaa                                                                       >  1:3599374/1‑93 (MQ=255)
                    atTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATaaa                                                                       >  1:3657961/1‑93 (MQ=255)
                    atTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATaaa                                                                       >  1:4730601/1‑93 (MQ=255)
                     tttCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATaaaa                                                                      <  1:1524049/93‑1 (MQ=255)
                     tttCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATaaaa                                                                      <  2:2950211/93‑1 (MQ=255)
                     tttCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATaaaa                                                                      >  1:4087903/1‑93 (MQ=255)
                       tCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAAt                                                                    <  2:2521523/93‑1 (MQ=255)
                        cAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATg                                                                   <  2:2383146/93‑1 (MQ=255)
                          gtgtTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGc                                                                 >  1:1579730/1‑93 (MQ=255)
                          gtgtTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGc                                                                 <  1:3912025/93‑1 (MQ=255)
                           tgtTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCg                                                                >  2:6495820/1‑93 (MQ=255)
                           tgtTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCg                                                                <  1:2582695/93‑1 (MQ=255)
                             ttCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGAt                                                              <  1:4615866/93‑1 (MQ=255)
                              tCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATc                                                             >  2:1035439/1‑93 (MQ=255)
                              tCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATc                                                             <  1:324901/93‑1 (MQ=255)
                              tCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATc                                                             <  1:1879709/93‑1 (MQ=255)
                               cTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATcc                                                            >  1:6418930/1‑93 (MQ=255)
                                tttGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCg                                                           <  1:922336/93‑1 (MQ=255)
                                tttGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCg                                                           <  2:3077104/93‑1 (MQ=255)
                                tttGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCg                                                           <  1:4376981/93‑1 (MQ=255)
                                 ttGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGt                                                          <  2:2667085/93‑1 (MQ=255)
                                 ttGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGt                                                          <  1:3081854/93‑1 (MQ=255)
                                 ttGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGt                                                          <  2:4997748/93‑1 (MQ=255)
                                   gcgcACGCTCTTTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGATCCGTTc                                                        <  2:3700241/93‑1 (MQ=255)
                                      cACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGc                                                     >  2:5386977/1‑93 (MQ=255)
                                        cgcgCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGcac                                                   >  2:4657193/1‑93 (MQ=255)
                                        cgcgCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGcac                                                   <  1:115986/93‑1 (MQ=255)
                                          cgcTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACAt                                                 >  1:5537101/1‑93 (MQ=255)
                                            ctctATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTg                                               >  1:3819944/1‑93 (MQ=255)
                                            ctctATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTg                                               >  1:2263080/1‑93 (MQ=255)
                                            ctctATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTg                                               >  2:6541029/1‑93 (MQ=255)
                                            ctctATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTg                                               >  1:4954674/1‑93 (MQ=255)
                                             tctATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGc                                              >  2:5765432/1‑93 (MQ=255)
                                              ctATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGcc                                             >  2:6311122/1‑93 (MQ=255)
                                               tataTTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCa                                            <  1:6533317/93‑1 (MQ=255)
                                                atatTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCAt                                           <  2:3902340/93‑1 (MQ=255)
                                                  atTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGc                                         <  2:1554666/93‑1 (MQ=255)
                                                  atTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGc                                         <  1:613057/93‑1 (MQ=255)
                                                   tttATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCACGCt                                        >  2:1902056/1‑93 (MQ=255)
                                                      aTGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGc                                     >  2:3725468/1‑93 (MQ=255)
                                                       tGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCg                                    >  1:4911983/1‑93 (MQ=255)
                                                           tAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTa                                >  1:5870119/1‑93 (MQ=255)
                                                           tAAAAACATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTa                                >  1:2510149/1‑93 (MQ=255)
                                                             aaaaCATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAAt                              <  2:1612758/93‑1 (MQ=255)
                                                             aaaaCATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAAt                              >  1:5390030/1‑93 (MQ=255)
                                                               aaCATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAAttt                            >  1:2018089/1‑93 (MQ=255)
                                                               aaCATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAAttt                            >  2:5918138/1‑93 (MQ=255)
                                                                aCATAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAAtttt                           >  2:6607113/1‑93 (MQ=255)
                                                                   tAATCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGtt                        >  2:85617/1‑93 (MQ=255)
                                                                    aaTCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGttt                       <  2:5841087/93‑1 (MQ=255)
                                                                    aaTCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGttt                       <  2:1218861/93‑1 (MQ=255)
                                                                     aTCTCTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGtttt                      >  1:4996511/1‑93 (MQ=255)
                                                                      tctcTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTa                     <  1:1438145/93‑1 (MQ=255)
                                                                      tctcTGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTa                     >  2:784729/1‑93 (MQ=255)
                                                                        tctGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACt                   >  2:400006/1‑93 (MQ=255)
                                                                        tctGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACt                   >  1:5910527/1‑93 (MQ=255)
                                                                          tGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACTGc                 <  2:1507020/93‑1 (MQ=255)
                                                                          tGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTc                                                        >  2:1601812/1‑54 (MQ=255)
                                                                          tGCGGGTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACTGc                 <  2:3828642/93‑1 (MQ=255)
                                                                             gggTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACTGCttg              >  2:443856/1‑93 (MQ=255)
                                                                             gggTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACTGCttg              >  1:2788444/1‑93 (MQ=255)
                                                                              ggTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACTGCttgt             <  2:6137869/93‑1 (MQ=255)
                                                                               gTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACTGCttgtt            <  2:4166128/93‑1 (MQ=255)
                                                                               gTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACTGCttgtt            <  1:1374598/93‑1 (MQ=255)
                                                                               gTGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACTGCttgtt            <  1:2198091/93‑1 (MQ=255)
                                                                                tGAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACTGCttgttg           <  2:2786956/93‑1 (MQ=255)
                                                                                 gAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACTGCTtgttgt          <  1:2363817/93‑1 (MQ=255)
                                                                                 gAAATTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACTGCTtgttgt          >  2:2762358/1‑93 (MQ=255)
                                                                                  aaaTTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACTGCTTGTTGTa         <  2:480570/93‑1 (MQ=255)
                                                                                    aTTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACTGCTTGTTGTAAt       >  2:5514847/1‑93 (MQ=255)
                                                                                    aTTTTACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGGAATTTTGTTTTACTGCTTGTTGTAAt       >  2:2110157/1‑93 (MQ=255)
                                                                                     ttttACGTTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACTGCTTGTTGTGATg      <  1:2913930/93‑1 (MQ=255)
                                                                                     ttttACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACTGCTTGTTGTAATg      >  1:1787749/1‑93 (MQ=255)
                                                                                     ttttACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACTGCTTGTTGTAATg      <  2:4517862/93‑1 (MQ=255)
                                                                                      tttACATTGATACTGCAAACCAATAAAAATGGCGATACGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACTGCTTGTTGTAATGa     <  1:973915/93‑1 (MQ=255)
                                                                                        tACATTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACTGCTTGTTGTAATGAtt   <  1:742999/93‑1 (MQ=255)
                                                                                           |                                                                                          
ATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGATCCGTTCCGCACATTGCCATGCTTGCGGGTAATTTTGTTTTACTGCTTGTTGTAATGATTC  >  NC_000913/3903336‑3903517

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