Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 741,817 C→T 20.4% G247D (GGC→GAC)  dtpD ← dipeptide and tripeptide permease D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913741,8170CT20.4% 200.7 / 42.2 98G247D (GGC→GAC) dtpDdipeptide and tripeptide permease D
Reads supporting (aligned to +/- strand):  ref base C (35/43);  new base T (9/11);  total (44/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

ACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACCAGCAGAACCAGCAGCCA  >  NC_000913/741719‑741914
                                                                                                  |                                                                                                 
acaacaTACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGcc                                                                                                  <  2:873441/100‑1 (MQ=255)
    caTACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCaata                                                                                              <  1:2073556/100‑1 (MQ=255)
    caTACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCaata                                                                                              <  1:1327049/100‑1 (MQ=255)
    caTACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCaata                                                                                              <  2:134076/100‑1 (MQ=255)
     aTACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTTCCAGTACACCCAGGCCaataa                                                                                             >  2:1679157/1‑100 (MQ=255)
     aTACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCaataa                                                                                             <  2:958004/100‑1 (MQ=255)
     aTACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCaataa                                                                                             <  1:2127525/100‑1 (MQ=255)
     aTACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCaataa                                                                                             <  1:204910/100‑1 (MQ=255)
     aTACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCaata                                                                                              <  1:2520532/99‑1 (MQ=255)
      tACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCaataat                                                                                            >  1:2293294/1‑100 (MQ=255)
       aCTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGTCTTTTCTCCCAGGCCAataata                                                                                           >  1:621855/1‑100 (MQ=255)
        cTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAg                                                                                          >  1:2259024/1‑100 (MQ=255)
        cTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAg                                                                                          >  2:1009843/1‑100 (MQ=255)
         tGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGt                                                                                         >  2:1626197/1‑100 (MQ=255)
           aagaagGTGAGCGTCACAATCTGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTcg                                                                                       <  1:1518035/100‑1 (MQ=255)
               agGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGAc                                                                                   >  2:1986229/1‑100 (MQ=255)
               agGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGGTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGAc                                                                                   <  1:963632/100‑1 (MQ=255)
                   gAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTTCCATTCTACCCAGGCCAATAATAGTCGCGACa                                                                                  >  1:324495/1‑97 (MQ=255)
                     gCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTaa                                                                             <  1:307980/100‑1 (MQ=255)
                       gTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAgg                                                                           <  1:81786/100‑1 (MQ=255)
                        tCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGc                                                                          >  1:1805059/1‑100 (MQ=255)
                         cacaATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCg                                                                         >  2:286978/1‑100 (MQ=255)
                         cacaATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCg                                                                         <  2:853803/100‑1 (MQ=255)
                           caATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGTCAATAATAGTCGCGACAATTAAGGCGta                                                                       >  1:1049471/1‑100 (MQ=255)
                             aTCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGtata                                                                     <  1:166958/100‑1 (MQ=255)
                               cAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGTCAATAATAGTCGCGACAATTAAGGCGTATAcc                                                                   >  2:1190345/1‑100 (MQ=255)
                                 gCCCCAGCTCCTTCCGCTGTTGCTGGTGTTCTGCTTTGCGATAAAGTGTTGCCAGTACACCCAGGCGAGTAATAGTCGCGACATTTAAGGCGTATACCGa                                                                 >  1:1509369/1‑100 (MQ=255)
                                  ccccAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGAc                                                                <  2:1612360/100‑1 (MQ=255)
                                   cccAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGTCAATAATAGTCGCGACAATTAAGGCGTATACCGAcc                                                               <  1:2550178/100‑1 (MQ=255)
                                   cccAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGAcc                                                               <  1:1327109/100‑1 (MQ=255)
                                   cccAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGAcc                                                               <  2:137489/100‑1 (MQ=255)
                                     cAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCAc                                                             >  1:516891/1‑100 (MQ=255)
                                          ccTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCttt                                                        <  1:2540752/100‑1 (MQ=255)
                                             tCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTTTTTTTTTTTTCAGGCCAATAATAGTCGCGACAATTAACTCGTATACCGACCACTCTTTCCa                                                     >  1:398196/1‑100 (MQ=255)
                                             tCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCa                                                     >  2:460026/1‑100 (MQ=255)
                                              ccGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAg                                                    <  2:711079/100‑1 (MQ=255)
                                                 cTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGTCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAAc                                                 <  2:2220693/100‑1 (MQ=255)
                                                 cTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAAc                                                 <  1:568807/100‑1 (MQ=255)
                                                   gTTTCTGGCGTTCTGCTTTGCGATCACTTTTTGCCAGTACACCCAGGCCAATAATCGTCGCGACCATGACTGCTTATACCGACCACTCCTTCCAg                                                    <  1:1587509/95‑1 (MQ=255)
                                                      ttttggtttCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTAt                                            <  1:2155960/94‑1 (MQ=255)
                                                      tCTGGTTTTCTTCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTAt                                            <  1:2153440/100‑1 (MQ=255)
                                                      tCTGGTTTTCTGCTTTGCGATAAATTTTTTTTTTTTTTTTCAGGTCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACATTa                                             >  1:1769461/1‑99 (MQ=255)
                                                       cTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATg                                           >  2:730929/1‑100 (MQ=255)
                                                       cTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATg                                           >  1:743033/1‑100 (MQ=255)
                                                        tGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGTCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATgg                                          <  2:1211112/100‑1 (MQ=255)
                                                         ggTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGTTCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGt                                         <  1:2512848/100‑1 (MQ=255)
                                                         ggTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGt                                         <  2:1173756/100‑1 (MQ=255)
                                                         ggTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGt                                         <  2:890750/100‑1 (MQ=255)
                                                         ggTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGt                                         >  2:1020169/1‑100 (MQ=255)
                                                          gTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGTCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTa                                        <  2:2501307/100‑1 (MQ=255)
                                                            tttCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACGACTCTTTCCAGAACAGTATgg                                          >  1:2225753/1‑96 (MQ=255)
                                                             ttCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATc                                     <  1:1240724/100‑1 (MQ=255)
                                                             ttCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATc                                     <  2:1102335/100‑1 (MQ=255)
                                                                 gCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTg                                 <  2:2066305/100‑1 (MQ=255)
                                                                     tGCGATAAATTTTTGCCAGTACACCCAGGTCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAgg                             <  1:574637/100‑1 (MQ=255)
                                                                     tGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTACGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAgg                             <  2:1921990/100‑1 (MQ=255)
                                                                      gCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGTGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGc                            >  1:723605/1‑100 (MQ=255)
                                                                       cGATAAATTTTTGCCAGTACACCCAGGTCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCg                           >  1:2541429/1‑100 (MQ=255)
                                                                       cGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAATGCAGGCg                           >  1:782391/1‑100 (MQ=255)
                                                                        gATAAATTTTTGCCAGTACACCCAGGTCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCg                           >  1:938697/1‑99 (MQ=255)
                                                                        gATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGt                          >  2:1201521/1‑100 (MQ=255)
                                                                        gATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGt                          >  1:1617840/1‑100 (MQ=255)
                                                                        gATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGt                          >  1:482874/1‑100 (MQ=255)
                                                                         aTAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGtt                         >  2:2409410/1‑100 (MQ=255)
                                                                         aTAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGtt                         >  2:1413885/1‑100 (MQ=255)
                                                                          tAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACATTATGGTAATCAGTGCAGGCGTTg                        >  2:2107794/1‑100 (MQ=255)
                                                                          tAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTg                        >  2:1230747/1‑100 (MQ=255)
                                                                           aaaTTTTTGCCAGTACACCCAGGTCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGc                       >  2:394122/1‑100 (MQ=255)
                                                                           aaaTTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGc                       <  2:2517765/100‑1 (MQ=255)
                                                                             aTTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGcca                     >  2:2438402/1‑100 (MQ=255)
                                                                             aTTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGcca                     >  2:1952771/1‑100 (MQ=255)
                                                                             aTTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGcca                     <  2:786525/100‑1 (MQ=255)
                                                                              tttttGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGccac                    <  2:1405068/100‑1 (MQ=255)
                                                                               ttttGCCAGTACACCCATGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGcca                     <  1:795650/98‑1 (MQ=255)
                                                                               ttttGCCAGTACACCCAGGTCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGccacc                   <  1:485857/100‑1 (MQ=255)
                                                                               ttttGCCAGTACACCCAGGTCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGccacc                   <  2:1907230/100‑1 (MQ=255)
                                                                               ttttGCCAGTACACCCAGGTCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGccacc                   <  2:250560/100‑1 (MQ=255)
                                                                                  tGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCTAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACcagc                <  1:1212876/100‑1 (MQ=255)
                                                                                   gCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCTCTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACcagca               >  1:1871574/1‑100 (MQ=255)
                                                                                    ccAGTACACCCAGGTCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACcagcag              >  1:725067/1‑100 (MQ=255)
                                                                                     cAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACCAGCAGa             <  1:2251347/100‑1 (MQ=255)
                                                                                      aGTACACCCAGGTCAATAATAGGCGCGACAATTAAGGCGTATAGCGACCACTCTTGCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACCAGCAGaa            <  2:760003/100‑1 (MQ=255)
                                                                                      aGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACCAGCAGaa            >  1:2012103/1‑100 (MQ=255)
                                                                                      aGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTAGACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGAGTTGCCACCAGCAGaa            >  1:1412617/1‑100 (MQ=255)
                                                                                         acacCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACCAGCAGAACca         <  1:1059306/100‑1 (MQ=255)
                                                                                          cacCCAGGCCAATAATAGTCGCGACAATTATGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACCAGCAGAACcag        <  1:548230/100‑1 (MQ=255)
                                                                                          cacCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACCAGCAGAACcag        >  1:1312139/1‑100 (MQ=255)
                                                                                          cacCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACCAGCAGAACcag        <  2:1278652/100‑1 (MQ=255)
                                                                                          cacCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACCAGCAGAACcag        >  1:1312190/1‑100 (MQ=255)
                                                                                          cacCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACCAGCAGAACcag        <  2:2228167/100‑1 (MQ=255)
                                                                                           acCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACCAGCAGAACcagc       >  1:1858607/1‑100 (MQ=255)
                                                                                            cccAGGTCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACCAGCAGAACcagca      >  1:1558765/1‑100 (MQ=255)
                                                                                            cccAGGTCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACCAGCAGAACcagca      >  1:1558084/1‑100 (MQ=255)
                                                                                             ccAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACCAGCAGAACcagcag     <  1:1454769/100‑1 (MQ=255)
                                                                                              cAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACCAGCAGAACCagcagc    <  2:2292026/100‑1 (MQ=255)
                                                                                              cAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACCAGCAGAACCagcagc    <  2:2292018/100‑1 (MQ=255)
                                                                                               aGGTCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACCAGCAGAACCAGCAGcc   <  1:740174/100‑1 (MQ=255)
                                                                                               aGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACCAGCAGAACCAGCAGcc   <  1:168505/100‑1 (MQ=255)
                                                                                                ggCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACCAGCAGAACCAGCAGCCa  <  1:1952761/100‑1 (MQ=255)
                                                                                                  |                                                                                                 
ACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTCGCGACAATTAAGGCGTATACCGACCACTCTTTCCAGAACAGTATGGTAATCAGTGCAGGCGTTGCCACCAGCAGAACCAGCAGCCA  >  NC_000913/741719‑741914

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: