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.7% 2259.7 / ‑2.5 597T578T (ACC→ACA) gspDgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (306/289);  total (307/290)
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.95e-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 604 total aligned reads displayed.

TACTATTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACCGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGT  >  NC_000913/3457963‑3458256
                                                                                                                                                   |                                                                                                                                                  
tACTATTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGa                                                                                                                                                  >  2:148512/1‑150 (MQ=255)
tACTATTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGa                                                                                                                                                  >  2:2864219/1‑150 (MQ=255)
     ttCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAgcgcg                                                                                                                                             <  1:765723/150‑1 (MQ=255)
     ttCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAgcgcg                                                                                                                                             <  2:7284286/150‑1 (MQ=255)
         aaaCGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTaa                                                                                                                                         <  1:714911/150‑1 (MQ=255)
            cGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACg                                                                                                                                      >  2:4414199/1‑150 (MQ=255)
              ccGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCa                                                                                                                                    <  1:4465552/150‑1 (MQ=255)
                     ggTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGAt                                                                                                                             <  2:7211520/150‑1 (MQ=255)
                      gTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATg                                                                                                                            <  1:8803694/150‑1 (MQ=255)
                      gTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATg                                                                                                                            <  1:2178235/150‑1 (MQ=255)
                      gTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATg                                                                                                                            >  2:3022917/1‑150 (MQ=255)
                       tCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATgg                                                                                                                           >  2:6057967/1‑150 (MQ=255)
                       tCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATgg                                                                                                                           >  2:3448787/1‑150 (MQ=255)
                          aaaCCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTAt                                                                                                                        <  1:6298338/150‑1 (MQ=255)
                           aaCCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTAtt                                                                                                                       >  1:260/1‑150 (MQ=255)
                            aCCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTAttt                                                                                                                      >  2:8236580/1‑150 (MQ=255)
                             ccGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCATCCTGATGGTATTTa                                                                                                                     >  2:9125736/1‑150 (MQ=255)
                              cGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTAt                                                                                                                    <  1:8060072/150‑1 (MQ=255)
                              cGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCACCACAGAGCGCGCTAAACGCAACCTGATGGTATTTAt                                                                                                                    >  2:8635111/1‑150 (MQ=255)
                                gTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTAtcc                                                                                                                  >  1:1701731/1‑150 (MQ=255)
                                    aaCGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTAtccgtcc                                                                                                              <  1:376673/150‑1 (MQ=255)
                                      cGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCcga                                                                                                            <  2:8337811/150‑1 (MQ=255)
                                       ggtggtCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCcgac                                                                                                           <  2:9310/150‑1 (MQ=255)
                                       ggtggtCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCcgac                                                                                                           <  1:9158073/150‑1 (MQ=255)
                                         tggtCCTGTGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCcgacga                                                                                                         <  2:6498187/150‑1 (MQ=255)
                                           gtCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGAtt                                                                                                       <  1:1776569/150‑1 (MQ=255)
                                           gtCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGAtt                                                                                                       >  1:6871653/1‑150 (MQ=255)
                                            tCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTa                                                                                                      >  2:254114/1‑150 (MQ=255)
                                              cTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATc                                                                                                    <  2:551911/150‑1 (MQ=255)
                                               tGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATcc                                                                                                   <  2:7110706/150‑1 (MQ=255)
                                                 ggCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGt                                                                                                 <  1:2349228/150‑1 (MQ=255)
                                                 ggCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGt                                                                                                 <  1:5542654/150‑1 (MQ=255)
                                                            ggATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGt                                                                                      <  2:7529069/150‑1 (MQ=255)
                                                                 atTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATc                                                                                 <  2:4627138/150‑1 (MQ=255)
                                                                  tttttCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCg                                                                                >  1:1433905/1‑150 (MQ=255)
                                                                   ttttCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGc                                                                               <  1:3605636/150‑1 (MQ=255)
                                                                      tCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCa                                                                            <  2:2338900/150‑1 (MQ=255)
                                                                      tCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCa                                                                            <  2:2060062/150‑1 (MQ=255)
                                                                        tAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACt                                                                          >  1:6624355/1‑150 (MQ=255)
                                                                        tAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACt                                                                          >  2:1350775/1‑150 (MQ=255)
                                                                        tAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACt                                                                          >  1:968003/1‑150 (MQ=255)
                                                                        tAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACt                                                                          <  1:3959618/150‑1 (MQ=255)
                                                                           agagCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTc                                                                       <  1:7864158/150‑1 (MQ=255)
                                                                           agagCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTc                                                                       <  1:1892702/150‑1 (MQ=255)
                                                                               cAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCaaaa                                                                   <  2:357648/150‑1 (MQ=255)
                                                                                 aGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAaga                                                                 >  1:3004130/1‑150 (MQ=255)
                                                                                     tCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGaaa                                                             >  1:2048262/1‑150 (MQ=255)
                                                                                     tCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGaaa                                                             <  1:5312809/150‑1 (MQ=255)
                                                                                      cAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAAt                                                            >  2:9101023/1‑150 (MQ=255)
                                                                                       aaaGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATa                                                           >  2:3234864/1‑150 (MQ=255)
                                                                                         aGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATaca                                                         >  1:852411/1‑150 (MQ=255)
                                                                                         aGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATaca                                                         <  1:2811332/150‑1 (MQ=255)
                                                                                           gTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACAcc                                                       >  1:3162797/1‑150 (MQ=255)
                                                                                                     ttGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTc                                             <  2:3653673/150‑1 (MQ=255)
                                                                                                       ggCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAg                                           <  2:3940063/150‑1 (MQ=255)
                                                                                                       ggCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACACAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAg                                           <  1:6058449/150‑1 (MQ=255)
                                                                                                         cGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGa                                         <  2:2681428/150‑1 (MQ=255)
                                                                                                          gATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAg                                        <  2:6871889/150‑1 (MQ=255)
                                                                                                             atTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGcaa                                     <  2:594727/150‑1 (MQ=255)
                                                                                                             atTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGcaa                                     >  2:1804374/1‑150 (MQ=255)
                                                                                                               tCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGcaaca                                   >  1:8678303/1‑150 (MQ=255)
                                                                                                                  tttAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAacaaca                                >  2:5280734/1‑150 (MQ=255)
                                                                                                                  tttAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAacaaca                                >  1:8192428/1‑150 (MQ=255)
                                                                                                                  tttAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAacaaca                                >  1:6790247/1‑150 (MQ=255)
                                                                                                                   ttAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAg                               >  1:598081/1‑150 (MQ=255)
                                                                                                                    tAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGc                              <  2:8914451/150‑1 (MQ=255)
                                                                                                                     aGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAATAACAGCg                             <  1:420418/150‑1 (MQ=255)
                                                                                                                       tGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGa                           <  1:6919385/150‑1 (MQ=255)
                                                                                                                       tGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGa                           >  2:4085248/1‑150 (MQ=255)
                                                                                                                       tGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGa                           <  1:5909804/150‑1 (MQ=255)
                                                                                                                         gggCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGGTGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCGACACCTGCGGATc                         >  1:5382976/1‑150 (MQ=255)
                                                                                                                         gggCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATc                         <  1:2512679/150‑1 (MQ=255)
                                                                                                                         gggCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATc                         >  2:6807490/1‑150 (MQ=255)
                                                                                                                         gggCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATc                         >  2:7910336/1‑150 (MQ=255)
                                                                                                                          ggCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCg                        >  1:6479711/1‑150 (MQ=255)
                                                                                                                           gtaaCTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATACGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGa                       >  1:5224137/3‑150 (MQ=255)
                                                                                                                           gCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGa                       >  2:3102152/1‑150 (MQ=255)
                                                                                                                           gCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGa                       >  1:2459198/1‑150 (MQ=255)
                                                                                                                           gCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGa                       >  2:1116057/1‑150 (MQ=255)
                                                                                                                            cAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGAc                      <  2:5205889/150‑1 (MQ=255)
                                                                                                                               ctctTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACggg                   >  1:3410244/1‑150 (MQ=255)
                                                                                                                                tctTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGa                  <  1:485093/150‑1 (MQ=255)
                                                                                                                                 ctTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGac                 >  1:5774079/1‑149 (MQ=255)
                                                                                                                                 ctTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGaa                 >  2:7696539/1‑150 (MQ=255)
                                                                                                                                  ttCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGaaa                <  2:5809772/150‑1 (MQ=255)
                                                                                                                                  ttCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGaaa                <  1:4099567/150‑1 (MQ=255)
                                                                                                                                   tCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGaact               >  1:9004645/1‑148 (MQ=255)
                                                                                                                                   tCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAAt               <  1:8374315/150‑1 (MQ=255)
                                                                                                                                   tCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAAt               >  1:4381288/1‑150 (MQ=255)
                                                                                                                                    ccGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTTTTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATc              <  2:8488806/150‑1 (MQ=255)
                                                                                                                                    ccGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATc              >  2:4847431/1‑150 (MQ=255)
                                                                                                                                     cGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCTGATCGACGGGAAATCa             <  1:5642591/150‑1 (MQ=255)
                                                                                                                                      gCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCaa            <  2:223031/150‑1 (MQ=255)
                                                                                                                                       caaaaCCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCAGACGATTATCCGTGACGATGATGTTTATCGCTCACTGACAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACTGGAAATaaaa           >  2:2820148/5‑150 (MQ=255)
                                                                                                                                        tataCCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCaaaa          >  2:5756441/1‑150 (MQ=255)
                                                                                                                                             cTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCt     <  2:6360845/150‑1 (MQ=255)
                                                                                                                                             cTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCt     >  1:4755938/1‑150 (MQ=255)
                                                                                                                                             cTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCt     >  1:10093/1‑150 (MQ=255)
                                                                                                                                             cTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCt     <  2:428332/150‑1 (MQ=255)
                                                                                                                                                cacaGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGt  >  1:1066507/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TACTATTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACCGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGT  >  NC_000913/3457963‑3458256

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