Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 993,255 A→C 15.6% intergenic (+21/+22) pepN → / ← ssuB aminopeptidase N/aliphatic sulfonate ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913993,2550AC15.6% 292.5 / 46.8 115intergenic (+21/+22)pepN/ssuBaminopeptidase N/aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (56/41);  new base C (5/13);  total (61/54)
Fisher's exact test for biased strand distribution p-value = 2.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 116 total aligned reads displayed.

CGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCG  >  NC_000913/993158‑993353
                                                                                                  |                                                                                                  
cggcgCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATa                                                                                                  <  1:1621200/100‑1 (MQ=255)
   cgcTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAgcg                                                                                               <  2:2342228/100‑1 (MQ=255)
    ncTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAgcgc                                                                                              >  1:2572534/2‑100 (MQ=255)
     cTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGcgcc                                                                                             >  2:751617/1‑100 (MQ=255)
      tGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGcgccg                                                                                            >  1:1291240/1‑100 (MQ=255)
        gAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGcgccg                                                                                            >  1:676497/1‑98 (MQ=255)
        gAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCgccgcc                                                                                          <  1:1538495/100‑1 (MQ=255)
         aaCAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCa                                                                                         >  1:166060/1‑100 (MQ=255)
           cAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATTCCCCGAATGGCGGCGCTATTGCCGCCAt                                                                                        <  1:1793757/100‑1 (MQ=25)
           cAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCGCTATTGCCGCCAtt                                                                                       <  1:191990/100‑1 (MQ=38)
           cAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCAtt                                                                                       >  1:776937/1‑100 (MQ=255)
           cAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCAtt                                                                                       >  1:390606/1‑100 (MQ=255)
            aGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCGCTATTGCCGCCATTc                                                                                      <  1:1784254/100‑1 (MQ=38)
            aGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCGCTATTGCCGCCATTc                                                                                      <  1:1856964/100‑1 (MQ=38)
            aGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCGCTATTGCCGCCATTc                                                                                      <  2:2046208/100‑1 (MQ=38)
            aGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTc                                                                                      >  1:2357578/1‑100 (MQ=255)
             gTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCGCTATTGCCGCCATTCg                                                                                     <  2:782470/100‑1 (MQ=38)
             gTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCg                                                                                     >  1:1988632/1‑100 (MQ=255)
             gTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCg                                                                                     <  2:591502/100‑1 (MQ=255)
             gTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCg                                                                                     <  2:590996/100‑1 (MQ=255)
             gTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCg                                                                                     <  2:1055324/100‑1 (MQ=255)
             gTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCg                                                                                     >  1:710461/1‑100 (MQ=255)
               tGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCGCTATTGCCGCCATTCggt                                                                                   <  2:783039/100‑2 (MQ=38)
               tGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCggg                                                                                   >  2:2042345/1‑100 (MQ=255)
               tGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCggg                                                                                   >  2:1622581/1‑100 (MQ=255)
               tGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCggg                                                                                   >  2:1945187/1‑100 (MQ=255)
               tGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCggg                                                                                   >  2:1622414/1‑100 (MQ=255)
                gAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCGCTATTGCCGCCATTCggtt                                                                                  <  2:2105918/100‑3 (MQ=37)
                gAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCgggg                                                                                  >  1:412654/1‑100 (MQ=255)
                gAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCgggg                                                                                  <  2:2067518/100‑1 (MQ=255)
                      gCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTa                                                                            >  1:1680372/1‑100 (MQ=255)
                       cTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAc                                                                           <  1:286240/100‑1 (MQ=255)
                        tGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACGGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAcc                                                                          <  1:869220/100‑1 (MQ=255)
                         ggAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAccc                                                                         >  1:2393960/1‑100 (MQ=255)
                         ggAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAccc                                                                         >  2:2131394/1‑100 (MQ=255)
                         ggAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAccc                                                                         <  2:673006/100‑1 (MQ=255)
                          gAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAcccc                                                                        >  1:2231131/1‑100 (MQ=255)
                             aaTCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGt                                                                     <  1:1129562/100‑1 (MQ=255)
                                ctctctGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGtttt                                                                  >  2:653489/1‑100 (MQ=255)
                                ctctctGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCCCCGCCATTCGGGGAATTTAcccc                                                                        >  1:2178735/1‑94 (MQ=255)
                                    ctGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCa                                                              <  1:736319/100‑1 (MQ=255)
                                     tGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAg                                                             <  2:1486157/100‑1 (MQ=255)
                                     tGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGCCAATAGCGCCGCCATTCGGGGAATTTACTCCTGTTTtctccg                                                             <  2:1755594/100‑3 (MQ=255)
                                      ggCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAgg                                                            <  2:1664773/100‑1 (MQ=255)
                                      ggCGATCTGTACGAGAAGATAACTAAAGCACTGGCTGGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAgg                                                            >  2:1825976/1‑100 (MQ=255)
                                       gCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGc                                                           >  2:488100/1‑100 (MQ=255)
                                            cTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGtt                                                      >  2:1715505/1‑100 (MQ=255)
                                             tGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGttt                                                     >  1:2164277/1‑100 (MQ=255)
                                              gTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTc                                                    >  2:2348645/1‑100 (MQ=255)
                                              gTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTc                                                    >  2:288031/1‑100 (MQ=255)
                                                aCGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAg                                                  <  2:166061/100‑1 (MQ=255)
                                                aCGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAg                                                  >  2:2109746/1‑100 (MQ=255)
                                                 cGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGa                                                 <  2:1988643/100‑1 (MQ=255)
                                                 cGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGa                                                 >  2:498982/1‑100 (MQ=255)
                                                 cGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGa                                                 <  2:739571/100‑1 (MQ=255)
                                                  gagaAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGAt                                                >  2:1955158/1‑100 (MQ=255)
                                                   agaagaTAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGAtt                                               <  2:1704807/100‑1 (MQ=255)
                                                      agaTAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcac                                            >  2:2318579/1‑100 (MQ=255)
                                                      agaTAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcac                                            >  2:593578/1‑100 (MQ=255)
                                                         tAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcaccac                                         >  2:1141935/1‑100 (MQ=255)
                                                           aCTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGt                                       <  1:888917/100‑1 (MQ=255)
                                                            cTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGtt                                      >  1:924159/1‑100 (MQ=255)
                                                                    cTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAcc                              <  2:2422065/100‑1 (MQ=255)
                                                                    cTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAcc                              <  2:2164291/100‑1 (MQ=255)
                                                                    cTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAc                               <  2:26861/99‑1 (MQ=255)
                                                                     tGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAccc                             <  2:1672826/100‑1 (MQ=255)
                                                                      ggCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCg                            >  2:1100298/1‑100 (MQ=255)
                                                                       gCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGt                           >  2:1723334/1‑100 (MQ=255)
                                                                        cTTGATAAATAA‑CCGAATGCCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGtt                          <  2:272191/100‑1 (MQ=255)
                                                                         ttGATAACTAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTg                         >  2:1355773/1‑100 (MQ=255)
                                                                           gATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCa                       <  1:498980/100‑1 (MQ=255)
                                                                             taaataa‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAac                     >  2:2173808/1‑100 (MQ=255)
                                                                             taaataa‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAac                     <  1:2001158/100‑1 (MQ=255)
                                                                             taaataa‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAac                     <  1:1464773/100‑1 (MQ=255)
                                                                              aaataa‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAaca                    <  1:751614/100‑1 (MQ=255)
                                                                                gatctccGAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAacac                   <  2:2198341/95‑1 (MQ=37)
                                                                                  aa‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTc                >  1:1625459/1‑100 (MQ=255)
                                                                                  aa‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTc                >  2:318595/1‑100 (MQ=255)
                                                                                  aa‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTc                <  2:1209778/100‑1 (MQ=255)
                                                                                   ccccGAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTc                >  2:799084/3‑100 (MQ=37)
                                                                                   a‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTcc               >  1:2574916/1‑100 (MQ=255)
                                                                                    cccGAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTcc               >  2:249970/2‑100 (MQ=38)
                                                                                     ccGAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAacacagatc              >  1:2198326/1‑95 (MQ=37)
                                                                                     ccGAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGAGTCACCACGTTGCATCACCCGTTGCAACACCTCCg              >  1:1160877/1‑100 (MQ=38)
                                                                                      cGAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGc             >  2:700868/1‑100 (MQ=38)
                                                                                      cGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGc             >  2:2487674/1‑100 (MQ=255)
                                                                                       gAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCt            <  1:933831/100‑1 (MQ=255)
                                                                                         aTGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCtt           >  1:56276/1‑99 (MQ=255)
                                                                                          tGGCGTCAATAGCGCCGCCATTCGGTGAATTTACCCCTGTTTTCTCATGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCtt           <  1:1249705/98‑1 (MQ=255)
                                                                                           ggcggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCGGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCa        >  1:1658721/1‑100 (MQ=255)
                                                                                           ggcggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTcc         <  1:2426994/99‑1 (MQ=255)
                                                                                           ggcggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCa        >  2:1534895/1‑100 (MQ=255)
                                                                                           ggcggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCa        <  2:1224312/100‑1 (MQ=255)
                                                                                           ggcggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCa        >  1:1524354/1‑100 (MQ=255)
                                                                                            gcggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACCCCTCCGCTTCCAg       >  2:2525007/1‑100 (MQ=255)
                                                                                            gcggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAg       >  1:706233/1‑100 (MQ=255)
                                                                                            gcggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAg       <  1:336882/100‑1 (MQ=255)
                                                                                             cggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGc      >  2:1839396/1‑100 (MQ=255)
                                                                                               gcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTc    >  2:16538/1‑100 (MQ=255)
                                                                                                 aaTAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCg  <  2:725300/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
CGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCG  >  NC_000913/993158‑993353

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