Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,161,946 A→G 100% A161A (GCA→GCG iceT → putative citrate/iron‑citrate/zinc‑citrate efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,161,9460AG95.9% 322.7 / 1.8 98A161A (GCA→GCGiceTputative citrate/iron‑citrate/zinc‑citrate efflux transporter
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base G (38/56);  total (39/59)
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.82e-01

GAGCAATATATGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGC  >  NC_000913/2161851‑2162043
                                                                                               |                                                                                                 
gAGCAATATATGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGcgtcgt                                                                                               <  1:2348866/100‑1 (MQ=255)
 aGCAATATCTGGCGGTGATGACCTTTGTCACGTTAGCCGGGCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTg                                                                                              <  2:2508056/100‑1 (MQ=255)
 aGCAATATATGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTg                                                                                              >  1:977692/1‑100 (MQ=255)
  gCAATATATGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTgg                                                                                             <  1:2463104/100‑1 (MQ=255)
  gCAATATATGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTgg                                                                                             >  1:1572721/1‑100 (MQ=255)
   cAATATATGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGc                                                                                            <  2:2070628/100‑1 (MQ=255)
   cAATATATGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGc                                                                                            <  1:1701434/100‑1 (MQ=255)
     atatatGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCAc                                                                                          <  1:1061454/100‑1 (MQ=255)
     atatatGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCAc                                                                                          >  1:462259/1‑100 (MQ=255)
       atatGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTg                                                                                        <  1:2305298/100‑1 (MQ=255)
       atatGGCGGCGATGACATTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCCGCTGGTGGAGTACGCGTCGTGGCACTg                                                                                        >  2:485396/1‑100 (MQ=255)
        tatGGCGGCGCTGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTgg                                                                                       >  2:1072109/1‑100 (MQ=255)
         atGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGa                                                                                      <  1:604268/100‑1 (MQ=255)
         atGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGa                                                                                      <  2:2317004/100‑1 (MQ=255)
         atGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGa                                                                                      <  1:1645663/100‑1 (MQ=255)
          tGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGAt                                                                                     >  1:2162591/1‑100 (MQ=255)
          tGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGAt                                                                                     <  1:1596580/100‑1 (MQ=255)
          tGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGAt                                                                                     >  1:1637417/1‑100 (MQ=255)
            gcggcgATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCt                                                                                   >  2:1673009/1‑100 (MQ=255)
                cgATGACCTTTGTCACGTTACCCGGTCGGTTTGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCttttt                                                                               <  2:185285/100‑1 (MQ=255)
                cgATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCttttt                                                                               <  2:1572710/100‑1 (MQ=255)
                 gATGACCTTTGTCACGTTGCCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCtttt                                                                                <  2:2497416/98‑1 (MQ=255)
                 gATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTg                                                                              <  1:133180/100‑1 (MQ=255)
                 gATGACCTTTGTCACGTTACCCGGGCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTg                                                                              <  2:1610031/100‑1 (MQ=255)
                  aTGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGa                                                                             <  1:1131518/100‑1 (MQ=255)
                   tGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGAt                                                                            <  1:2055790/100‑1 (MQ=255)
                     aCCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCa                                                                          <  1:2248997/100‑1 (MQ=255)
                      ccTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCaa                                                                         <  1:1566655/100‑1 (MQ=255)
                      ccTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCaa                                                                         >  2:163148/1‑100 (MQ=255)
                      ccTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCaa                                                                         >  2:1096244/1‑100 (MQ=255)
                      ccTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCaa                                                                         >  2:1054986/1‑100 (MQ=255)
                      ccTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCaa                                                                         >  2:166458/1‑100 (MQ=255)
                      ccTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCaa                                                                         >  1:2334579/1‑100 (MQ=255)
                      ccTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCa                                                                          <  2:108992/99‑1 (MQ=255)
                        tttGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATTTTTTTGATCAACa                                                                       <  2:1586532/100‑1 (MQ=255)
                           gTCACGTTACCCGGTCAGGTCGGTCAGCTGCTCGGTCCGGCGCTCGGCGGGCTGCTGGTGGAGTACGCGTCGTTGAGCCGGACCTTTTTg                                                                              >  2:2344897/1‑90 (MQ=255)
                            tgccGTTCCCCGGGCAGGGCGGTCCGCTGCTCGGGCCGGCGCGCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACa                                                                       <  2:2313584/93‑1 (MQ=255)
                              aCGTTACGCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCgg                                                                 >  2:1343141/1‑100 (MQ=255)
                              aCGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCgg                                                                 >  2:2491491/1‑100 (MQ=255)
                               cGTTACCCGGTCAGTTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTTGTGGAGTACGCGTCGTGGCACTGGATCTTTTTTATCAACATTCCgg                                                                 <  2:1584381/99‑1 (MQ=255)
                               cGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGt                                                                >  2:599598/1‑100 (MQ=255)
                                   aCCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTgggg                                                            >  2:1718490/1‑100 (MQ=255)
                                   aCCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTgggg                                                            <  1:130742/100‑1 (MQ=255)
                                    cccGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTgggg                                                            >  1:2187173/1‑99 (MQ=255)
                                     ccGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGAt                                                          >  2:1475737/1‑100 (MQ=255)
                                     ccGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGGCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGAt                                                          <  1:1314390/100‑1 (MQ=255)
                                       ggTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTa                                                        <  1:1475748/100‑1 (MQ=255)
                                       ggTCAGGTCGGGGCGCTGGTCGGTGCGGCGCGCGGCGGTCTGGTGGTGGAGTGCGCGTCGGGGCACTGGATCGTTTTGATCAACATTCCGGTggg                                                             <  2:1993516/95‑1 (MQ=255)
                                           aGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGGACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCgg                                                    <  2:2019632/100‑1 (MQ=255)
                                            ggtcggtcCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGt                                                   >  2:1081686/1‑100 (MQ=255)
                                             gtcggtcCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACTCGTCGTGGCACTGGATCTTTTTGAGCAACATTCCGGTGGGGATTATCGGTg                                                  <  1:1826515/100‑1 (MQ=255)
                                             gtcggtcCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTg                                                  >  1:137489/1‑100 (MQ=255)
                                             gtcggtcCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTg                                                  <  1:2187510/100‑1 (MQ=255)
                                              tcggtcCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGc                                                 <  1:1007098/100‑1 (MQ=255)
                                                ggtcCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGa                                               >  1:1145235/1‑100 (MQ=255)
                                                  tcCGCTGCGCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCg                                                <  2:749832/97‑1 (MQ=255)
                                                   ccGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCg                                            >  2:763173/1‑100 (MQ=255)
                                                     gctgctCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGcc                                          <  1:1527913/100‑1 (MQ=255)
                                                     gctgctCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGcc                                          <  1:358439/100‑1 (MQ=255)
                                                        gctCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCcaca                                       <  2:2063661/100‑1 (MQ=255)
                                                         ctCGGTCCGGCGCTCGGCTGTCTGCTGGTGGAGTACGCGTCGTGGCACCGGATCTTGTTGAGCAACATTCCGGTGGGTATTATCGGTGCGATCGCcactt                                      <  2:500945/100‑3 (MQ=255)
                                                           cGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTg                                    <  1:651180/100‑1 (MQ=255)
                                                           cGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTg                                    >  1:38065/1‑100 (MQ=255)
                                                             gTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCt                                  >  2:1854276/1‑100 (MQ=255)
                                                             gTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCt                                  >  1:2502650/1‑100 (MQ=255)
                                                              tCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTg                                 <  1:1525454/100‑1 (MQ=255)
                                                              tCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTg                                 <  2:2159056/100‑1 (MQ=255)
                                                               ccGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGt                                <  1:493983/100‑1 (MQ=255)
                                                                 ggCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTa                              <  1:1956781/100‑1 (MQ=255)
                                                                     cTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGc                          <  1:722185/100‑1 (MQ=255)
                                                                        ggcggTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGa                       <  2:2162577/100‑1 (MQ=255)
                                                                        ggcggTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGa                       <  1:231587/100‑1 (MQ=255)
                                                                           ggTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACt                    <  1:1718502/100‑1 (MQ=255)
                                                                           ggTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACt                    <  1:1545996/100‑1 (MQ=255)
                                                                            gTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTa                   >  1:656527/1‑100 (MQ=255)
                                                                              ctgctgGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTaca                 >  1:714197/1‑100 (MQ=255)
                                                                              ctgctgGTGGAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTaca                 <  1:581129/100‑1 (MQ=255)
                                                                               tgctgGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTacac                <  1:1589791/100‑1 (MQ=255)
                                                                                gctgGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACAcc               <  1:763180/100‑1 (MQ=255)
                                                                                gctgGTGGAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTacaag               >  2:581122/1‑98 (MQ=255)
                                                                                 ctgGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCa              <  1:499747/100‑1 (MQ=255)
                                                                                 ctgGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCa              >  1:1119196/1‑100 (MQ=255)
                                                                                   ggtggAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATg            >  1:298052/1‑100 (MQ=255)
                                                                                   ggtggAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATg            <  1:1816784/100‑1 (MQ=255)
                                                                                   gggggAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTCTCGGTGCGATCGCCACATTGCTGTTAATGCCGCACTACACCCTg            >  2:2208094/4‑100 (MQ=255)
                                                                                    gtggAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCTATCGCCACATTGCTGTTAACGGCGAACTACACCATGc           <  2:1074151/100‑1 (MQ=255)
                                                                                    gtggAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGc           <  2:1428765/100‑1 (MQ=255)
                                                                                    gtggAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGc           <  1:1148959/100‑1 (MQ=255)
                                                                                     tggAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGTTTCGATTGCTGAATCGGTGTTCATCCAGAACTACACCATGCa          <  2:903659/100‑1 (MQ=255)
                                                                                     tggAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCa          <  1:1744859/100‑1 (MQ=255)
                                                                                      ggAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAg         >  1:2070989/1‑100 (MQ=255)
                                                                                      ggAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAg         >  1:654363/1‑100 (MQ=255)
                                                                                       gAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGa        <  2:1901627/100‑1 (MQ=255)
                                                                                       gAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGa        >  2:930464/1‑100 (MQ=255)
                                                                                        aGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGAc       >  1:106166/1‑100 (MQ=255)
                                                                                        aGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGAc       >  1:111325/1‑100 (MQ=255)
                                                                                         gTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCAACATTGCTGTTAATGCCGAACTCCACCATGCAGAcg      >  2:1188467/1‑100 (MQ=255)
                                                                                            cGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACgcgg   <  1:232682/100‑1 (MQ=255)
                                                                                             gCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACgcggc  <  1:590342/100‑1 (MQ=255)
                                                                                               |                                                                                                 
GAGCAATATATGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGC  >  NC_000913/2161851‑2162043

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: