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.

CGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTA  >  NC_000913/993161‑993356
                                                                                                  |                                                                                                  
cgcTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAgcg                                                                                                  <  1:1619466/100‑1 (MQ=255)
 ncTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAgcgc                                                                                                 >  1:2572534/2‑100 (MQ=255)
  cTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGcgcc                                                                                                >  2:751617/1‑100 (MQ=255)
   tGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGcgccg                                                                                               >  1:1291240/1‑100 (MQ=255)
     gAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGcgccg                                                                                               >  1:676497/1‑98 (MQ=255)
      aaCAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCa                                                                                            >  1:166060/1‑100 (MQ=255)
        cAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATTCCCCGAATGGCGGCGCTATTGCCGCCAt                                                                                           <  1:1793757/100‑1 (MQ=25)
        cAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCAtt                                                                                          <  1:191990/100‑1 (MQ=38)
        cAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCAtt                                                                                          >  1:776937/1‑100 (MQ=255)
        cAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCAtt                                                                                          >  1:390606/1‑100 (MQ=255)
         aGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCATTc                                                                                         <  1:1784254/100‑1 (MQ=38)
         aGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCATTc                                                                                         <  2:2046208/100‑1 (MQ=38)
         aGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCATTc                                                                                         <  1:1856964/100‑1 (MQ=38)
         aGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTc                                                                                         >  1:2357578/1‑100 (MQ=255)
          gTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCATTCg                                                                                        <  2:782470/100‑1 (MQ=38)
          gTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCg                                                                                        >  1:710461/1‑100 (MQ=255)
          gTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCg                                                                                        <  2:1055324/100‑1 (MQ=255)
          gTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCg                                                                                        <  2:591502/100‑1 (MQ=255)
          gTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCg                                                                                        >  1:1988632/1‑100 (MQ=255)
            tGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCATTCggt                                                                                      <  2:783039/100‑2 (MQ=38)
            tGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCggg                                                                                      >  2:1945187/1‑100 (MQ=255)
            tGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCggg                                                                                      >  2:168556/1‑100 (MQ=255)
            tGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCggg                                                                                      >  2:2042345/1‑100 (MQ=255)
            tGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCggg                                                                                      >  2:1622581/1‑100 (MQ=255)
            tGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCggg                                                                                      >  2:1622414/1‑100 (MQ=255)
             gagaGGGGTTGAAAATCTCTCTGGCTAGCTTTACGAGGAGATAACTAAAGCACTGGCTTGATAATTAAGCGAATGGCGGCGCTATTGCCGCCATTCggt                                                                                      <  1:522510/96‑2 (MQ=16)
             gAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCATTCggtt                                                                                     <  2:2105918/100‑3 (MQ=37)
             gAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCgggg                                                                                     <  2:2067518/100‑1 (MQ=255)
                   gCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCATTCggttatttat                                                                               <  1:2165882/100‑9 (MQ=25)
                   gCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTa                                                                               >  1:1680372/1‑100 (MQ=255)
                    cTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAc                                                                              <  1:286240/100‑1 (MQ=255)
                     tGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACGGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAcc                                                                             <  1:869220/100‑1 (MQ=255)
                      ggAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCATTCggttatttat                                                                               <  1:771409/97‑9 (MQ=21)
                      ggAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAccc                                                                            <  2:673006/100‑1 (MQ=255)
                      ggAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAccc                                                                            >  1:2393960/1‑100 (MQ=255)
                      ggAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAccc                                                                            >  2:2131394/1‑100 (MQ=255)
                       gAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAcccc                                                                           >  1:2231131/1‑100 (MQ=255)
                          aaTCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGt                                                                        <  1:1129562/100‑1 (MQ=255)
                             ctctctGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGtttt                                                                     >  2:653489/1‑100 (MQ=255)
                             ctctctGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCCCCGCCATTCGGGGAATTTAcccc                                                                           >  1:2178735/1‑94 (MQ=255)
                                  tGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAg                                                                <  2:1486157/100‑1 (MQ=255)
                                  tGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGCCAATAGCGCCGCCATTCGGGGAATTTACTCCTGTTTtctccg                                                                <  2:1755594/100‑3 (MQ=255)
                                   ggCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAgg                                                               <  2:1664773/100‑1 (MQ=255)
                                   ggCGATCTGTACGAGAAGATAACTAAAGCACTGGCTGGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAgg                                                               >  2:1825976/1‑100 (MQ=255)
                                    gCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGc                                                              >  2:488100/1‑100 (MQ=255)
                                           gTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTc                                                       >  2:288031/1‑100 (MQ=255)
                                           gTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTc                                                       >  2:2348645/1‑100 (MQ=255)
                                            tACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCa                                                      <  2:1500020/100‑1 (MQ=255)
                                             aCGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAg                                                     <  2:166061/100‑1 (MQ=255)
                                              cGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGa                                                    <  2:739571/100‑1 (MQ=255)
                                              cGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGa                                                    <  2:1988643/100‑1 (MQ=255)
                                              cGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGa                                                    >  2:498982/1‑100 (MQ=255)
                                                agaagaTAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGAtt                                                  <  2:1704807/100‑1 (MQ=255)
                                                   agaTAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcac                                               >  2:593578/1‑100 (MQ=255)
                                                   agaTAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcac                                               >  2:2318579/1‑100 (MQ=255)
                                                      tAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcaccac                                            >  2:1141935/1‑100 (MQ=255)
                                                       aaCTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACg                                           <  2:1941552/100‑1 (MQ=255)
                                                        aCTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGt                                          <  1:888917/100‑1 (MQ=255)
                                                        aCTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGt                                          >  2:503645/1‑100 (MQ=255)
                                                         cTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGtt                                         >  1:924159/1‑100 (MQ=255)
                                                                 cTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAcc                                 <  2:2164291/100‑1 (MQ=255)
                                                                 cTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAcc                                 <  2:2422065/100‑1 (MQ=255)
                                                                  tGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAccc                                <  2:1672826/100‑1 (MQ=255)
                                                                   ggCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCg                               >  2:1100298/1‑100 (MQ=255)
                                                                    gCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGt                              >  2:1723334/1‑100 (MQ=255)
                                                                     cTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGtt                             <  2:1061196/100‑1 (MQ=255)
                                                                     cTTGATAAATAACCGAATGCCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGtt                             <  2:272191/100‑1 (MQ=255)
                                                                       tGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGc                           >  2:2547291/1‑100 (MQ=255)
                                                                       tGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGc                           <  2:221047/100‑1 (MQ=255)
                                                                       tGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGc                           >  2:2088833/1‑100 (MQ=255)
                                                                        gATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCa                          <  1:498980/100‑1 (MQ=255)
                                                                          taaataaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAac                        >  2:2173808/1‑100 (MQ=255)
                                                                          taaataaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAac                        <  1:1464773/100‑1 (MQ=255)
                                                                          taaataaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAac                        <  1:2001158/100‑1 (MQ=255)
                                                                           aaataaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAaca                       <  1:751614/100‑1 (MQ=255)
                                                                             gatctccGAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAacac                      <  2:2198341/95‑1 (MQ=37)
                                                                             ataaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACAcc                     >  1:2521850/1‑100 (MQ=255)
                                                                               aaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTc                   >  1:1625459/1‑100 (MQ=255)
                                                                               aaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTc                   <  2:1209778/100‑1 (MQ=255)
                                                                               aaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTc                   >  2:318595/1‑100 (MQ=255)
                                                                                ccccGAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTc                   >  2:799084/3‑100 (MQ=37)
                                                                                aCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTcc                  >  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)
                                                                                                  |                                                                                                  
CGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTA  >  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: