Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 993,258 G→T 15.9% intergenic (+24/+19) 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,2580GT15.9% 318.1 / 43.1 113intergenic (+24/+19)pepN/ssuBaminopeptidase N/aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (56/39);  new base T (5/13);  total (61/52)
Fisher's exact test for biased strand distribution p-value = 2.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01

Only 100 of 117 total aligned reads displayed.

CGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTA  >  NC_000913/993161‑993356
                                                                                                  |                                                                                                  
cgcTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAgcg                                                                                                  <  1:1619466/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)
      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                                                                                         <  2:2046208/100‑1 (MQ=38)
         aGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCGCTATTGCCGCCATTc                                                                                         <  1:1856964/100‑1 (MQ=38)
         aGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTc                                                                                         >  1:2357578/1‑100 (MQ=255)
          gTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCGCTATTGCCGCCATTCg                                                                                        <  2:782470/100‑1 (MQ=38)
          gTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCg                                                                                        >  1:710461/1‑100 (MQ=255)
          gTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCg                                                                                        <  2:1055324/100‑1 (MQ=255)
          gTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCg                                                                                        <  2:591502/100‑1 (MQ=255)
          gTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCg                                                                                        >  1:1988632/1‑100 (MQ=255)
            tGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCGCTATTGCCGCCATTCggt                                                                                      <  2:783039/100‑2 (MQ=38)
            tGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCggg                                                                                      >  2:1945187/1‑100 (MQ=255)
            tGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCggg                                                                                      >  2:168556/1‑100 (MQ=255)
            tGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCggg                                                                                      >  2:2042345/1‑100 (MQ=255)
            tGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCggg                                                                                      >  2:1622581/1‑100 (MQ=255)
            tGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCggg                                                                                      >  2:1622414/1‑100 (MQ=255)
             gagaGGGGTTGAAAATCTCTCTGGCTAGCTTTACGAGGAGATAACTAAAGCACTGGCTTGATAATTAA‑GCGAATGGCGGCGCTATTGCCGCCATTCggt                                                                                      <  1:522510/96‑2 (MQ=16)
             gAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCGCTATTGCCGCCATTCggtt                                                                                     <  2:2105918/100‑3 (MQ=37)
             gAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCgggg                                                                                     <  2:2067518/100‑1 (MQ=255)
                   gCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCGCTATTGCCGCCATTCggttatttat                                                                               <  1:2165882/100‑9 (MQ=25)
                   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‑CCGAATGGCGGCGCTATTGCCGCCATTCggttatttat                                                                               <  1:771409/97‑9 (MQ=21)
                      ggAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAccc                                                                            <  2:673006/100‑1 (MQ=255)
                      ggAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAccc                                                                            >  1:2393960/1‑100 (MQ=255)
                      ggAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAccc                                                                            >  2:2131394/1‑100 (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)
                                  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)
                                           gTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTc                                                       >  2:288031/1‑100 (MQ=255)
                                           gTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTc                                                       >  2:2348645/1‑100 (MQ=255)
                                            tACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCa                                                      <  2:1500020/100‑1 (MQ=255)
                                             aCGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAg                                                     <  2:166061/100‑1 (MQ=255)
                                              cGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGa                                                    <  2:739571/100‑1 (MQ=255)
                                              cGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGa                                                    <  2:1988643/100‑1 (MQ=255)
                                              cGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGa                                                    >  2:498982/1‑100 (MQ=255)
                                                agaagaTAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGAtt                                                  <  2:1704807/100‑1 (MQ=255)
                                                   agaTAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcac                                               >  2:593578/1‑100 (MQ=255)
                                                   agaTAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcac                                               >  2:2318579/1‑100 (MQ=255)
                                                      tAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcaccac                                            >  2:1141935/1‑100 (MQ=255)
                                                       aaCTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACg                                           <  2:1941552/100‑1 (MQ=255)
                                                        aCTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGt                                          <  1:888917/100‑1 (MQ=255)
                                                        aCTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGt                                          >  2:503645/1‑100 (MQ=255)
                                                         cTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGtt                                         >  1:924159/1‑100 (MQ=255)
                                                                 cTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAcc                                 <  2:2164291/100‑1 (MQ=255)
                                                                 cTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAcc                                 <  2:2422065/100‑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‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGtt                             <  2:1061196/100‑1 (MQ=255)
                                                                     cTTGATAAATAA‑CCGAATGCCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGtt                             <  2:272191/100‑1 (MQ=255)
                                                                       tGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGc                           >  2:2547291/1‑100 (MQ=255)
                                                                       tGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGc                           <  2:221047/100‑1 (MQ=255)
                                                                       tGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGc                           >  2:2088833/1‑100 (MQ=255)
                                                                        gATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCa                          <  1:498980/100‑1 (MQ=255)
                                                                          taaataa‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAac                        >  2:2173808/1‑100 (MQ=255)
                                                                          taaataa‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAac                        <  1:1464773/100‑1 (MQ=255)
                                                                          taaataa‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAac                        <  1:2001158/100‑1 (MQ=255)
                                                                           aaataa‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAaca                       <  1:751614/100‑1 (MQ=255)
                                                                             gatctccGAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAacac                      <  2:2198341/95‑1 (MQ=37)
                                                                             ataa‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACAcc                     >  1:2521850/1‑100 (MQ=255)
                                                                               aa‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTc                   >  1:1625459/1‑100 (MQ=255)
                                                                               aa‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTc                   <  2:1209778/100‑1 (MQ=255)
                                                                               aa‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTc                   >  2:318595/1‑100 (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)
                                                                                      aTGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCtt              >  1:56276/1‑99 (MQ=255)
                                                                                       tGGCGTCAATAGCGCCGCCATTCGGTGAATTTACCCCTGTTTTCTCATGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCtt              <  1:1249705/98‑1 (MQ=255)
                                                                                        ggcggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCGGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCa           >  1:1658721/1‑100 (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)
                                                                                             cAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTcc      <  2:710464/100‑1 (MQ=255)
                                                                                                tAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCt   >  2:1665258/1‑100 (MQ=255)
                                                                                                  gcgcCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTa  <  1:1825966/99‑1 (MQ=255)
                                                                                                  |                                                                                                  
CGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAA‑CCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTA  >  NC_000913/993161‑993356

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