Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,458,110 C→A 100% T578T (ACC→ACA)  gspD → general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,458,1100CA99.2% 421.7 / ‑5.4 123T578T (ACC→ACA) gspDgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (61/60);  total (62/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 127 total aligned reads displayed.

ACTATTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACCGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGT  >  NC_000913/3457964‑3458256
                                                                                                                                                  |                                                                                                                                                  
aCTATTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACagag                                                                                                                                                 <  1:1494563/150‑1 (MQ=255)
 cTATTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAgc                                                                                                                                                <  1:1187905/150‑1 (MQ=255)
 atATTCAAACCGCCGTGCTGGTCAAAACCGGTGTATCGTTTGTCCTGCGCGGATTGCTTGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACCGAgc                                                                                                                                                <  2:729326/149‑1 (MQ=255)
   aTTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTTAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCTACTCTTCCGCTATACCTCCACAGAgcgc                                                                                                                                              >  1:265143/1‑150 (MQ=255)
     tCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAgcgcgc                                                                                                                                            >  2:2101156/1‑150 (MQ=255)
      cAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCt                                                                                                                                           >  2:1281718/1‑150 (MQ=255)
      cAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCt                                                                                                                                           >  2:675887/1‑150 (MQ=255)
      cAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTATAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCt                                                                                                                                           >  2:1500484/1‑150 (MQ=255)
        aaaCGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAAATCTTCCGCTATACCTACACAGAGCGCGCTaa                                                                                                                                         >  2:1369944/1‑150 (MQ=255)
           cGCCGTGCTTGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACg                                                                                                                                      <  1:1836907/150‑1 (MQ=255)
              cGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATGTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCaa                                                                                                                                   <  1:1347981/150‑1 (MQ=255)
                tGCTGTGCAAAACCTGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGATCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAAcc                                                                                                                                 <  1:593271/150‑1 (MQ=255)
                     gTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCTCTTAACGCAACCTGATg                                                                                                                            >  2:1148988/1‑150 (MQ=255)
                     gTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATg                                                                                                                            >  1:464867/1‑150 (MQ=255)
                      tCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATgg                                                                                                                           <  2:762778/150‑1 (MQ=255)
                      tCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATgg                                                                                                                           >  2:1677884/1‑150 (MQ=255)
                       cAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTGAGTGGGGCATCTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGt                                                                                                                          >  2:916130/1‑150 (MQ=255)
                           aCCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTAttt                                                                                                                      >  2:1766483/1‑150 (MQ=255)
                            ccGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTa                                                                                                                     >  1:580558/1‑150 (MQ=255)
                            ccGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTa                                                                                                                     >  2:1164739/1‑150 (MQ=255)
                              ggTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTAtc                                                                                                                   <  2:531345/150‑1 (MQ=255)
                              ggTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTAtc                                                                                                                   >  1:1253958/1‑150 (MQ=255)
                               gTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTAtcc                                                                                                                  >  1:698093/1‑150 (MQ=255)
                                  aaaCGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCTAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTAtccgtc                                                                                                               >  2:324398/1‑150 (MQ=255)
                                   aaCGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTAtccgtcc                                                                                                              >  1:1853842/1‑150 (MQ=255)
                                         ggtCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGAt                                                                                                        >  2:1952694/1‑150 (MQ=255)
                                             cTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATc                                                                                                    >  1:1687185/1‑150 (MQ=255)
                                               gggCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCg                                                                                                  >  2:1171662/1‑150 (MQ=255)
                                                ggCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGt                                                                                                 <  1:131700/150‑1 (MQ=255)
                                                 gCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTg                                                                                                <  2:845588/150‑1 (MQ=255)
                                                 gCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTg                                                                                                <  1:748891/150‑1 (MQ=255)
                                                   ggATTGCTGGCTGATTTTTCTAAAGAGCAAGTGTCAAAGGGTCCTCTGCTTGGCGTTATTCCTTTAGTGGGGCAACTCTTCCGCTAGACCTCCCCAGAGCGCGCTCACCGCAACCTGATGGTTTTTATCCGTCCGACGATTATCCGGGAc                                                                                              <  1:1511638/150‑1 (MQ=255)
                                                         cTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACgatgat                                                                                        <  1:1597056/150‑1 (MQ=255)
                                                             ctgatTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGttt                                                                                    <  2:890256/149‑1 (MQ=255)
                                                             atgatTTTTCTAAAGAGCAAGTGTCAAAGGGTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGttt                                                                                    <  2:74179/150‑1 (MQ=255)
                                                               gatTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATTTTTAt                                                                                  >  2:1862891/1‑150 (MQ=255)
                                                               gatTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTAt                                                                                  <  1:1045265/150‑1 (MQ=255)
                                                                atTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATc                                                                                 >  1:411009/1‑150 (MQ=255)
                                                                 tttttCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCg                                                                                <  1:2025747/150‑1 (MQ=255)
                                                                 tttttCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAAAGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCg                                                                                >  1:1018562/1‑150 (MQ=255)
                                                                  ttttCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGc                                                                               <  2:424858/150‑1 (MQ=255)
                                                                    ttCTAAAGAGCAAGTGTCAAAGGTGCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTc                                                                             <  1:1454596/150‑1 (MQ=255)
                                                                      cTAAAGAGCAAGTTTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGTGAAACTCTTCCGCTATACCTCCACAGAGAGCGCTAAACGCAACCTGATGGTATTTATCCGTCAGACGATTATCAGTGACGATGATGTTTATCGCTCAc                                                                           >  2:624166/1‑150 (MQ=255)
                                                                        aaaGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATCCCTCAACAGAGCGCGCTCAACGCAACCTGATGGGATTTATCCGTCCGACGATTATCCGTGCCGAGGATGTTTATCGCTCCCTg                                                                         >  2:1002135/1‑150 (MQ=255)
                                                                        aaaGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTg                                                                         >  2:1718799/1‑150 (MQ=255)
                                                                         aaGAGCAAGTGTCAAAGGTTCCTCTTCTTGTCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGt                                                                        <  2:348681/150‑1 (MQ=255)
                                                                             gCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCAGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCaaa                                                                    >  2:2028806/1‑150 (MQ=255)
                                                                             gCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTCCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCaaa                                                                    <  1:429365/150‑1 (MQ=255)
                                                                              cAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCaaaa                                                                   >  1:1587483/1‑150 (MQ=255)
                                                                                 gtgtCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGCTGGTATTTATCCGTCCGCCGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAagag                                                                >  1:1245997/1‑150 (MQ=255)
                                                                                 gtgtCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAagag                                                                >  1:2028103/1‑150 (MQ=255)
                                                                                  tgtCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAgaga                                                               <  1:806777/150‑1 (MQ=255)
                                                                                  tgtCAAAGGTTCCTCAGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGACCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAgaga                                                               <  2:97922/150‑1 (MQ=255)
                                                                                   gtCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGaa                                                              >  1:295193/1‑150 (MQ=255)
                                                                                   gtCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGGGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGaa                                                              <  1:1259118/150‑1 (MQ=255)
                                                                                    tCAAAGGTTCCTCTGCTTGTCGATAATCCTTTATTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTCTCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGaaa                                                             <  2:1570945/150‑1 (MQ=255)
                                                                                     cAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAAt                                                            >  2:266671/1‑150 (MQ=255)
                                                                                     cAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAAt                                                            <  1:1901498/150‑1 (MQ=255)
                                                                                       aaGGGTCCTCTCCTGGGCGATATTCCTTTCCGGGGGCAACTCTTCCGCTATACCTCCACAGCCCGCGCTAACCGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATac                                                          <  2:422103/150‑1 (MQ=255)
                                                                                        aGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATaca                                                         <  1:1689489/150‑1 (MQ=255)
                                                                                           ttgctctGCATGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCCCAGAGCGCGCTAAACGCAACCTGATGGTTTTTATCCGTCCGACGAATATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACAccc                                                      <  2:945886/147‑1 (MQ=255)
                                                                                            tCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAACCGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGCCGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCg                                                     >  1:74224/1‑150 (MQ=255)
                                                                                            tCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCg                                                     <  2:520625/150‑1 (MQ=255)
                                                                                            tCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCATCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCg                                                     <  1:1298076/150‑1 (MQ=255)
                                                                                             ccTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGt                                                    <  1:864192/150‑1 (MQ=255)
                                                                                              ctctGCTTGGCGATATTCCTTTAGTGGGTCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGtt                                                   <  1:1122140/150‑1 (MQ=255)
                                                                                               tctGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAACGAGAAATACACACGTTa                                                  >  1:887833/1‑150 (MQ=255)
                                                                                                  gCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCg                                               >  1:1713976/1‑150 (MQ=255)
                                                                                                    ttGGCGATATTCCTTTAGTGGGGCCACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGtg                                             >  1:376973/1‑149 (MQ=255)
                                                                                                      ggCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAg                                           <  1:1281718/150‑1 (MQ=255)
                                                                                                       gCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAgg                                          >  2:1960690/1‑150 (MQ=255)
                                                                                                       gCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAgg                                          <  2:919226/150‑1 (MQ=255)
                                                                                                       gCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAgg                                          <  1:642421/150‑1 (MQ=255)
                                                                                                        cGATATTCCTTTAGTGGGTCAACTCTTCCGCTATAACTCCACAGAGCGCGCTAAACGCAACCTTATGTTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGa                                         <  2:666539/150‑1 (MQ=255)
                                                                                                        cGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGa                                         <  2:657665/150‑1 (MQ=255)
                                                                                                         gATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAg                                        <  2:1972348/150‑1 (MQ=255)
                                                                                                         gATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAg                                        <  2:1964505/150‑1 (MQ=255)
                                                                                                              tCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGcaaca                                   <  1:213098/150‑1 (MQ=255)
                                                                                                              tCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGcaaca                                   >  2:1820746/1‑150 (MQ=255)
                                                                                                              tCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGcaaca                                   <  2:1602340/150‑1 (MQ=255)
                                                                                                               ccTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGcaacaa                                  >  2:1998486/1‑150 (MQ=255)
                                                                                                                 tttAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGCCGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCaacaacc                                >  2:1412630/1‑149 (MQ=255)
                                                                                                                   tAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACACCAGc                              >  2:1894004/1‑150 (MQ=255)
                                                                                                                   tAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGc                              <  1:28930/150‑1 (MQ=255)
                                                                                                                      tGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGa                           <  1:1952694/150‑1 (MQ=255)
                                                                                                                      aggggCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGACCGACGATTATCCGTGACGATGATGTGTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGa                           >  2:1612896/2‑150 (MQ=255)
                                                                                                                            aaCTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACg                     <  2:411009/150‑1 (MQ=255)
                                                                                                                            aaCTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACg                     <  2:214057/150‑1 (MQ=255)
                                                                                                                                    cGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCa             >  2:765945/1‑150 (MQ=255)
                                                                                                                                        attccTCCACAGAGCTCGCTAAAGGGAACCTGATGGTATTTTTCCGTCCGATGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGTAATACACCCGTTTCCGTCAGGAGCAACAACAGCGTATCGACGGGAAATCAAAAg         <  2:974604/147‑1 (MQ=255)
                                                                                                                                        ataCCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAg         <  2:1737888/150‑1 (MQ=255)
                                                                                                                                         taCCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAgc        >  1:203973/1‑150 (MQ=255)
                                                                                                                                         taCCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAgc        <  1:941703/150‑1 (MQ=255)
                                                                                                                                           ccTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACCCCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGTAATCAAAAgcgc      >  1:1052945/1‑150 (MQ=255)
                                                                                                                                           ccTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGCCGGGAAATCAAAAgcgc      >  1:1704242/1‑150 (MQ=255)
                                                                                                                                           ccTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAgcgc      >  1:54070/1‑150 (MQ=255)
                                                                                                                                            cTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCt     >  1:1675743/1‑150 (MQ=255)
                                                                                                                                              ccACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCTAAAGAGAAATACACCCGTTACCGTCAGGAGCAAAAACAGCGGATCGACGGGAATTCAAAAGCGCTgg   >  2:350951/1‑150 (MQ=255)
                                                                                                                                              ccACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTgg   >  2:435469/1‑150 (MQ=255)
                                                                                                                                               cacaGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGt  <  1:1565152/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ACTATTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACCGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGT  >  NC_000913/3457964‑3458256

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑