Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,051,259 G→T 10.1% noncoding (224/245 nt) csrC → small regulatory RNA CsrC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,051,2590GT10.1% 591.6 / 36.2 182noncoding (224/245 nt)csrCsmall regulatory RNA CsrC
Reads supporting (aligned to +/- strand):  ref base G (75/88);  new base T (13/6);  total (88/94)
Fisher's exact test for biased strand distribution p-value = 8.87e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 189 total aligned reads displayed.

TGTCAGGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACA  >  NC_000913/4051110‑4051406
                                                                                                                                                     |                                                                                                                                                   
tGTCAGGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACagag                                                                                                                                                     <  2:5712/150‑1 (MQ=255)
 gTCAGGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGt                                                                                                                                                    >  2:1091398/1‑150 (MQ=255)
 gTCAGGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGt                                                                                                                                                    >  1:326863/1‑150 (MQ=255)
   cAGGAACACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTaa                                                                                                                                                  >  2:669627/1‑150 (MQ=255)
    aGGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAAt                                                                                                                                                 >  2:47787/1‑150 (MQ=255)
     ggAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATc                                                                                                                                                >  1:1356646/1‑150 (MQ=255)
     ggAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCAGGCTAAAGGAAAAACTGGGTGTTTTGGCGGCCTGCAATGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATc                                                                                                                                                <  1:586745/150‑1 (MQ=255)
       aaGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTg                                                                                                                                              <  2:969023/150‑1 (MQ=255)
         gACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTc                                                                                                                                            >  1:457306/1‑150 (MQ=255)
            aaaCTTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGCCt                                                                                                                                        <  1:610393/150‑1 (MQ=255)
            aaaCGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                                                                         <  1:303244/150‑1 (MQ=255)
             aaCGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCt                                                                                                                                        >  1:1295946/1‑150 (MQ=255)
               cGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCttt                                                                                                                                      >  1:79705/1‑150 (MQ=255)
                   cGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTc                                                                                                                                  <  2:873903/150‑1 (MQ=255)
                    ggAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCATTTTTCt                                                                                                                                 <  1:1007340/150‑1 (MQ=255)
                       gaCGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                                                                              <  2:257558/150‑1 (MQ=255)
                       gaCGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                                                                              <  2:437227/150‑1 (MQ=255)
                        aCGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgc                                                                                                                             <  2:232002/150‑1 (MQ=255)
                         cGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgct                                                                                                                            >  2:200046/1‑150 (MQ=255)
                         cGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgct                                                                                                                            <  2:386715/150‑1 (MQ=255)
                         cGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgct                                                                                                                            >  1:1398228/1‑150 (MQ=255)
                          gTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgctt                                                                                                                           <  1:827359/150‑1 (MQ=255)
                             aTTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTtgcttgct                                                                                                                        <  2:1301369/150‑1 (MQ=255)
                             aTTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTgcttgctt                                                                                                                       <  1:769424/150‑1 (MQ=255)
                              ttAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTgcttgctt                                                                                                                       <  1:669627/150‑1 (MQ=255)
                              ttAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTgcttgctt                                                                                                                       >  1:638383/1‑150 (MQ=255)
                                aaaCGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTTGCCTTTTTTCTTTGCTTGCTTTc                                                                                                                     >  1:898690/1‑150 (MQ=255)
                                aaaCGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTc                                                                                                                     <  2:79705/150‑1 (MQ=255)
                                     gAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                                                                <  2:1398228/150‑1 (MQ=255)
                                     gAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                                                                <  1:47787/150‑1 (MQ=255)
                                     gAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                                                                <  1:505838/150‑1 (MQ=255)
                                        aTGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGAt                                                                                                             <  1:1182918/150‑1 (MQ=255)
                                                aaCACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCcgccgc                                                                                                     >  2:82919/1‑150 (MQ=255)
                                                     ggATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATc                                                                                                <  1:431008/150‑1 (MQ=255)
                                                     ggATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCGGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATc                                                                                                <  2:1022979/150‑1 (MQ=255)
                                                       aTTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCta                                                                                              <  2:1390038/150‑1 (MQ=255)
                                                       aTTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCta                                                                                              <  1:629758/150‑1 (MQ=255)
                                                          gTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATAc                                                                                           <  2:416252/150‑1 (MQ=255)
                                                          gTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATAc                                                                                           >  1:820397/1‑150 (MQ=255)
                                                              cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGaa                                                                                       <  1:1379403/150‑1 (MQ=255)
                                                                gCTAACTGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATa                                                                                     <  2:401940/150‑1 (MQ=255)
                                                                    aaGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAAc                                                                                 <  2:1054579/150‑1 (MQ=255)
                                                                      ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACgg                                                                               >  1:1132178/1‑150 (MQ=255)
                                                                         aaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTa                                                                            <  1:698583/150‑1 (MQ=255)
                                                                           aaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTaaa                                                                          >  1:280840/1‑150 (MQ=255)
                                                                             cAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAg                                                                        <  2:574115/150‑1 (MQ=255)
                                                                             cAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAg                                                                        <  2:947875/150‑1 (MQ=255)
                                                                               gggTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTg                                                                                         <  2:396264/131‑1 (MQ=255)
                                                                                 gtgtgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGa                                                                    <  2:79219/150‑1 (MQ=255)
                                                                                  tgtgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGaa                                                                   <  1:388265/150‑1 (MQ=255)
                                                                                     gtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCa                                                                <  1:200046/150‑1 (MQ=255)
                                                                                     gtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCa                                                                >  2:773175/1‑150 (MQ=255)
                                                                                       tGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCAtt                                                              <  1:851417/150‑1 (MQ=255)
                                                                                       tGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCAtt                                                              <  1:679584/150‑1 (MQ=255)
                                                                                        ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa                                                             >  2:1111516/1‑150 (MQ=255)
                                                                                        ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa                                                             <  2:543616/150‑1 (MQ=255)
                                                                                        ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa                                                             <  2:588732/150‑1 (MQ=255)
                                                                                        ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa                                                             <  1:75547/150‑1 (MQ=255)
                                                                                          cggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg                                                           <  2:1179820/150‑1 (MQ=255)
                                                                                          cggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg                                                           <  2:1295946/150‑1 (MQ=255)
                                                                                             ccTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGaaa                                                        >  1:397306/1‑150 (MQ=255)
                                                                                              cTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAc                                                       <  1:1179486/150‑1 (MQ=255)
                                                                                                gCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACca                                                     >  2:949785/1‑150 (MQ=255)
                                                                                                 cAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcat                                                    <  2:337694/150‑1 (MQ=255)
                                                                                                  aaGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAATAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatc                                                   <  2:1339740/150‑1 (MQ=255)
                                                                                                  aaGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa                                                             <  1:532769/140‑1 (MQ=255)
                                                                                                  aaGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa                                                             >  2:532769/1‑140 (MQ=255)
                                                                                                  aaGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatc                                                   >  2:1251421/1‑150 (MQ=255)
                                                                                                      aTTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCt                                               >  1:85701/1‑150 (MQ=255)
                                                                                                        tGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTc                                             <  1:1022977/150‑1 (MQ=255)
                                                                                                        tGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTc                                             <  1:1091398/150‑1 (MQ=255)
                                                                                                          tAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCaa                                           <  2:845275/150‑1 (MQ=255)
                                                                                                            aGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAAc                                         >  1:180395/1‑150 (MQ=255)
                                                                                                             gACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACt                                        <  1:415935/150‑1 (MQ=255)
                                                                                                             gACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACt                                        >  1:1066808/1‑150 (MQ=255)
                                                                                                                  gTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAg                                                                        <  2:358258/113‑1 (MQ=255)
                                                                                                                   ttAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACgca                                  <  2:179511/150‑1 (MQ=255)
                                                                                                                       gggTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCaa                              >  1:1331130/1‑150 (MQ=255)
                                                                                                                           tATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGt                          >  2:739808/1‑150 (MQ=255)
                                                                                                                            aTGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTc                         <  2:1356646/150‑1 (MQ=255)
                                                                                                                               aGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATg                      <  1:949221/150‑1 (MQ=255)
                                                                                                                               aGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATg                      <  1:1389259/150‑1 (MQ=255)
                                                                                                                                gcaaagAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGc                     >  2:1005536/6‑150 (MQ=255)
                                                                                                                                gTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGc                     >  1:1334105/1‑150 (MQ=255)
                                                                                                                                   aagAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCaaa                  >  1:667847/3‑150 (MQ=255)
                                                                                                                                   aagAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCaaa                  >  2:746206/3‑150 (MQ=255)
                                                                                                                                   aGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCaaa                  >  2:59577/1‑150 (MQ=255)
                                                                                                                                   aGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCaaa                  <  1:20853/150‑1 (MQ=255)
                                                                                                                                    ggAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCaaaa                 >  1:443593/1‑150 (MQ=255)
                                                                                                                                    ggAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCaaaa                 >  1:490240/1‑150 (MQ=255)
                                                                                                                                    ggAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCaaaa                 >  2:165031/1‑150 (MQ=255)
                                                                                                                                     gAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAg                >  1:837989/1‑150 (MQ=255)
                                                                                                                                         aaaGGCGACAGATTACTCTGTCGCCTTTTTTATTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGcgcg            >  2:1021607/1‑150 (MQ=255)
                                                                                                                                          aaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGt           >  2:416939/1‑150 (MQ=255)
                                                                                                                                          aaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGt           >  2:386436/1‑150 (MQ=255)
                                                                                                                                            ggCGACAGATTACTCTGTCGCATTTTTTCTTTGCTTGCTTTCTCTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGAAGAACCATTATCAAACCATCATATACAAACTCACGCAGCAAAGGTCATCCAAAAGcgcgacc         >  2:388717/1‑147 (MQ=255)
                                                                                                                                             gCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTc          >  2:708474/1‑148 (MQ=255)
                                                                                                                                              cGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGaa       >  2:1264136/1‑150 (MQ=255)
                                                                                                                                               gACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGaaa      <  1:1251421/150‑1 (MQ=255)
                                                                                                                                               caCAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGaaa      <  1:340193/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGTCAGGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACA  >  NC_000913/4051110‑4051406

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