Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,411,263 A→C 7.6% A30A (GCT→GCG yfbS ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,411,2630AC 7.6% 2889.9 / 43.6 994A30A (GCT→GCGyfbSputative transporter
Reads supporting (aligned to +/- strand):  ref base A (593/297);  new base C (2/101);  total (596/398)
Fisher's exact test for biased strand distribution p-value = 3.26e-41
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1057 total aligned reads displayed.

CCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGA  >  NC_000913/2411114‑2411411
                                                                                                                                                      |                                                                                                                                                    
  cATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGc                                                                                                                                                    >  1:6224212/1‑150 (MQ=255)
  cATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGc                                                                                                                                                    <  2:7726126/150‑1 (MQ=255)
   aTTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGCCACTAACAAACAAAGCg                                                                                                                                                   >  1:713943/1‑150 (MQ=255)
    ttACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGa                                                                                                                                                  <  2:3933416/150‑1 (MQ=255)
       cGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCg                                                                                                                                               >  1:2834467/1‑150 (MQ=255)
         gTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCa                                                                                                                                             >  2:2375109/1‑150 (MQ=255)
             cAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCa                                                                                                                                         >  1:10157488/1‑150 (MQ=255)
              aaCACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCAt                                                                                                                                        >  1:6557701/1‑150 (MQ=255)
               acacCTGTACGGACAAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGAAAACGTTAGGATCATAAAAGCCGTAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGAAAACGCGACAATAACAAACAAAGCGACCGCATCCATa                                                                                                                                       <  2:8266615/150‑1 (MQ=255)
               acacCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGTACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATa                                                                                                                                       <  2:2742156/150‑1 (MQ=255)
               acacCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATa                                                                                                                                       >  2:1538145/1‑150 (MQ=255)
               acacCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATa                                                                                                                                       >  2:655451/1‑150 (MQ=255)
                ccccGGTACGGACCAAACCAGCGCCAATAATAAACAAGGCGGCAATCAGGACAAGGTTAGGATCAGAAAAGCCGGAAAATAGTTCTGGGACTGTCTGCTTTCCGATTACTGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATAc                                                                                                                                      <  2:7915855/148‑1 (MQ=255)
                   cGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCa                                                                                                                                   >  1:5387222/1‑150 (MQ=255)
                   cGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCa                                                                                                                                   >  2:1967029/1‑150 (MQ=255)
                   cGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCa                                                                                                                                   >  2:6131363/1‑150 (MQ=255)
                    ggTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCAc                                                                                                                                  >  1:5867831/1‑150 (MQ=255)
                      tACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCActc                                                                                                                                >  1:1335717/1‑150 (MQ=255)
                         gtaCCAAACATTCGCAAATATTAAACATGGCGTCATTAGGGAGATCTTTAGTTTCAAAATAGCGGTAAAATATTTCTGGGAGTTTCAGCTTCGGTTTTACTGCAACCTCTACAATAACAAACAACGCGACCGCATCCATACGCACTCTGc                                                                                                                             <  2:4127730/148‑1 (MQ=255)
                          gACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGcc                                                                                                                            >  1:11096791/1‑150 (MQ=255)
                           aCCAAACCATCGCCAATAATACACACGGCGGCATTCCGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTACTGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGccc                                                                                                                           <  2:9059061/150‑1 (MQ=255)
                           aCCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCACACGCGACAATAACACACAAAGCGACCGCATCCATACGCACTCTGccc                                                                                                                           >  1:9945726/1‑150 (MQ=255)
                                  cATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCaa                                                                                                                    >  2:4476304/1‑150 (MQ=255)
                                        aatttataaACCCGGCGGAAATCTGCACACCTTTCGGATCATAACAGGCGGAAAATACTTGGGGTACTGTCCGGGTTCGGCTTATGGCAAACGCTAAAATCACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGa                                                                                                              <  2:970544/145‑1 (MQ=255)
                                        aaataatAAACATGGCGGCAATCAGGACAACTTTTGGATCATAAAAGCCGGAAAATACTTCTGGTACTGTCAGCGTTCCGCTTAATGCAAACGCGAAAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGa                                                                                                              <  2:7061287/149‑1 (MQ=255)
                                         aataatAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGAc                                                                                                             >  2:4810104/1‑150 (MQ=255)
                                           taataaACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGa                                                                                                           >  1:4927400/1‑150 (MQ=255)
                                             ataaACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATg                                                                                                         >  1:6076281/1‑150 (MQ=255)
                                               aaaCAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGc                                                                                                       >  2:1652183/1‑150 (MQ=255)
                                               aaaCAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGc                                                                                                       >  1:21768/1‑150 (MQ=255)
                                                 aCAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCaa                                                                                                     >  2:10406110/1‑150 (MQ=255)
                                                 aCAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCaa                                                                                                     >  1:7671056/1‑150 (MQ=255)
                                                      gcggcAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGcc                                                                                                >  2:11123131/1‑150 (MQ=255)
                                                        ggcAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAg                                                                                              >  2:6277799/1‑150 (MQ=255)
                                                            aTCTGGACAAGGTTAGTATCAGAAAAGCCGTAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAAt                                                                                          <  2:8922118/150‑1 (MQ=255)
                                                             acAGGACAACTTTAGCAAAACAAAAGCCGTAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATCAAAACGCGAAAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATg                                                                                         <  2:9901473/149‑1 (MQ=255)
                                                              cAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGa                                                                                        >  1:5277623/1‑150 (MQ=255)
                                                                ggACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGaaa                                                                                      <  2:9317020/150‑1 (MQ=255)
                                                                  aCAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGa                                                                                    >  2:426184/1‑150 (MQ=255)
                                                                      ctttAGCACCAAAAAAGCCGCAAAATAATACTGTGCCTGTCCGCGTTCCTATTAATCAAAACGCGAAAAAAAAAAAAAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAAccc                                                                                <  1:5685232/148‑1 (MQ=255)
                                                                       gTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCa                                                                               >  2:3807195/1‑150 (MQ=255)
                                                                        tttggATCATAAACGCCGTAAAATACTTCTGGGACTTTCAGGTTTCCGCTTAATGAAAACGCGACAATAACAAAAAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCaa                                                                              <  2:3411767/147‑1 (MQ=255)
                                                                        ttAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCaa                                                                              >  1:784609/1‑150 (MQ=255)
                                                                          tttataaGAAAAGCCTTAAAATACTTCTGGGACGGTCAGCGTTCCTATTAATGAAAACGCGAAAATAAAAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAAt                                                                            <  2:4946714/144‑1 (MQ=255)
                                                                            gATCATAAACGCCGTAAAATACTTCTGGGAGTGTCTGGTTTCCGATTACTGAAAACGCTAAAATAAAAAACAACGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCa                                                                          <  1:927895/150‑1 (MQ=255)
                                                                            gATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCa                                                                          >  2:7527543/1‑150 (MQ=255)
                                                                            gATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCa                                                                          >  2:10343440/1‑150 (MQ=255)
                                                                              tCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAAt                                                                        >  1:4737486/1‑150 (MQ=255)
                                                                               cAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAAtt                                                                       >  2:3291585/1‑150 (MQ=255)
                                                                               cAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAAtt                                                                       >  1:3360296/1‑150 (MQ=255)
                                                                                aGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTc                                                                      >  2:10811439/1‑150 (MQ=255)
                                                                                 gAAAAGCCGTAAATTACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCa                                                                     <  2:1831568/150‑1 (MQ=255)
                                                                                 gAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCa                                                                     >  1:5124122/1‑150 (MQ=255)
                                                                                   aaaGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCAcc                                                                   >  1:1782901/1‑150 (MQ=255)
                                                                                       ccGGAAAATACTTCTGGGACTGTCAGCGATCCGCTTAATGCAAACGCGACCATAACAAACAAAGCGACCGCATCCATAAGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGAAATGAAAGAACCCAAATCAAATCACCGTTc                                                               >  1:3027404/1‑150 (MQ=255)
                                                                                        cgtaaaaTTCTTCTGGTAGTGTCAGGGTTCCTATTAATGCAAACGCTACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTca                                                              <  1:2543706/147‑1 (MQ=255)
                                                                                         ggAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTcac                                                             >  1:10411827/1‑150 (MQ=255)
                                                                                          gAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTcaca                                                            >  2:5914200/1‑150 (MQ=255)
                                                                                           aaaaTACTTCTGGTACTGTCAGCGTTCCGATTAATCAAAACGCGAAAATAACAAAAAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACaa                                                           <  1:10709843/150‑1 (MQ=255)
                                                                                                 cTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATc                                                     <  2:3540697/150‑1 (MQ=255)
                                                                                                 cTTCTGGGACTGTCAGCGTTCCGCTTAATGAAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATc                                                     <  1:8370741/150‑1 (MQ=255)
                                                                                                 cTTCTGGGACGGTCAGCTTTCCGCTAAATGCAAACCCCACAATAACAACCAACGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTacc                                                     <  1:3855315/150‑3 (MQ=255)
                                                                                                    cTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCtt                                                  >  2:8621081/1‑150 (MQ=255)
                                                                                                    cTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCtt                                                  >  2:10124069/1‑150 (MQ=255)
                                                                                                    cTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCAACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCtt                                                  >  2:2519931/1‑150 (MQ=255)
                                                                                                    cTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGAAAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCtt                                                  <  1:10011370/150‑1 (MQ=255)
                                                                                                      gggACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGt                                                <  2:3222948/150‑1 (MQ=255)
                                                                                                       ggcccgccAGCGTTCCCCTCCCTCCAACCGCGACACCCACACCCCACGCGACCGCATCCATCCGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCACATCGATTCACCGTTCACAACTTATCCTTGtt                                               <  2:6826343/143‑1 (MQ=255)
                                                                                                       ggACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGtt                                               >  1:2318492/1‑150 (MQ=255)
                                                                                                             tCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGa                                         >  2:5148693/1‑150 (MQ=255)
                                                                                                              ctgctttCCGATTAATGCAAACGCTAAAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAg                                        <  2:8040290/145‑1 (MQ=255)
                                                                                                              ccgCGGTCCGCTTAATGCAACCGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAg                                        <  2:8806497/148‑1 (MQ=255)
                                                                                                               aGCGTTCCGATTAATGCAAACGCTACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAgg                                       <  1:1643709/150‑1 (MQ=255)
                                                                                                                gCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAggg                                      >  1:806128/1‑150 (MQ=255)
                                                                                                                gCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAggg                                      <  1:6484978/150‑1 (MQ=255)
                                                                                                                   ttCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGa                                   >  1:4213965/1‑150 (MQ=255)
                                                                                                                   ttCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGa                                   >  1:404409/1‑150 (MQ=255)
                                                                                                                         cTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTCATTGAGGGGGAGGACt                              >  2:8423476/1‑150 (MQ=255)
                                                                                                                         cTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACt                              >  1:3607102/1‑150 (MQ=255)
                                                                                                                         cTTAATGCAAACGCGAAAATAAAAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACt                              <  2:7147095/150‑1 (MQ=255)
                                                                                                                            aaTGCAACCGCCAAAATAACAACCAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGa                           <  1:1055450/150‑1 (MQ=255)
                                                                                                                            aaTGCAAACGCGACAATAACAAACAACGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGa                           <  2:5677468/150‑1 (MQ=255)
                                                                                                                             aTGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTggt                          <  2:10620149/150‑3 (MQ=255)
                                                                                                                             aTGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGAt                          >  1:8151145/1‑150 (MQ=255)
                                                                                                                              tGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGAtt                         >  2:2173556/1‑150 (MQ=255)
                                                                                                                               gCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTc                        <  2:739833/150‑1 (MQ=255)
                                                                                                                                 aaaCGCGAAAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTg                      <  1:8712934/150‑1 (MQ=255)
                                                                                                                                    cgcgACAATACCAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCa                   <  1:7366664/150‑1 (MQ=255)
                                                                                                                                       gACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATaa                <  1:6886768/150‑1 (MQ=255)
                                                                                                                                          aaTAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATaaaaa             >  2:5282809/1‑150 (MQ=255)
                                                                                                                                          aaTAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATaaaaa             >  2:554737/1‑150 (MQ=255)
                                                                                                                                           aaaaaaaacaacGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATaaaaaa            <  2:2861332/145‑1 (MQ=255)
                                                                                                                                             aacaaacaCAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGc          <  2:11015732/150‑1 (MQ=255)
                                                                                                                                             aacaaacaAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGc          >  2:951194/1‑150 (MQ=255)
                                                                                                                                              acaaacaaAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGcc         >  2:6868679/1‑150 (MQ=255)
                                                                                                                                                aaacaacGCGACCGCATCCATACGCACTCGGCCCGTCGCAAACACGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGcccc       <  1:9776175/150‑1 (MQ=255)
                                                                                                                                                 aacaaaGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCcg      >  1:94780/1‑150 (MQ=255)
                                                                                                                                                  acaaaGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCTTTCACAACTTATCCTTGTTATTTGTGGGGGATGACTTGATTCTGCAATAAAAAAGCCCcga     <  1:9054972/150‑1 (MQ=255)
                                                                                                                                                   caacGCGACCGAATCAATAGGCACTCTGCCCGTCGAAAACAATAGGATGGCAAGGGCCTTTAATGAAATAACCCAAATCAATTCACCTTTCACAAATTATCCTTTTTATTTGTGGGGTATGACTTTATTCTGCCATAAAAAAGCCCcgac    <  1:4536010/150‑1 (MQ=255)
                                                                                                                                                   caaaGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCcgac    >  2:1309734/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGA  >  NC_000913/2411114‑2411411

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