Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 83,587 A→C 16.6% intergenic (‑58/+35) leuA ← / ← leuL 2‑isopropylmalate synthase/leu operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091383,5870AC16.6% 246.4 / 46.6 109intergenic (‑58/+35)leuA/leuL2‑isopropylmalate synthase/leu operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (39/50);  new base C (8/10);  total (47/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.85e-01

Only 100 of 113 total aligned reads displayed.

GTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGA  >  NC_000913/83488‑83686
                                                                                                   |                                                                                                    
gTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCa                                                                                                      >  1:909955/1‑100 (MQ=255)
gTCGCGCAATGTGGTATCGAAAATAATGACTTCCTGGCTTATGGCTTGGGTCCTTGTCTCTTTTAGAGCGCCTCTCTTCAGGAATAAAAAAACCCGCGCa                                                                                                      <  2:947152/100‑1 (MQ=255)
  cgcgCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAAt                                                                                                    >  2:1848082/1‑100 (MQ=255)
  cgcgCAATGTGGTATCGAAAATAATGACTTGCTGGCCCATGGTTTGGGTCATTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCa                                                                                                      >  2:854241/1‑98 (MQ=255)
     gcAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGgc                                                                                                 >  1:663076/1‑100 (MQ=255)
        aTGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGcgcg                                                                                              <  1:1406464/100‑1 (MQ=255)
         tgtgGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCgg                                                                                             <  2:893235/100‑1 (MQ=255)
          gtgtTATCGAAAATAATTACTTGCTGGCTCATGGTTTGGGTCCTTGTATATTTTAAAACGAATAAATTAAAAAATAAAAAAACCCGCGCAATGGCGCggg                                                                                            <  2:151362/100‑1 (MQ=255)
          gtgGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCggg                                                                                            <  2:237640/100‑1 (MQ=255)
          gtgGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCggg                                                                                            >  2:373134/1‑100 (MQ=255)
           tgGTATCTAAAATAATGACTTGCTGGCACATGGTTTGGGTCCTTGTCTCTTTTCGAGCGCCGCGCTTCGGGCAT‑AAAAAACCCGCGCCATTGCGCGGG‑tt                                                                                         <  2:2179904/100‑1 (MQ=37)
           tgGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑t                                                                                          <  1:1299849/100‑1 (MQ=255)
           tgGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑t                                                                                          <  1:591622/100‑1 (MQ=255)
           tgGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGCCGCGGG‑t                                                                                          <  1:1253410/100‑1 (MQ=255)
            ggTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑tt                                                                                         <  1:1435665/100‑1 (MQ=255)
            ggTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑tt                                                                                         <  2:891492/100‑1 (MQ=255)
               aTCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑ttttt                                                                                      <  1:514183/100‑1 (MQ=255)
               aTCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑ttttt                                                                                      >  1:268519/1‑100 (MQ=255)
               aTCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑ttttt                                                                                      >  1:2080806/1‑100 (MQ=255)
                tCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑tttttt                                                                                     <  1:291156/100‑1 (MQ=255)
                tCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑tttttt                                                                                     <  1:1974606/100‑1 (MQ=255)
                 cGAAAATAATGACTTGCTGGCTCATTGTTTGGGTCCTTGTTTCTTTTAGAGCGCCTCTCTTCGGGCATAAAAAAAGCCGCGCAATGGCGCGGG‑TTCtttg                                                                                    >  1:495195/1‑100 (MQ=255)
                 cGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTtttg                                                                                    <  1:947158/100‑1 (MQ=255)
                 cGAAAATAATGACTTGCTGACTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCTGGCATAAAAAAACCCGCGCAATGGAGCCGG‑ttttt                                                                                      >  2:1153892/1‑98 (MQ=255)
                  gAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑ttttttta                                                                                   <  1:1768021/100‑3 (MQ=255)
                  gAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCACTGGCGCGGG‑TTTtttgt                                                                                   >  1:22381/1‑100 (MQ=255)
                   aaaaTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTtttgtt                                                                                  <  1:1798226/100‑1 (MQ=255)
                   aaaaTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TCTtttgtt                                                                                  <  1:1252430/100‑1 (MQ=255)
                   aaaaTAATGACGTGCTGGGTCATGGGTTGGGTCCTTGTCTCTTTTAGAGCGCGTCGCTTTTGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTtttgtt                                                                                  <  2:1558376/100‑1 (MQ=255)
                    aaaTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTATAGCGCCTCGCTTCGGGCATAAAACAACCCGCGCAATGGCCCGGG‑TTTtttgttt                                                                                 <  2:487394/100‑1 (MQ=255)
                     aataatGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTtttgtttg                                                                                <  2:1856683/100‑1 (MQ=255)
                      ataatGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑tttttttatgcc                                                                               <  2:312888/100‑7 (MQ=37)
                       taatGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGAc                                                                              <  2:1708201/100‑1 (MQ=255)
                        aatGACTTGCTGGCTCATGGTTTGAGTCCTTGTCCCTTTTAGAGCGCACCGCTTCGTGCATCAACAAACCCGCGCAATGACGCGGG‑TTTCTTGTTCGAc                                                                              <  2:348982/99‑1 (MQ=255)
                         atGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTg                                                                            >  2:1562707/1‑100 (MQ=255)
                          tGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGAc                                                                              <  2:1189803/97‑1 (MQ=255)
                          tGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTg                                                                            <  2:1685228/99‑1 (MQ=255)
                           gACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCg                                                                          <  2:420992/100‑1 (MQ=255)
                                 cTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTgg                                                                    <  2:1311937/100‑1 (MQ=255)
                                  tGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGCGCTgga                                                                   <  1:908063/100‑2 (MQ=255)
                                   ggCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCt                                                                  >  2:2089024/1‑100 (MQ=255)
                                    gCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCtt                                                                 <  2:255099/100‑1 (MQ=255)
                                     cTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTa                                                                >  1:698159/1‑100 (MQ=255)
                                     cTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTa                                                                >  2:1106315/1‑100 (MQ=255)
                                     cTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTa                                                                <  2:1584838/100‑1 (MQ=255)
                                     cTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTa                                                                >  1:886539/1‑100 (MQ=255)
                                     cTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTa                                                                >  2:1410614/1‑100 (MQ=255)
                                       cATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAAt                                                              >  1:580730/1‑100 (MQ=255)
                                       cATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAAt                                                              <  2:1318555/100‑1 (MQ=255)
                                         tGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGc                                                            >  1:589978/1‑100 (MQ=255)
                                         tGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGc                                                            >  1:452669/1‑100 (MQ=255)
                                         tGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGc                                                            <  2:1798271/100‑1 (MQ=255)
                                          ggTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTACTGGCTTAATGCt                                                           >  1:1386226/1‑100 (MQ=255)
                                             ttGGGTCCTTGTCTCCTTTAGAGCGACTTGCTTTCGTCAGAAAAAAGCCCGCGCAATGGCGGGGG‑TTTTTTGTTTGATGGCGTGCTGGCTTAATGCTg                                                          <  2:279135/98‑1 (MQ=255)
                                               gtctCCTGGCCTCTTGTAGAGCGCCTCGCTTCGGGCAGAAAACGACCCGCGCAATGGCGCGGG‑TTTGTCGTTTGACTGCCTGCTGGCTTACTGCTggcgg                                                      <  2:1202061/97‑4 (MQ=255)
                                                  tCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCAGA‑TTTTTTGTTTAACAGCGTGCTGGCTTAATGCTGGATGACg                                                   <  2:943888/100‑1 (MQ=255)
                                                   ccTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGc                                                  <  2:1200042/100‑1 (MQ=255)
                                                     ttGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTc                                                <  1:304729/100‑1 (MQ=255)
                                                       gTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCAc                                              >  1:681004/1‑100 (MQ=255)
                                                       gTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCAc                                              <  1:701174/100‑1 (MQ=255)
                                                         ctctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTc                                            <  1:775684/100‑1 (MQ=255)
                                                          tctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCg                                           <  2:1358418/100‑1 (MQ=255)
                                                          tctTTTAGAGCACCTCTCTTCGGGCATAAAAAGACCCGCGCATTGGCGCAGG‑TTTTTTGTTTGACTGCGTGCTTGCTTAATGCTGGAAACAGCTCACTCg                                           <  2:521429/100‑1 (MQ=255)
                                                           ctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGt                                          >  2:1179137/1‑100 (MQ=255)
                                                              tgagagCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTa                                       >  1:1614744/3‑100 (MQ=255)
                                                               tAGAGCGCGTCGCTTCGGGCATAAAAATTCCCGCGCTATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCt                                        <  2:2053468/98‑1 (MQ=255)
                                                               tAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTAc                                      >  1:2175756/1‑100 (MQ=255)
                                                                 gagCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTAcc                                     >  2:610665/1‑99 (MQ=255)
                                                                   gcgcCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACcgcg                                  >  1:2099340/1‑100 (MQ=255)
                                                                     gcCTCGCTTCGGGCATAAAAAATCCCGCGCAATGTCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCa                                <  2:1830740/100‑1 (MQ=255)
                                                                     gcCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCa                                <  2:2307151/100‑1 (MQ=255)
                                                                     gcCTCGCTTCGGGCATAAAAAAACCCGCGCAACGGCCCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCa                                <  1:269934/100‑1 (MQ=255)
                                                                      ccTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCaa                               <  1:2089041/100‑1 (MQ=255)
                                                                          gCTTCGGGCATAAACAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaagaa                           >  1:715946/1‑100 (MQ=255)
                                                                          gCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaagaa                           >  1:1426624/1‑100 (MQ=255)
                                                                              cgatcttAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATg                       <  2:1903521/94‑1 (MQ=255)
                                                                               gggCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaagaag                          <  2:568311/96‑1 (MQ=255)
                                                                                ggCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCg                     <  2:20120/100‑1 (MQ=255)
                                                                                  aaacaaaaaaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGtt                   >  2:2249510/5‑100 (MQ=38)
                                                                                   aacaaaaaaCCCGCGCCATTGCGCGGGTTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCACAGAAGATGCGt                    >  1:1632433/4‑99 (MQ=37)
                                                                                   aTAAAAAAACCCGCGCTATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGttt                  <  1:1028369/100‑1 (MQ=255)
                                                                                    tAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATgagatcg                 >  1:1903537/1‑94 (MQ=255)
                                                                                    tAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTta                 >  1:2248337/1‑100 (MQ=255)
                                                                                    acaaaaaaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTta                 >  1:545486/3‑100 (MQ=255)
                                                                                     caaaaaaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCCAATAAAATGCGTTtag                >  1:2102888/2‑100 (MQ=38)
                                                                                     caaaaaaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtag                >  2:1516115/2‑100 (MQ=255)
                                                                                     caaaaaaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtag                >  2:2081841/2‑100 (MQ=255)
                                                                                     aaaaaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtag                >  1:1787660/1‑100 (MQ=255)
                                                                                     aaaaaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtag                >  2:1096987/1‑100 (MQ=255)
                                                                                      aaaaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagt               <  1:400527/100‑1 (MQ=255)
                                                                                        aaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtag             <  1:922204/100‑1 (MQ=255)
                                                                                        aaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtag             <  1:1487778/100‑1 (MQ=255)
                                                                                        aaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtag             <  2:874938/100‑1 (MQ=255)
                                                                                          aaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagta           <  1:306438/100‑1 (MQ=255)
                                                                                            cccGCGCCATTGCGCGGG‑TTTTTTGTTTGACTCCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtagt         >  1:758549/1‑100 (MQ=255)
                                                                                            cccGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtagt         >  1:1676187/1‑100 (MQ=255)
                                                                                            cccGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtagt         >  1:2384401/1‑100 (MQ=255)
                                                                                            cccGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGACTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtag          <  2:750955/99‑1 (MQ=255)
                                                                                                cgcAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGAc     >  1:456896/1‑100 (MQ=255)
                                                                                                  cAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCg   >  1:932696/1‑100 (MQ=255)
                                                                                                   |                                                                                                    
GTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGA  >  NC_000913/83488‑83686

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: ‑