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% 473.9 / ‑4.5 148T318A (ACT→GCT) aroH3‑deoxy‑D‑arabinoheptulosonate‑7‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (52/95);  total (53/95)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 153 total aligned reads displayed.

ATCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGC  >  NC_000913/1787320‑1787499
                                                                                          |                                                                                         
aTCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCcgctcgct                                                                                         >  2:5174900/1‑93 (MQ=255)
aTCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCcgctcgct                                                                                         <  2:2403102/93‑1 (MQ=255)
 tCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCtt                                                                                        >  2:5479475/1‑93 (MQ=255)
  ccGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTa                                                                                       <  1:4578306/93‑1 (MQ=255)
   cGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTAc                                                                                      >  2:6065326/1‑93 (MQ=255)
    gCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACg                                                                                     <  2:1561242/93‑1 (MQ=255)
    gCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACg                                                                                     >  2:151769/1‑93 (MQ=255)
    gCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACg                                                                                     >  2:5898857/1‑93 (MQ=255)
    gCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGAAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACg                                                                                     <  2:6297682/93‑1 (MQ=255)
       aTGGCTTTACGGCTATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACTCAAAAAATCGTCGGCAGCCAGCCGCTCGCTTACGGTc                                                                                  <  1:6477763/93‑1 (MQ=255)
         ggCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCaa                                                                                <  2:1074360/93‑1 (MQ=255)
          gCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAAt                                                                               >  2:5355220/1‑93 (MQ=255)
          gCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAAt                                                                               <  1:2938127/93‑1 (MQ=255)
           ctctACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATc                                                                              <  2:3142661/93‑1 (MQ=255)
            tctACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATcc                                                                             >  2:351929/1‑93 (MQ=255)
            tctACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATcc                                                                             >  1:2324747/1‑93 (MQ=255)
             ctACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCa                                                                            >  1:6123883/1‑93 (MQ=255)
              tACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCAt                                                                           <  1:6432947/93‑1 (MQ=255)
                cGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTa                                                                         <  2:2101253/93‑1 (MQ=255)
                cGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTa                                                                         <  2:6169800/93‑1 (MQ=255)
                  gCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTAcc                                                                       <  2:4247025/93‑1 (MQ=255)
                  gCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTAcc                                                                       >  1:561546/1‑93 (MQ=255)
                     aTTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGAc                                                                    >  2:5050429/1‑93 (MQ=255)
                     aTTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGAc                                                                    >  2:4548230/1‑93 (MQ=255)
                     aTTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGAc                                                                    >  2:1646559/1‑93 (MQ=255)
                      ttGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGAcc                                                                   >  2:2327990/1‑93 (MQ=255)
                      ttGCTGGAATTATGGCCGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGAcc                                                                   <  2:5437818/93‑1 (MQ=255)
                        gCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCg                                                                 <  2:5827032/93‑1 (MQ=255)
                         cTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCgt                                                                >  2:1363202/1‑93 (MQ=255)
                          tGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCgtg                                                               >  1:5215955/1‑93 (MQ=255)
                            gAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTc                                                             <  1:3356428/93‑1 (MQ=255)
                              aTTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctg                                                           >  2:271882/1‑93 (MQ=255)
                              aTTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctg                                                           >  1:2519813/1‑93 (MQ=255)
                               ttATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgg                                                          <  1:6150589/93‑1 (MQ=255)
                               ttATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgg                                                          >  1:1121245/1‑93 (MQ=255)
                               ttATGGCGAAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgg                                                          <  1:2739309/93‑1 (MQ=255)
                                 aTGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggc                                                        >  1:4888762/1‑93 (MQ=255)
                                 aTGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggc                                                        <  1:4005297/93‑1 (MQ=255)
                                 aTGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggc                                                        <  1:4672765/93‑1 (MQ=255)
                                 aTGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggc                                                        <  1:1761518/93‑1 (MQ=255)
                                 aTGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggc                                                        >  2:5892554/1‑93 (MQ=255)
                                 aTGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggc                                                        >  1:3001107/1‑93 (MQ=255)
                                 aTGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggc                                                        <  2:347154/93‑1 (MQ=255)
                                 aTGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggc                                                        >  1:3523907/1‑93 (MQ=255)
                                  tGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggct                                                       >  1:3467293/1‑93 (MQ=255)
                                   ggCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggctg                                                      <  1:4694445/93‑1 (MQ=255)
                                     cGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggctggg                                                    <  2:5872009/93‑1 (MQ=255)
                                      ggAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGgga                                                   <  2:4596137/93‑1 (MQ=255)
                                      ggAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGgga                                                   >  1:1652504/1‑93 (MQ=255)
                                       gAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGCCCCGTGTCTGGGCTGggag                                                  <  1:1970866/93‑1 (MQ=255)
                                         aaGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGggagga                                                >  1:2957230/1‑93 (MQ=255)
                                         aaGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGggagga                                                <  2:5739347/93‑1 (MQ=255)
                                           gTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATa                                              >  1:1424432/1‑93 (MQ=255)
                                            tttCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATAc                                             >  2:2273136/1‑93 (MQ=255)
                                             ttCTTGCGCGAGGAAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATAcc                                            <  1:69278/93‑1 (MQ=255)
                                               ccTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGa                                          >  1:2397896/1‑93 (MQ=255)
                                               ccTGCGCGAAGCAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGa                                          <  1:39414/93‑1 (MQ=255)
                                                cTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGaa                                         <  2:3094334/93‑1 (MQ=255)
                                                cTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGaa                                         <  2:2773044/93‑1 (MQ=255)
                                                 tGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAAc                                        <  2:4161109/93‑1 (MQ=255)
                                                 tGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAAc                                        <  1:1891416/93‑1 (MQ=255)
                                                    gcgAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGcc                                     >  2:5876279/1‑93 (MQ=255)
                                                      gAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTg                                   >  1:3589068/1‑93 (MQ=255)
                                                      gAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTg                                   >  1:199569/1‑93 (MQ=255)
                                                        cggAACGCAAGAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCGATCCATTGCCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGt                                 <  2:6733760/92‑1 (MQ=255)
                                                        aGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGt                                 <  2:1709771/93‑1 (MQ=255)
                                                        aGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGt                                 >  1:438100/1‑93 (MQ=255)
                                                        aGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGt                                 >  1:483564/1‑93 (MQ=255)
                                                         cgAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTc                                >  1:3322850/2‑93 (MQ=255)
                                                          gAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCg                               <  1:3836064/93‑1 (MQ=255)
                                                           aaCGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGa                              <  1:3795575/93‑1 (MQ=255)
                                                           aaCGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGa                              <  1:135986/93‑1 (MQ=255)
                                                               cAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGaaaaa                          >  2:1361905/1‑93 (MQ=255)
                                                               cAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGaaaaa                          <  1:3186324/93‑1 (MQ=255)
                                                               aaaaaaaTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGaaaaa                          <  2:5239774/92‑1 (MQ=255)
                                                                aaaaaaTCGTCGGCAGTCAGCCGCTCGCTTACGGTCACTCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAAc                         <  2:5021912/93‑1 (MQ=255)
                                                                aaaaaaTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAAc                         <  2:4303253/93‑1 (MQ=255)
                                                                    agtcgtcgGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGc                     <  2:4653786/91‑1 (MQ=255)
                                                                     aTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGcc                    <  1:5890032/93‑1 (MQ=255)
                                                                     aTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGcc                    <  1:3605677/93‑1 (MQ=255)
                                                                      tcgtcgGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCct                   <  2:2060481/93‑1 (MQ=255)
                                                                       cgtcgGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCctc                  <  1:6359996/93‑1 (MQ=255)
                                                                       cgtcgGCAGTCAGCCCCTCGCTTAGGGTCATTCCATTACCGACCCGTGTCGGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCctc                  <  1:23785/93‑1 (MQ=255)
                                                                          cgGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGc               <  2:604553/93‑1 (MQ=255)
                                                                             cAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGt            <  1:4919912/93‑1 (MQ=255)
                                                                             cAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGt            <  1:4842434/93‑1 (MQ=255)
                                                                              aGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAATCTCGCCTCTGCGGTa           <  1:4218924/93‑1 (MQ=255)
                                                                               gTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAg          >  1:871200/1‑93 (MQ=255)
                                                                               gTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAg          >  1:4578470/1‑93 (MQ=255)
                                                                               gTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAg          >  1:3189179/1‑93 (MQ=255)
                                                                                tCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGa         <  1:957685/93‑1 (MQ=255)
                                                                                tCAGCCGCTCGCTTACGATCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGa         <  1:5947945/93‑1 (MQ=255)
                                                                                 cAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGAt        <  2:2449515/93‑1 (MQ=255)
                                                                                 cAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGAt        <  1:271399/93‑1 (MQ=255)
                                                                                    ccGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATAcc     <  1:4400901/93‑1 (MQ=255)
                                                                                    ccGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATAcc     <  2:2559877/93‑1 (MQ=255)
                                                                                     cgctcgctTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATAccc    >  1:2765503/1‑93 (MQ=255)
                                                                                       cTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGc  >  1:2187514/1‑93 (MQ=255)
                                                                                       cTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGc  >  1:1889288/1‑93 (MQ=255)
                                                                                       cTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGc  <  1:1175001/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGC  >  NC_000913/1787320‑1787499

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