Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,787,410 A→G T318A (ACT→GCT)  aroH → 3‑deoxy‑D‑arabinoheptulosonate‑7‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,787,4100AG99.3% 436.4 / NA 136T318A (ACT→GCT) aroH3‑deoxy‑D‑arabinoheptulosonate‑7‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (51/84);  minor base T (0/1);  total (51/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 148 total aligned reads displayed.

ATCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCG  >  NC_000913/1787320‑1787498
                                                                                          |                                                                                        
aTCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCcgctcgct                                                                                        >  1:5541400/1‑93 (MQ=255)
aTCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCcgctcgct                                                                                        >  1:5861318/1‑93 (MQ=255)
aTCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCcgctcgct                                                                                        <  2:5643368/93‑1 (MQ=255)
aTCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCcgctcgct                                                                                        <  1:4235154/93‑1 (MQ=255)
aTCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCcgctcgct                                                                                        <  1:4274996/93‑1 (MQ=255)
 tCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCtt                                                                                       >  1:1208496/1‑93 (MQ=255)
 tCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCtt                                                                                       <  2:5679952/93‑1 (MQ=255)
  ccGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTa                                                                                      <  1:6227651/93‑1 (MQ=255)
   cGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTAc                                                                                     >  2:5470225/1‑93 (MQ=255)
   cGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTAc                                                                                     >  1:3728199/1‑93 (MQ=255)
    gCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTTGGCAGTCAGCAGCTCGCTTNCg                                                                                    >  1:3353792/1‑93 (MQ=255)
    gCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACg                                                                                    >  2:2884352/1‑93 (MQ=255)
    gCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACg                                                                                    >  2:3699723/1‑93 (MQ=255)
     cAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACgg                                                                                   >  1:5016133/1‑93 (MQ=255)
     cAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACgg                                                                                   >  1:409545/1‑93 (MQ=255)
      aaTGGCTCTACGGCGATTTCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGt                                                                                  <  2:1861568/93‑1 (MQ=255)
      aaTGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGt                                                                                  >  2:2796934/1‑93 (MQ=255)
      aaTGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGt                                                                                  >  1:1430983/1‑93 (MQ=255)
       aTGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTc                                                                                 >  1:4113521/1‑93 (MQ=255)
         ggCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCaa                                                                               >  1:1952079/1‑93 (MQ=255)
          gCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAAt                                                                              >  1:378818/1‑93 (MQ=255)
          gCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAAt                                                                              >  1:1128283/1‑93 (MQ=255)
           ctctACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATc                                                                             <  2:1881590/93‑1 (MQ=255)
            tctACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATcc                                                                            >  2:4757635/1‑93 (MQ=255)
             ctACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCa                                                                           <  1:3539507/93‑1 (MQ=255)
              tACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCAt                                                                          <  1:4300832/93‑1 (MQ=255)
              tACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCAt                                                                          <  2:3168978/93‑1 (MQ=255)
               aCGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCAtt                                                                         >  2:3411245/1‑93 (MQ=255)
               aCGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCAtt                                                                         >  2:257424/1‑93 (MQ=255)
               aCGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCAtt                                                                         <  1:1808303/93‑1 (MQ=255)
                gggCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTa                                                                        <  2:1687046/92‑1 (MQ=255)
                cGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTa                                                                        <  1:1277077/93‑1 (MQ=255)
                 ggCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTAc                                                                       <  1:4907476/93‑1 (MQ=255)
                 ggCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTAc                                                                       >  1:950062/1‑93 (MQ=255)
                 ggCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTAc                                                                       <  2:1842352/93‑1 (MQ=255)
                 ggCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTAc                                                                       <  2:6034540/93‑1 (MQ=255)
                 ggCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTAc                                                                       <  2:44106/93‑1 (MQ=255)
                 ggCGATTGCTGAAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTAc                                                                       <  2:5439246/93‑1 (MQ=255)
                  gCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTAcc                                                                      >  1:3332417/1‑93 (MQ=255)
                     aTTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGAc                                                                   >  2:911249/1‑93 (MQ=255)
                     aTTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGAc                                                                   >  1:2295798/1‑93 (MQ=255)
                     aTTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGAc                                                                   >  1:3906302/1‑93 (MQ=255)
                      ttGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGAcc                                                                  >  1:5003745/1‑93 (MQ=255)
                      ttGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGAcc                                                                  <  2:6078056/93‑1 (MQ=255)
                          tGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCgtg                                                              >  2:3648187/1‑93 (MQ=255)
                          tGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCgtg                                                              >  2:6235043/1‑93 (MQ=255)
                          tGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCgtg                                                              <  2:2381851/93‑1 (MQ=255)
                           ggAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCgtgt                                                             >  1:4148672/1‑93 (MQ=255)
                             aaTTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTct                                                           <  2:3427761/93‑1 (MQ=255)
                              aTTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctg                                                          >  2:3283100/1‑93 (MQ=255)
                              aTTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctg                                                          >  2:329131/1‑93 (MQ=255)
                              aTTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctg                                                          <  2:2531586/93‑1 (MQ=255)
                               ttATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgg                                                         <  1:3158774/93‑1 (MQ=255)
                               ttATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgg                                                         <  2:4187927/93‑1 (MQ=255)
                                   ggCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggctg                                                     >  1:939032/1‑93 (MQ=255)
                                    gCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggctgg                                                    <  2:4888643/93‑1 (MQ=255)
                                     cGGAAAGTTTCCTACGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggctggg                                                   <  1:3602188/93‑1 (MQ=255)
                                      ggAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGgga                                                  >  2:160761/1‑93 (MQ=255)
                                       gAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGggag                                                 <  2:2934589/93‑1 (MQ=255)
                                       gAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGggag                                                 >  1:1038918/1‑93 (MQ=255)
                                         aaGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGggagga                                               <  2:1765507/93‑1 (MQ=255)
                                          aGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGAt                                              >  1:248250/1‑93 (MQ=255)
                                            tttCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATAc                                            <  2:895061/93‑1 (MQ=255)
                                             ttCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATAcc                                           <  1:321496/93‑1 (MQ=255)
                                             tcccTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATAcc                                           <  2:4228228/91‑1 (MQ=255)
                                              tCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCTCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCg                                          <  1:2602456/93‑1 (MQ=255)
                                              tCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCg                                          >  2:2171331/1‑93 (MQ=255)
                                              tCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCg                                          <  2:1396738/93‑1 (MQ=255)
                                              tCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCg                                          >  1:449404/1‑93 (MQ=255)
                                                cTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGaa                                        <  2:5950091/93‑1 (MQ=255)
                                                  gcgcGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACg                                      <  1:3586047/93‑1 (MQ=255)
                                                       aaGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTgg                                 >  1:4609714/1‑93 (MQ=255)
                                                        aGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGt                                <  2:456675/93‑1 (MQ=255)
                                                        aGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGt                                <  2:2743442/93‑1 (MQ=255)
                                                        aGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGt                                <  2:4612362/93‑1 (MQ=255)
                                                           aaCGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGa                             <  1:1778710/93‑1 (MQ=255)
                                                                aaaaaaTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAAc                        <  2:5209246/93‑1 (MQ=255)
                                                                aaaaaaTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAAc                        >  1:1036249/1‑93 (MQ=255)
                                                                 aaaaaTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACt                       <  1:5981398/93‑1 (MQ=255)
                                                                 aaaaaTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACt                       <  2:6275264/93‑1 (MQ=255)
                                                                  aaaaTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTc                      <  2:4451288/93‑1 (MQ=255)
                                                                  aaaaTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTc                      <  1:3877238/93‑1 (MQ=255)
                                                                  aaaaTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTc                      <  1:3726064/93‑1 (MQ=255)
                                                                  aaaaGCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTc                      <  1:5421615/93‑1 (MQ=255)
                                                                    aaTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGc                    <  1:3574862/93‑1 (MQ=255)
                                                                     aTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGcc                   <  1:1682807/93‑1 (MQ=255)
                                                                       cgtcgGCAGTCAGCCGCTCGCTTTCGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCctc                 >  2:5070875/1‑93 (MQ=255)
                                                                       cgtcgGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCctc                 <  1:6229954/93‑1 (MQ=255)
                                                                        gtcgGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCctct                <  1:2543661/93‑1 (MQ=255)
                                                                           ggCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCg             <  1:2675034/93‑1 (MQ=255)
                                                                             cAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGt           <  2:3131182/93‑1 (MQ=255)
                                                                             cAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGt           <  2:4799297/93‑1 (MQ=255)
                                                                              cgTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTa          <  2:1016983/92‑1 (MQ=255)
                                                                               gTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAg         >  1:4456362/1‑93 (MQ=255)
                                                                                tCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGa        >  2:1078710/1‑93 (MQ=255)
                                                                                 cAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGAt       <  1:2770170/93‑1 (MQ=255)
                                                                                    ccGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATAcc    <  2:5981735/93‑1 (MQ=255)
                                                                                    ccGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATAcc    <  2:5379010/93‑1 (MQ=255)
                                                                                     cgctcgctTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATAccc   <  1:1014878/93‑1 (MQ=255)
                                                                                      gCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCg  <  2:4028971/93‑1 (MQ=255)
                                                                                          |                                                                                        
ATCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCG  >  NC_000913/1787320‑1787498

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑