Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 469,244 A→T 32.0% D125V (GAT→GTT)  mdlA → putative multidrug ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913469,2440AT32.0% 112.8 / 73.0 87D125V (GAT→GTT) mdlAputative multidrug ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (33/26);  new base T (18/10);  total (51/36)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

CGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTG  >  NC_000913/469147‑469343
                                                                                                 |                                                                                                   
cGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATc                                                                                                   >  1:387437/1‑100 (MQ=255)
 gTGAAGATTATTACCGTCAGCTAAGCCGTCAGCATCGTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATcg                                                                                                  >  2:2468801/1‑100 (MQ=255)
 gTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGATTTTTACCAGCGACATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATcg                                                                                                  <  2:378629/100‑1 (MQ=255)
 gTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATcg                                                                                                  >  2:1799227/1‑100 (MQ=255)
     aGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGcgtc                                                                                              >  1:1045976/1‑100 (MQ=255)
      gATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGcgtcg                                                                                             <  2:1439018/100‑1 (MQ=255)
      gATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGcgtcg                                                                                             >  1:2264927/1‑100 (MQ=255)
       attattACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGTTCGcgtcgt                                                                                            <  2:2495315/100‑1 (MQ=255)
          attaCCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGtt                                                                                         >  2:1730070/1‑100 (MQ=255)
          attaCCGTCAGCTAAGCCGGAAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGtt                                                                                         >  2:827054/1‑100 (MQ=255)
            taCCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTg                                                                                       <  2:905455/100‑1 (MQ=255)
             aCCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATTGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTgc                                                                                      >  1:1859250/1‑100 (MQ=255)
              ccGTCAGCTAAGCGGGCAGCATCCTGAGTTTTACCTGCTTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGGGTTTgc                                                                                      >  2:1605434/1‑99 (MQ=255)
               cGTCCGCTAAGCCGGCAGCATCCTGAGATTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGCCGATCGCGTCGTGTTTgc                                                                                      >  1:563584/1‑98 (MQ=255)
               cGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGTTCGCGTCGTGTTTgccg                                                                                    >  2:1946872/1‑100 (MQ=255)
               cGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTgccg                                                                                    >  2:959497/1‑100 (MQ=255)
                gTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGTTCGCGTCGTGTTTgccgc                                                                                   >  1:996843/1‑100 (MQ=255)
                  cAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGTTCGCGTCGTGTTTGccgccg                                                                                 >  2:1853183/1‑100 (MQ=255)
                  cAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGccgccg                                                                                 >  1:125112/1‑100 (MQ=255)
                   aGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCgg                                                                                <  2:2318750/100‑1 (MQ=255)
                    gCTAAGCCGGCAGCATGCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGga                                                                               <  2:382284/100‑1 (MQ=255)
                    gCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGTTCGCGTCGTGTTTGCCGCCgg                                                                                >  2:2491712/1‑99 (MQ=255)
                    gCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGTTCGCGTCGTGTTTGCCGCCGga                                                                               >  2:2491653/1‑100 (MQ=255)
                    gCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGga                                                                               <  1:1122001/100‑1 (MQ=255)
                    gCTAAGCCGGCAGCATCCTGAGTTTTACCTGCATCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCgg                                                                                <  2:1013638/99‑1 (MQ=255)
                      tAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGgaga                                                                             <  1:2302373/100‑1 (MQ=255)
                      tAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGgaga                                                                             >  2:159962/1‑100 (MQ=255)
                        aGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAg                                                                           >  2:1477127/1‑100 (MQ=255)
                         gCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAgg                                                                          >  1:195512/1‑100 (MQ=255)
                           cGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAgggg                                                                        <  1:1976730/100‑1 (MQ=255)
                             gcagcaTCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGTTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTg                                                                      <  2:673733/100‑1 (MQ=255)
                             gcagcaTCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTg                                                                      <  2:2264914/100‑1 (MQ=255)
                              cagcaTCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGc                                                                     >  1:2415851/1‑100 (MQ=255)
                               agcaTCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGTTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCt                                                                    <  1:2042175/100‑1 (MQ=255)
                               agcaTCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACCTCGTTCGCGTCGTGTTTGCCGCCGGGGAAGGGGGGc                                                                     >  2:997678/1‑99 (MQ=255)
                                gcaTCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGTTCGCGTCGGGTTTGCCGCCGGAGCAGGGGTGCTg                                                                   >  2:2243745/1‑100 (MQ=255)
                                   tCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGTTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACg                                                                >  1:273692/1‑100 (MQ=255)
                                    ccTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGTTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGc                                                               <  2:1956757/100‑1 (MQ=255)
                                     cTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGTTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCt                                                              >  2:1228404/1‑100 (MQ=255)
                                     cTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGTTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCt                                                              >  2:537217/1‑100 (MQ=255)
                                     cTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCt                                                              >  2:117703/1‑100 (MQ=255)
                                     cTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCt                                                              >  2:2113184/1‑100 (MQ=255)
                                         gTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGTTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCtggtg                                                          <  2:1306921/100‑1 (MQ=255)
                                         gTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCtggtg                                                          >  1:2228040/1‑100 (MQ=255)
                                          ttttACCTGCGTCATCGCACCGGTGACCTCACGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCtggtg                                                          <  2:2372102/99‑1 (MQ=255)
                                            ttACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGAt                                                       <  2:2122835/100‑1 (MQ=255)
                                               ccTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGTTCGCGTCGTGTTTGCCGCCGGCGAAGGGGTGCTGACGCTGGTGGATTCa                                                    >  2:1982697/1‑100 (MQ=255)
                                                cTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCAc                                                   <  1:2354488/100‑1 (MQ=255)
                                                cTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCAc                                                   <  1:1477135/100‑1 (MQ=255)
                                                cTGCGTCATCGCACCGGTGACCTCATGGATCGTGCGACCAATGAAGTCGATCGCGTCGTGTTTTCCGCCGGAGAAGGGGTGCGGAGGCTGGTGGATTCAc                                                   >  2:592204/1‑100 (MQ=255)
                                                 tGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGTTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACt                                                  >  1:566209/1‑100 (MQ=255)
                                                  gCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGTTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTg                                                 >  1:1751282/1‑100 (MQ=255)
                                                  gCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTg                                                 <  2:195510/100‑1 (MQ=255)
                                                   cGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTgg                                                >  1:2159769/1‑100 (MQ=255)
                                                     tCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTg                                              <  1:794089/100‑1 (MQ=255)
                                                            aCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGTTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCt                                       >  2:962716/1‑100 (MQ=255)
                                                                gTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTgcgc                                   <  2:377360/100‑1 (MQ=255)
                                                                 tcaCCTCATTGCTCCTGCGACTAATTATCTCGATCGCGTCGTGTTTGCCGCCGGAGAATGGGTTCTGACGTTTGTGTTTTCACTGGTGATGGGCTgcg                                    <  2:16300/96‑1 (MQ=255)
                                                                     cTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTgtgt                              <  2:1980513/100‑1 (MQ=255)
                                                                      tCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCATTGGTGATGGGCTGCGCTGTGtt                             <  2:1511161/100‑1 (MQ=255)
                                                                      tCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGtt                             <  1:1185586/100‑1 (MQ=255)
                                                                      tCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGtt                             >  1:1534033/1‑100 (MQ=255)
                                                                       cATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTg                            >  1:2153665/1‑100 (MQ=255)
                                                                       cATGGCTCGTGCGACCAACGGCGTCGTTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTg                            <  2:2236186/100‑1 (MQ=255)
                                                                        aTGGCTCGTGCGACCAATGCCGTCGTTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGGGGATTCACTGGTGATGGGCTGCGCGGTGTTg                            >  2:684749/1‑99 (MQ=255)
                                                                         tGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCTCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGAt                          >  2:1421667/1‑100 (MQ=255)
                                                                         tGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGAt                          >  1:839957/1‑100 (MQ=255)
                                                                          ggCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGAGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGAtt                         <  2:584161/100‑1 (MQ=255)
                                                                              cGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTatga                     <  2:887144/100‑1 (MQ=255)
                                                                               gTGCGACCAATGACGTCGTTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTatgat                    >  1:328964/1‑100 (MQ=255)
                                                                               gTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTatgat                    <  2:810392/100‑1 (MQ=255)
                                                                                tGCGACCAATGACGTCGTTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTatgatg                   >  2:555724/1‑100 (MQ=255)
                                                                                tGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTatgatg                   >  1:153482/1‑100 (MQ=255)
                                                                                tGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTatgatg                   <  1:2113195/100‑1 (MQ=255)
                                                                                tGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGATGTGTTGATTatgatg                   >  2:1519523/1‑100 (MQ=255)
                                                                                   gACCAATGACGTCGTTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCt                <  2:1751273/100‑1 (MQ=255)
                                                                                   gACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCt                >  1:101398/1‑100 (MQ=255)
                                                                                    aCCAATGACGTCGTTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTa               <  2:2415703/100‑1 (MQ=255)
                                                                                     ccAATGACGTCGTTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTAc              <  2:2414871/100‑1 (MQ=255)
                                                                                      cAATGACGTCGTTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACg             <  2:860619/100‑1 (MQ=255)
                                                                                       aaTGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGc            <  1:777364/100‑1 (MQ=255)
                                                                                         tGAGGTCGTTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCaa          >  1:1802909/1‑100 (MQ=255)
                                                                                         tGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCaa          >  1:247030/1‑100 (MQ=255)
                                                                                            cGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAAtt       >  2:111339/1‑100 (MQ=255)
                                                                                            cGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAAtt       >  1:1591335/1‑100 (MQ=255)
                                                                                            cGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAAtt       >  1:97500/1‑100 (MQ=255)
                                                                                                gATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCt   >  2:1382200/1‑100 (MQ=255)
                                                                                                 aTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTg  <  2:1356815/100‑1 (MQ=255)
                                                                                                 aTCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTg  <  2:1045971/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
CGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTG  >  NC_000913/469147‑469343

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: ‑