Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,344,036 G→C 16.7% T104S (ACT→AGT)  yciT ← global regulator of transcription; DeoR family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,344,0360GC16.7% 233.2 / 42.1 96T104S (ACT→AGT) yciTglobal regulator of transcription; DeoR family
Reads supporting (aligned to +/- strand):  ref base G (41/39);  new base C (8/8);  total (49/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01

CCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCG  >  NC_000913/1343937‑1344135
                                                                                                   |                                                                                                   
ccACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCagag                                                                                                     <  1:2541759/100‑1 (MQ=255)
 cACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGt                                                                                                    >  2:739948/1‑100 (MQ=255)
   ccGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTcc                                                                                                  >  1:575175/1‑100 (MQ=255)
    cGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCg                                                                                                 >  1:2215287/1‑100 (MQ=255)
       gCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGc                                                                                              >  1:23212/1‑100 (MQ=255)
        cAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCg                                                                                             <  1:447560/100‑1 (MQ=255)
           aaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTACTCGCCCGGACTCTGGGCGAgc                                                                                          <  1:1927887/100‑2 (MQ=255)
           aaTAACTGCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGt                                                                                          <  1:1141617/100‑1 (MQ=255)
            aTAACTTCACACGGCGCGTCTTTCAGCAAATGCGCTATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTACTCGCCCGGACTCTGGGCGAgc                                                                                          <  1:2139693/99‑2 (MQ=38)
            aTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCGGACTCTGGGCGAgca                                                                                         <  1:1425379/100‑3 (MQ=255)
            aTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCGGACTCTGGGCGAgca                                                                                         <  1:439496/100‑3 (MQ=255)
            aTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTa                                                                                         <  2:2567163/100‑1 (MQ=255)
            aTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTa                                                                                         <  2:1026784/100‑1 (MQ=255)
             tAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCGGACTCTGGGCGAGCAg                                                                                        <  2:2559945/100‑1 (MQ=255)
             tAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAg                                                                                        >  2:1654776/1‑100 (MQ=255)
             tAACGTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCCCGCCCGGACTCTGGGCGAGCAg                                                                                        <  1:1923501/100‑1 (MQ=255)
             tAAATTAACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATGGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAg                                                                                        <  1:2062643/100‑1 (MQ=255)
              aaCTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAgg                                                                                       >  2:146366/1‑100 (MQ=255)
              aaCTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTACTCGCCCGGACTCTGGGCGAGCAga                                                                                       <  1:2493567/100‑2 (MQ=38)
               aCTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAggg                                                                                      <  2:887090/100‑1 (MQ=255)
                cgtCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGc                                                                                     <  1:1690293/98‑1 (MQ=255)
                cTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGc                                                                                     >  2:1955562/1‑100 (MQ=255)
                cTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTACTCGCCCGGACTCTGGGCGAgcagaag                                                                                     <  2:1704315/100‑7 (MQ=37)
                 ttCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCa                                                                                    >  1:270760/1‑100 (MQ=255)
                   cacaAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCAtt                                                                                  >  2:604392/1‑100 (MQ=255)
                    acaAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtg                                                                                 >  2:1656742/1‑100 (MQ=255)
                     caAGGCGCGTCTTTCAGCACATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgc                                                                                >  2:908475/1‑100 (MQ=255)
                      aaGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCATGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgc                                                                                <  1:1812478/99‑1 (MQ=255)
                       aGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgctg                                                                              <  1:2438302/100‑1 (MQ=255)
                        ggCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgctgc                                                                             <  1:1972822/100‑1 (MQ=255)
                        ggCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgctgc                                                                             <  1:2201487/100‑1 (MQ=255)
                         gcgcGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgctgct                                                                            >  2:201115/1‑100 (MQ=255)
                          cgcgTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgctgctg                                                                           <  1:983515/100‑1 (MQ=255)
                          cgcgTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgctgctg                                                                           <  2:342677/100‑1 (MQ=255)
                           gcgTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCGTCAGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgctgctgc                                                                          >  1:1814882/1‑100 (MQ=255)
                               cTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCAtt                                                                      >  2:2464609/1‑100 (MQ=255)
                                tttCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACAGTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCAttt                                                                     <  1:912023/100‑1 (MQ=255)
                                  tCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGGTAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTc                                                                   <  1:190982/100‑1 (MQ=255)
                                  tCAGCAAATGCGCGAGGTAGCTGCTGACCGTGAGGATAGTGACATTTTTCTTCTGCTCGCCTAGCGTCCGGGCGAGTAGGGCATTGCTGCTGCCAtttgc                                                                   >  1:2406084/1‑98 (MQ=255)
                                   cAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCg                                                                  >  2:948788/1‑100 (MQ=255)
                                    aGCAAATGTGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCg                                                                  >  1:2290291/1‑99 (MQ=255)
                                      cAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATa                                                               >  2:2051901/1‑100 (MQ=255)
                                       aaaTGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATaa                                                              >  2:2388961/1‑100 (MQ=255)
                                        aaTGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATaaa                                                             <  1:1461667/100‑1 (MQ=255)
                                          tGCGCGATGTAGCTGCTCACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCACAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTCTCGATAACGa                                                           >  1:403973/1‑100 (MQ=255)
                                           gcgcGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGAt                                                          <  2:735696/100‑1 (MQ=255)
                                            cgcgATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATg                                                         >  1:560907/1‑100 (MQ=255)
                                                 tGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTc                                                    >  1:1687627/1‑100 (MQ=255)
                                                 tGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTc                                                    <  1:2023923/100‑1 (MQ=255)
                                                       tgctgACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCgg                                              <  1:44894/100‑1 (MQ=255)
                                                       tgctgACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCgg                                              >  2:644558/1‑100 (MQ=255)
                                                        gctgACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCgg                                              >  1:1914111/1‑99 (MQ=255)
                                                        gctgACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGt                                             >  1:1553380/1‑100 (MQ=255)
                                                          tgACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTg                                           <  1:1934876/100‑1 (MQ=255)
                                                           gACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGa                                          >  1:1665495/1‑100 (MQ=255)
                                                            aCCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGaa                                         >  1:683198/1‑100 (MQ=255)
                                                            aCCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGaa                                         <  2:1451468/100‑1 (MQ=255)
                                                             ccGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAAc                                        >  2:371391/1‑100 (MQ=255)
                                                                   tgatAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGa                                  <  1:902854/100‑1 (MQ=255)
                                                                   tgatAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGa                                  >  2:159900/1‑100 (MQ=255)
                                                                      tAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGAcgc                               <  1:2479193/100‑1 (MQ=255)
                                                                       aGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGAcgcg                              <  1:1218186/100‑1 (MQ=255)
                                                                        gTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCgg                             >  2:1732283/1‑100 (MQ=255)
                                                                         tGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGc                            >  1:2397671/1‑100 (MQ=255)
                                                                          gACATTTTTCGTCTGCTCGCCCAGAGTCCGAGCGATGAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGc                            <  1:2020613/99‑1 (MQ=255)
                                                                           aCATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCaa                          <  2:1871328/100‑1 (MQ=255)
                                                                            cATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCaaa                         >  2:2183068/1‑100 (MQ=255)
                                                                            cATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCaaa                         >  1:106675/1‑100 (MQ=255)
                                                                                  tcttctGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGc                   <  1:877330/100‑1 (MQ=255)
                                                                                  tcttctGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGc                   <  1:131517/100‑1 (MQ=255)
                                                                                   gccctacTCGCCCGGACTCTGGGCGAGCAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCAACGCGGCAAACTCGGCg                  >  2:2196992/7‑100 (MQ=37)
                                                                                   cttctGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCg                  <  1:1875567/100‑1 (MQ=255)
                                                                                     tctGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAg                <  2:122557/100‑1 (MQ=255)
                                                                                     tctGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAg                <  2:1687636/100‑1 (MQ=255)
                                                                                     cctacTCGCCCGGACTCTGGGCGAGCAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAg                >  2:1965942/5‑100 (MQ=38)
                                                                                       tGCTCGCCCCGAGCTCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGACAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGtt              <  1:185526/100‑1 (MQ=255)
                                                                                       tGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGtt              >  2:648639/1‑100 (MQ=255)
                                                                                       tGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGtt              >  2:539890/1‑100 (MQ=255)
                                                                                       tGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGtt              <  2:395035/100‑1 (MQ=255)
                                                                                        gCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTc             >  2:1807338/1‑100 (MQ=255)
                                                                                        acTCGCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTc             >  1:1345734/2‑100 (MQ=255)
                                                                                        acTCGCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTc             >  1:388067/2‑100 (MQ=255)
                                                                                        acTCGCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTc             >  1:1064305/2‑100 (MQ=255)
                                                                                        acTCGCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCCGCAAACTCGGCGAGTTc             >  1:1611310/2‑100 (MQ=255)
                                                                                         cTCGCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCa            >  2:175197/1‑100 (MQ=255)
                                                                                         cTCGCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCa            >  2:1371009/1‑100 (MQ=255)
                                                                                         cTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCa            >  2:2374513/1‑100 (MQ=255)
                                                                                         cTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCa            <  2:575177/100‑1 (MQ=255)
                                                                                          tCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACGCGGCGAGTTCAc           <  1:908472/100‑1 (MQ=255)
                                                                                           cGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACg          >  1:1462900/1‑100 (MQ=255)
                                                                                           cGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACg          >  2:835260/1‑100 (MQ=255)
                                                                                             cccAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCt        >  2:2141518/1‑100 (MQ=255)
                                                                                             cccAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCt        <  2:1553389/100‑1 (MQ=255)
                                                                                             cccAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCt        <  2:1696721/100‑1 (MQ=255)
                                                                                              ccAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCtt       >  2:1929325/1‑100 (MQ=255)
                                                                                               cAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTc      <  1:2186619/100‑1 (MQ=255)
                                                                                               cAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTc      <  1:1798689/100‑1 (MQ=255)
                                                                                                 gagTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAg    >  2:1062894/1‑100 (MQ=255)
                                                                                                   gTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCg  >  1:1866085/1‑100 (MQ=255)
                                                                                                   gTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCg  <  2:601494/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCG  >  NC_000913/1343937‑1344135

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: