Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,666,470 C→A 18.9% L386I (CTA→ATA)  clcB → H(+)/Cl(‑) exchange transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,666,4700CA18.9% 392.2 / 77.0 158L386I (CTA→ATA) clcBH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base C (60/68);  new base A (15/15);  total (75/83)
Fisher's exact test for biased strand distribution p-value = 8.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 162 total aligned reads displayed.

TTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGC  >  NC_000913/1666371‑1666568
                                                                                                   |                                                                                                  
ttccTCGTATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGc                                                                                                    <  1:3489459/97‑1 (MQ=255)
   ntCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                                 >  1:2929800/2‑100 (MQ=255)
    tCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACt                                                                                                >  1:372749/1‑100 (MQ=255)
     cGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACt                                                                                                <  1:1174564/99‑1 (MQ=255)
        aTTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTcccc                                                                                            <  1:3904636/100‑1 (MQ=255)
        aTTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTcccc                                                                                            <  2:1262718/100‑1 (MQ=255)
         ttGTCCGGGATGGCGCCACTGTTGGCGGCGACCACGCATGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCg                                                                                           <  1:3572558/100‑1 (MQ=255)
          tGACTGGGATCCCGACACTGTTGTCGTCGACCACGCACGCGCCGATGATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                          <  1:3900869/100‑1 (MQ=37)
          tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCgg                                                                                          <  1:1155517/100‑1 (MQ=255)
          tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCgg                                                                                          <  1:2327048/100‑1 (MQ=255)
          tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCgg                                                                                          <  2:373770/100‑1 (MQ=255)
          tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                          <  2:3656435/100‑1 (MQ=39)
          tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                          <  1:3228419/100‑1 (MQ=39)
          tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                          <  1:2781912/100‑1 (MQ=39)
          tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                          <  2:4381550/100‑1 (MQ=255)
          tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                          <  2:830051/100‑1 (MQ=39)
           gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                          <  2:3999520/99‑1 (MQ=39)
            aCCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGtt                                                                                        >  2:232772/1‑100 (MQ=255)
            aCCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGtt                                                                                        <  2:1419121/100‑1 (MQ=255)
            aCCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGtt                                                                                        <  2:516661/100‑1 (MQ=255)
             ccGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGGAGATGATGATATGTGAGATGACCGGGGGGAATCAGCTACTCCCCGGgtt                                                                                       >  1:2969928/1‑100 (MQ=255)
             ccGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGttt                                                                                       >  2:3148576/1‑100 (MQ=255)
             ccGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGttt                                                                                       <  1:2238951/100‑1 (MQ=255)
             ccGGGATGGCGACACTGTTGGCGGCGACCACACACGCGCCGATTATGTCGACGTTTATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGtc                                                                                        <  1:1347488/99‑2 (MQ=38)
               gggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTat                                                                                     >  2:1904952/1‑100 (MQ=255)
                ggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTatt                                                                                    <  2:3414531/100‑1 (MQ=255)
                ggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTatt                                                                                    <  2:1296003/100‑1 (MQ=255)
                ggATGGCGACACTGTCGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGtcattt                                                                                    <  2:1210423/100‑6 (MQ=37)
                 gATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattg                                                                                   <  2:4339634/100‑1 (MQ=255)
                 gATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattg                                                                                   <  2:1600393/100‑1 (MQ=255)
                 gATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattg                                                                                   <  2:2929816/100‑1 (MQ=255)
                     gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgattg                                                                               <  1:484066/100‑1 (MQ=255)
                         cacTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTTATGATATGTGATATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTg                                                                           <  1:1678967/100‑1 (MQ=255)
                         aacTGTTGGCGGAGACCACGCATGCGCCGATTATGTCGGCTTTGATGTTATGTGAAATGTCCGGGGAGTATCGGCTACTCCTCGGTTTATAGATTGCCt                                                                            >  1:3041537/2‑99 (MQ=255)
                           cTGTTGGCGGCGACCACGCCCGCGCCGTTTATGTCGACTTTGATGATATGTGAAATGACCGGGGAGTATCAGCTATTCCCCGGTTTATTGATTGCCTGCg                                                                         <  1:2863512/100‑1 (MQ=255)
                           cTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCg                                                                         >  1:952531/1‑100 (MQ=255)
                           cTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCg                                                                         >  2:3907289/1‑100 (MQ=255)
                             ggtGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTa                                                                       >  2:2962568/3‑100 (MQ=255)
                             gTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTa                                                                       <  2:752367/100‑1 (MQ=255)
                               tGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAt                                                                     <  2:1563833/100‑1 (MQ=255)
                               tGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAt                                                                     >  1:3270249/1‑100 (MQ=255)
                                ggcggcGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAtt                                                                    >  2:2236396/1‑100 (MQ=255)
                                ggcggcGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAtt                                                                    <  1:3689149/100‑1 (MQ=255)
                                 gcggcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTg                                                                   >  2:4131384/1‑100 (MQ=255)
                                   ggcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCg                                                                 <  2:1144492/100‑1 (MQ=255)
                                    gcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCGACTCCCCGGTTTATTGATTGCCTGCGTAATTg                                                                   <  1:1444067/97‑1 (MQ=255)
                                     cgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGAAAttaccga                                                               >  2:1294436/1‑95 (MQ=255)
                                      gACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCg                                                              >  2:1556277/1‑100 (MQ=255)
                                       aCCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCgg                                                             <  1:1904946/100‑1 (MQ=255)
                                        ccACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGt                                                            <  2:368654/100‑1 (MQ=255)
                                        ccACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGt                                                            <  2:1244411/100‑1 (MQ=255)
                                         cacgcacgCGCCGATTATGTCGACGTTGATGATCTGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTGTATTGATTGCCTGCGTAATTGAGTCGGTa                                                           >  1:2649881/1‑100 (MQ=255)
                                         cacgcacgCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTa                                                           >  2:2666070/1‑100 (MQ=255)
                                                  gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCgg                                                  <  2:4366550/100‑1 (MQ=255)
                                                   ccGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGa                                                 >  1:3060089/1‑100 (MQ=255)
                                                    cGATTATGTCGACGTTGATGATATGTGGAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGAc                                                <  1:956660/100‑1 (MQ=255)
                                                     gATTATGTCGACGTTGATGATCTGTGAAATGACCGGGGAGTATCAGCTACGCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACg                                               <  1:284443/100‑1 (MQ=255)
                                                     gATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACg                                               <  2:1393697/100‑1 (MQ=255)
                                                         atgcagccGTTGACGATATGTGAAATGAACGGGGAGACTCAGCTACTCCCCGGTTTATAGATTACCTGCGTAAATGCGCGGGCGCTTTCGCGGACGTTa                                            >  1:184615/8‑99 (MQ=255)
                                                          tGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTaca                                          <  1:53565/100‑1 (MQ=255)
                                                            tCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACAcc                                        <  2:274765/100‑1 (MQ=255)
                                                             cGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCg                                       >  1:3044046/1‑100 (MQ=255)
                                                               aCGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTg                                     <  2:3770984/100‑1 (MQ=255)
                                                               aCGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTg                                     <  2:3078912/100‑1 (MQ=255)
                                                                 gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGCTAATTTCGCGGACGTTACACCGTGAc                                   >  1:234745/1‑100 (MQ=255)
                                                                   tgatgaTATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGActc                                 <  2:4085946/100‑1 (MQ=255)
                                                                     atgatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTCGGTTCGCGGACGTTACACCGTGActct                                >  1:4103367/1‑99 (MQ=255)
                                                                        atatGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatct                            <  2:4162573/100‑1 (MQ=255)
                                                                         tatGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatcta                           <  1:2536706/100‑1 (MQ=255)
                                                                         tatGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatcta                           >  2:4751/1‑100 (MQ=255)
                                                                          caatcaataaaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAc                          >  2:2141762/11‑100 (MQ=25)
                                                                           tgtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAcc                         <  2:2350863/100‑1 (MQ=255)
                                                                           tgtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAcc                         >  2:2369273/1‑100 (MQ=255)
                                                                           tgtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAcc                         <  2:4222251/100‑1 (MQ=255)
                                                                            gtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCg                        >  2:4312878/1‑100 (MQ=255)
                                                                             tgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGc                       >  2:3695012/1‑100 (MQ=255)
                                                                              gAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                      <  1:3581193/100‑1 (MQ=255)
                                                                              gAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                      >  1:4072741/1‑100 (MQ=255)
                                                                                aaTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                      >  1:4012901/1‑98 (MQ=255)
                                                                                 aTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGc                   <  1:2420152/100‑1 (MQ=255)
                                                                                  tGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGca                  >  1:2610092/1‑100 (MQ=255)
                                                                                  tGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGca                  >  2:2224451/1‑100 (MQ=255)
                                                                                     ccGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac               >  2:2496939/1‑100 (MQ=255)
                                                                                     ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac               >  1:1032035/1‑100 (MQ=255)
                                                                                     ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac               >  2:1398040/1‑100 (MQ=39)
                                                                                     ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac               >  2:1752756/1‑100 (MQ=39)
                                                                                     ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac               >  1:2447190/1‑100 (MQ=255)
                                                                                     ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac               >  1:3794172/1‑100 (MQ=39)
                                                                                      cGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac               >  1:2206997/1‑99 (MQ=39)
                                                                                       ggggAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg             >  2:1558033/1‑100 (MQ=255)
                                                                                        gggAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAgatcggaag            >  1:3970385/1‑92 (MQ=255)
                                                                                            gTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag        >  2:2648225/1‑100 (MQ=255)
                                                                                            gTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag        >  2:2863548/1‑100 (MQ=255)
                                                                                            gTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag        <  2:860983/100‑1 (MQ=255)
                                                                                            gTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACACTATCTACCGCCAGCACACTgcgcc         >  2:187494/1‑98 (MQ=255)
                                                                                              aTCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGCAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcagca      >  1:1123555/1‑100 (MQ=255)
                                                                                              aTCAGCTACGCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCCCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcagca      >  1:1739011/1‑100 (MQ=255)
                                                                                               tCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCAt     >  2:9340/1‑100 (MQ=255)
                                                                                                cAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATa    <  2:650968/100‑1 (MQ=255)
                                                                                                 aGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAg   <  1:2369262/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
TTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGC  >  NC_000913/1666371‑1666568

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: