Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 88,208 A→T 100% T61S (ACA→TCA)  cra → transcriptional repressor‑activator for carbon metabolism

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091388,2080AT97.6% 262.6 / ‑2.0 85T61S (ACA→TCA) cratranscriptional repressor‑activator for carbon metabolism
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (34/48);  total (36/49)
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.95e-01

CCGTGTGAGCGACAAAACCGTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCACACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCTTGAACGCCAGGCGCGGCAACGGGGTTA  >  NC_000913/88111‑88302
                                                                                                 |                                                                                              
ccGTGTGAGCGACAAAACCGTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGctca                                                                                              <  1:1606215/100‑1 (MQ=255)
ccGTGTGAGCGACAAAACCGTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGcgca                                                                                              >  2:2173416/1‑100 (MQ=255)
  gtgtGAGCGACAAAACCGTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACg                                                                                            <  2:300547/100‑1 (MQ=255)
   tgtgAGCGACAAAACCGTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGt                                                                                           >  1:1987312/1‑100 (MQ=255)
     tgAGCGACAAAACCGTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTc                                                                                         >  2:1817591/1‑100 (MQ=255)
         cGACAAAACCGTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTAtt                                                                                     <  2:901214/100‑1 (MQ=255)
         cGACAAAACCGTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTAtt                                                                                     <  1:2335079/100‑1 (MQ=255)
          gACAAAACCGTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTg                                                                                    <  2:317870/100‑1 (MQ=255)
           aCAAAACCGTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTgg                                                                                   <  2:2397616/100‑1 (MQ=255)
           aCAAAACCGTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGGTCTATTgg                                                                                   >  2:405604/1‑100 (MQ=255)
            cAAAACCGTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGt                                                                                  <  1:2001575/100‑1 (MQ=255)
             aaaaCCGTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTc                                                                                 <  2:2558573/100‑1 (MQ=255)
               aaCCGTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCtt                                                                               >  1:2485038/1‑100 (MQ=255)
               aaCCGTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCtt                                                                               >  2:1029569/1‑100 (MQ=255)
                agcGTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCGGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTtg                                                                              <  2:484105/98‑1 (MQ=255)
                aCCGTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTtg                                                                              >  1:1458027/1‑100 (MQ=255)
                aCCGTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTtg                                                                              <  1:842771/100‑1 (MQ=255)
                 ccGTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTtgt                                                                             >  1:1671520/1‑100 (MQ=255)
                   gTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGa                                                                           >  1:1586220/1‑100 (MQ=255)
                   gTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGa                                                                           >  1:451523/1‑100 (MQ=255)
                     tGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATc                                                                         <  2:1234906/100‑1 (MQ=255)
                     tGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATc                                                                         <  2:1235319/100‑1 (MQ=255)
                      gAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATcc                                                                        >  2:549787/1‑100 (MQ=255)
                       aaaaaGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATccc                                                                       >  1:2187132/1‑100 (MQ=255)
                        aaaaGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATcccc                                                                      >  1:80411/1‑100 (MQ=255)
                         aaaGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCg                                                                     >  1:957015/1‑100 (MQ=255)
                           aGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGAt                                                                   >  1:2510991/1‑100 (MQ=255)
                              cATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTg                                                                >  1:1224775/1‑100 (MQ=255)
                               aTGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGGTGGACGCACACGTTCTATTGGTCTTGTGATCCCCGATCTgg                                                               <  2:1398720/100‑1 (MQ=255)
                               aTGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTgg                                                               <  2:1141970/100‑1 (MQ=255)
                               aTGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTgg                                                               <  2:806292/100‑1 (MQ=255)
                                tGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGga                                                              <  1:1381073/100‑1 (MQ=255)
                                tGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGga                                                              <  2:1326955/100‑1 (MQ=255)
                                tGGCTGTCGTGGGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCGCACGTTCTATTGGTCTTGTGATCCCCGATCTGga                                                              >  2:1344434/1‑100 (MQ=255)
                                 ggCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGgag                                                             <  1:1258416/100‑1 (MQ=255)
                                 ggCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGgag                                                             <  2:1017270/100‑1 (MQ=255)
                                 ggCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGgag                                                             <  1:2426956/100‑1 (MQ=255)
                                 ggCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGgag                                                             <  1:362395/100‑1 (MQ=255)
                                  gCTGGGCTGCGTGATCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGgaga                                                            <  2:2339580/100‑1 (MQ=255)
                                   cTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGaa                                                           >  1:508069/1‑100 (MQ=255)
                                   cTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGaa                                                           >  2:1223752/1‑100 (MQ=255)
                                    tgtgGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAac                                                          <  1:660898/100‑1 (MQ=255)
                                          gCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGc                                                    <  1:1229721/100‑1 (MQ=255)
                                          gCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGc                                                    <  2:832236/100‑1 (MQ=255)
                                           cGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCt                                                   >  1:2478219/1‑100 (MQ=255)
                                                 cacaATTAGCACCTGAACGCCGTTGCAGCTGGGTTTCGTGCTGGACGCTCACGTTCTATTGTTCTTGTGATCCCCGATCTGGAGAACACCAGCTATAccc                                             <  2:865836/100‑1 (MQ=255)
                                                  acaATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCg                                            <  2:300269/100‑1 (MQ=255)
                                                  acaATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCCTGTGATCCCCGATCTGGAGAACACCAGCTATACCCg                                            >  2:869275/1‑100 (MQ=255)
                                                   caATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGc                                           >  1:622205/1‑100 (MQ=255)
                                                    aaTTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCa                                          <  1:1442617/100‑1 (MQ=255)
                                                    aaTTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTACTGGACGCTCACGTTCTATTGGTCTTGTGATGCCCGATCTGGAGAACACCAGCTATACCCGCa                                          >  2:1989812/1‑100 (MQ=255)
                                                        accaccCGAACGCCGTGGCAGCTGGGCTTCGTGCTTAACTCTCACGGTCTATTGGCCTTGTGAACACCGATCTGGAGAACACCTGCTCTACCCGCAtacc                                      >  2:1077363/1‑97 (MQ=255)
                                                          caccCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTa                                    <  1:1322706/100‑1 (MQ=255)
                                                          caccCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTa                                    <  1:1487047/100‑1 (MQ=255)
                                                           accCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTaa                                   >  1:2158168/1‑100 (MQ=255)
                                                           accCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTaa                                   >  2:1008502/1‑100 (MQ=255)
                                                           accCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTaa                                   <  1:2285078/100‑1 (MQ=255)
                                                            cccGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAAc                                  >  2:359787/1‑100 (MQ=255)
                                                            cccGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCACACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTaa                                   >  2:1929479/1‑99 (MQ=255)
                                                              cGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTa                                <  1:1966121/100‑1 (MQ=255)
                                                                 aCGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACTAGCTATACCCGCATCGCTAACTATCt                             >  2:2378524/1‑100 (MQ=255)
                                                                 aCGCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCt                             <  1:1880812/100‑1 (MQ=255)
                                                                  cGCCGTTGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCtt                            <  1:716246/100‑1 (MQ=255)
                                                                  cGCCGTGGCAGCTGGGCTTCGTGCTTGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCtt                            <  2:1784566/100‑1 (MQ=255)
                                                                   gCCGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCTTg                           <  1:398090/100‑1 (MQ=255)
                                                                    ccGTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCTTGa                          <  1:2074330/100‑1 (MQ=255)
                                                                     cttGGCAGCTGGGCTTCGTGGTGGGCGCTCACGTTCTATTTGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAAc                                  <  2:2087116/89‑1 (MQ=255)
                                                                      gTGGCAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCTTGAAc                        <  2:1198250/100‑1 (MQ=255)
                                                                          cAGCTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCTTGAACGCCa                    >  1:2148854/1‑100 (MQ=255)
                                                                             cTGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCTTGAACGCCAGgc                 <  1:427945/100‑1 (MQ=255)
                                                                              tGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCTTGAACGCCAGgcg                <  2:48131/100‑1 (MQ=255)
                                                                              tGGGCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCTTGAACGCCAGgcg                >  2:731099/1‑100 (MQ=255)
                                                                               gggCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCTTGAACGCCAGgcgc               >  1:113320/1‑100 (MQ=255)
                                                                                 gCTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCTTGAACGCCAGGCGCgg             >  1:1020266/1‑100 (MQ=255)
                                                                                  cTTCGTGCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCTTGAACGCCAGGCGCGGc            <  1:227017/100‑1 (MQ=255)
                                                                                     cGTGCTGGACGCTCACGTTCTATAGGTCTTGTGATCCCCGATCTGAGAACACCAGCTATACCCGCATCGCTAACTATCTTGAACGCCAGACGCGGCAAc         >  2:299934/1‑99 (MQ=255)
                                                                                        gCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTATCTATCTTGAACGCCAGGCGCGGCAACggg      >  1:765151/1‑100 (MQ=255)
                                                                                        gCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCTTGAACGCCAGGCGCGGCAACggg      >  2:327464/1‑100 (MQ=255)
                                                                                        gCTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCTTGAACGCCAGGCGCGGCAACggg      <  2:1328251/100‑1 (MQ=255)
                                                                                         cTGGACGCTCGCGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCAGCGCTAACTTTCTTGAACGCCAGGCTCGGCAACgtgg     <  2:2484961/100‑1 (MQ=255)
                                                                                         cTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCTTGAACGCCAGGCGCGGCAACgggg     <  1:1461655/100‑1 (MQ=255)
                                                                                         cTGGACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCTTGAACGCCAGGCGCGGCAACgggg     <  1:1224073/100‑1 (MQ=255)
                                                                                          tcgccGCTCACGTTCTTTTTGTCTTATGATTCCCGATCTGGAGACCACCAGCTATACCCGCATCGCTATCTATCTTGAATGCTATTCGCGGCAACGGGGt    <  2:246044/96‑1 (MQ=255)
                                                                                           ggACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCTTGAACGCCAGGCGCGGCAACGGGGtt   <  1:405608/100‑1 (MQ=255)
                                                                                            gACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCTTGAACGCCAGGCGCGGCAACGGGGTTa  <  1:196802/100‑1 (MQ=255)
                                                                                            gACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCTTGAACGCCAGGCGCGGCAACGGGGTTa  <  1:248221/100‑1 (MQ=255)
                                                                                            gACGCTCACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCTTGAACGCCAGGCGCGGCAACGGGGTTa  >  1:2422752/1‑100 (MQ=255)
                                                                                                 |                                                                                              
CCGTGTGAGCGACAAAACCGTTGAAAAAGTCATGGCTGTGGTGCGTGAGCACAATTACCACCCGAACGCCGTGGCAGCTGGGCTTCGTGCTGGACGCACACGTTCTATTGGTCTTGTGATCCCCGATCTGGAGAACACCAGCTATACCCGCATCGCTAACTATCTTGAACGCCAGGCGCGGCAACGGGGTTA  >  NC_000913/88111‑88302

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: