Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,722,301 A→C 5.8% intergenic (+196/+27) sokA → / ← glyS pseudogene, sRNA antisense regulator blocking mokA/hokA translation/glycine tRNA synthetase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,722,3010AC 5.8% 567.4 / 18.6 187intergenic (+196/+27)sokA/glySpseudogene, sRNA antisense regulator blocking mokA/hokA translation/glycine tRNA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (90/86);  new base C (7/4);  total (97/90)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.25e-01

Only 100 of 191 total aligned reads displayed.

CCCGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCGAAGAAAGC  >  NC_000913/3722152‑3722450
                                                                                                                                                      |                                                                                                                                                     
cccGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGa                                                                                                                                                       >  2:1292666/1‑150 (MQ=255)
cccGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGACTGCCAGa                                                                                                                                                       >  1:568074/1‑150 (MQ=255)
   gTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATgg                                                                                                                                                    <  2:1136559/150‑1 (MQ=255)
    tGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGTCAGATGGc                                                                                                                                                   <  1:979668/150‑1 (MQ=255)
     ggCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGACAGATGGCa                                                                                                                                                  >  1:259511/1‑150 (MQ=255)
        tATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGc                                                                                                                                               <  2:923731/150‑1 (MQ=255)
         aTTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCt                                                                                                                                              >  1:1502957/1‑150 (MQ=255)
          ttCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTa                                                                                                                                             >  1:1543682/1‑150 (MQ=255)
           tCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTAATCAGAAGAATATTAAAGGTTCCTGAATTATATAGCGATAAAAAAAGCCTGCCAGAAGTCAGGCAAt                                                                                                                                            >  1:120754/1‑150 (MQ=255)
             tAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCATCTGGCAGGCTtttt                                                                                                                                          <  1:444349/150‑1 (MQ=255)
              aGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTta                                                                                                                                         >  2:793241/1‑150 (MQ=255)
               gCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTtaa                                                                                                                                        >  2:674453/1‑150 (MQ=255)
               gCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTtaa                                                                                                                                        >  1:1373880/1‑150 (MQ=255)
               gCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTtaa                                                                                                                                        <  2:1574992/150‑1 (MQ=255)
               gCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTtaa                                                                                                                                        <  2:446288/150‑1 (MQ=255)
                 ttAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCATCTGGCAGGCTttttttat                                                                                                                                      <  2:636097/150‑5 (MQ=255)
                    ataaGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCATCTGGCAGGCTttttttatcgc                                                                                                                                   <  2:1384533/150‑8 (MQ=255)
                      aaGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGc                                                                                                                                 >  1:884097/1‑150 (MQ=255)
                          ttgtttCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGtta                                                                                                                             >  2:179941/1‑150 (MQ=255)
                           tgtttCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGttat                                                                                                                            <  2:53739/150‑1 (MQ=255)
                           tgtttCTTATACTTGATATAAACCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGTTTGATGTATTATTTATAGATAAAAAAAGCATGCCATATGGAATGCTATATAATAACTGCTttat                                                                                                                            >  2:1367007/1‑150 (MQ=255)
                            gtttCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGttatt                                                                                                                           >  2:1672417/1‑150 (MQ=255)
                             tttCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGAGATAAAAAAATCCTGCCAGATGGCAGGCTATTTAATAACGGCGttatta                                                                                                                          >  1:982102/1‑150 (MQ=255)
                              ttcttcTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTtattat                                                                                                                         <  2:581841/150‑1 (MQ=255)
                               tcttcttctTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACgg                                                                                                                                  >  2:1422703/1‑140 (MQ=255)
                               tcttctCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACTGCGTTattatt                                                                                                                        >  1:1029115/1‑150 (MQ=255)
                                 ttctCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGGTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGc                                                                                                                      <  2:169227/150‑1 (MQ=255)
                                 ttctCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGc                                                                                                                      <  2:1560039/150‑1 (MQ=255)
                                 ttctCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGc                                                                                                                      >  2:1479391/1‑150 (MQ=255)
                                 ttctCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGc                                                                                                                      <  1:1644780/150‑1 (MQ=255)
                                 ttctCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGc                                                                                                                      <  1:236739/150‑1 (MQ=255)
                                 ttctCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGAGGGCAGGCTATTTAATAACGGCGTTATTATTGc                                                                                                                      <  1:119390/150‑1 (MQ=255)
                                  tctcCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCa                                                                                                                     <  2:1172655/150‑1 (MQ=255)
                                  tctcCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCa                                                                                                                     <  2:303494/150‑1 (MQ=255)
                                   ctcCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGATGATTATTAAAGGTTGGTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATACAGGCGTTATTATTGCaa                                                                                                                    >  1:113000/1‑150 (MQ=255)
                                   ctcCTTGATATAATCCTAAAAAAATCTCATAAAATAAATATATGAGAAAATCTATATTCAGCAGAAGATTATTAAAGTTTGCTATATTATTTAGCTATAAAAAAAGCCTGACAGATGGCAGGCTATTTAATAACGGCGTTATTATTGaaa                                                                                                                    >  1:199400/1‑150 (MQ=255)
                                   atcCTTGATATAATACTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGCCGttattaaatagc                                                                                                                    <  2:1321622/149‑7 (MQ=255)
                                     ccTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACa                                                                                                                  >  2:715944/1‑150 (MQ=255)
                                     ccTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACa                                                                                                                  >  2:829596/1‑150 (MQ=255)
                                     ccTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACa                                                                                                                  >  1:99846/1‑150 (MQ=255)
                                     ccTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACa                                                                                                                  <  1:1176511/150‑1 (MQ=255)
                                     ccTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACa                                                                                                                  >  1:368175/1‑150 (MQ=255)
                                     ccTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTAAGAAGAAGATTATTAAAGGTTGCTGTATTATATAGCGATAAAAAAAGCCTGCCAGATTGCAGGCTATTTAATAACTGCGTTTTTATTGCAACa                                                                                                                  >  2:1331316/1‑150 (MQ=255)
                                        tGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCTATAAAAAAAGCATGCCAGATGGCAGGCTATTTAATAACTGCGTTATTATTGCAACAGcc                                                                                                               >  1:540831/1‑149 (MQ=255)
                                        tGATATAATCCTAAAAAAATCACATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCTATAAAAAAAGCCTGCCAGATGACAGGCTATTTAATAACGGCGTTATTATTGAACAGCGa                                                                                                              >  2:1607054/1‑150 (MQ=255)
                                                tCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATcc                                                                                                       <  1:698506/150‑1 (MQ=255)
                                                  cTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGc                                                                                                     <  1:179941/150‑1 (MQ=255)
                                                  cTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGc                                                                                                     >  2:1205790/1‑150 (MQ=255)
                                                   tAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCa                                                                                                    <  1:1479391/150‑1 (MQ=255)
                                                     aaaaaaTCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGATGTATTATTTAGCGATAAAAAAAAGCATGCCAGATGGAAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCaa                                                                                                   >  2:913310/1‑150 (MQ=255)
                                                         aaTCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACgcgc                                                                                              <  2:1007313/150‑1 (MQ=255)
                                                            ctcATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAgg                                                                                           >  1:1081555/1‑150 (MQ=255)
                                                              cATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGaa                                                                                         >  1:1002317/1‑150 (MQ=255)
                                                                 aaaaTTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAg                                                                                      <  2:442042/150‑1 (MQ=255)
                                                                     ttAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTcg                                                                                  <  2:368175/150‑1 (MQ=255)
                                                                        atatatGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCgcgc                                                                               >  1:742150/1‑150 (MQ=255)
                                                                            atGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGtt                                                                           <  1:1193323/150‑1 (MQ=255)
                                                                              gagaTAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGtttg                                                                         >  2:1157845/1‑149 (MQ=255)
                                                                                  tAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCa                                                                     <  2:1296308/150‑1 (MQ=255)
                                                                                        ttATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGG                                                               <  2:1649855‑M1/150‑1 (MQ=255)
                                                                                          aTTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCAAGATGGCAGGCTATTTAATAACGGCGTTATTATTTCAACAGCGAAATATCCGCAACGCGCAGTAAAAGTTCGCGCAGTTTCTCCATCATGGTC                                                             >  1:1043384‑M1/1‑150 (MQ=255)
                                                                                           ttCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCA                                                            >  1:1589242‑M1/1‑150 (MQ=255)
                                                                                           ttCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCA                                                            <  1:1597946‑M1/150‑1 (MQ=255)
                                                                                                     gaTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGAt                                                  <  2:553861/150‑1 (MQ=255)
                                                                                                     gaTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGAt                                                  <  1:1685566/150‑1 (MQ=255)
                                                                                                     gaTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGAt                                                  <  1:1659614/150‑1 (MQ=255)
                                                                                                                  tGCTGTATTATTTAGCGATAAAAATAGCCTGCCAGATGGCAGTCTATTTAATAACTGCGTTATTATTGCAACAGCGAAATATCAGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTg                                     <  2:1226334/150‑1 (MQ=255)
                                                                                                                  tGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTAATACGCAATTCTTTg                                     <  2:972533/150‑1 (MQ=255)
                                                                                                                    cTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGtc                                   >  1:1230528/1‑150 (MQ=255)
                                                                                                                    cTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGGCAGACGGTTGATACGCAATTCTTTGtc                                   <  2:205596/150‑1 (MQ=255)
                                                                                                                      gTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGGAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGtcat                                 >  1:133304/1‑150 (MQ=255)
                                                                                                                      gTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGtcat                                 >  1:1053702/1‑150 (MQ=255)
                                                                                                                        attattTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGtcatca                               <  1:968159/150‑1 (MQ=255)
                                                                                                                           attTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAAcc                            <  2:936055/150‑1 (MQ=255)
                                                                                                                              tAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGAGAAATATCCGCAACGCGCAGGAAAAGTTCTCGCAGTTTCTCCAGCATTGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCAtg                         <  1:1039390‑M2/150‑3 (MQ=255)
                                                                                                                                gCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCAtgac                       >  2:1549993‑M2/1‑146 (MQ=255)
                                                                                                                                gCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCAtgac                       >  2:259511‑M2/1‑146 (MQ=255)
                                                                                                                                 cGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCAtgacc                      >  2:344903‑M2/1‑145 (MQ=255)
                                                                                                                                  gATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAtggttgatgacaaagaattgcgtatcaaccgtttgacca                     >  1:390427‑M1/1‑111 (MQ=255)
                                                                                                                                  gATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCAtgacca                     <  1:1555132‑M2/150‑7 (MQ=255)
                                                                                                                                  gATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCAtgacca                     >  2:1496831‑M2/1‑144 (MQ=255)
                                                                                                                                  gATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCAtgacca                     >  1:674453‑M2/1‑144 (MQ=255)
                                                                                                                                  gATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTACAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCAtgacca                     <  2:1082959‑M2/150‑7 (MQ=255)
                                                                                                                                   aTAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAtggttgatgacaaagaattgcgtatcaaccgtctgaccat                    <  2:1448888‑M1/150‑41 (MQ=255)
                                                                                                                                    tAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGCAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATc                   >  1:1159162/1‑150 (MQ=255)
                                                                                                                                    tAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATc                   <  1:100395/150‑1 (MQ=255)
                                                                                                                                      ttaaataGCCTGCCATCTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCa                  >  2:491392/3‑150 (MQ=255)
                                                                                                                                       aaaaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCAc                 >  2:1317528/1‑150 (MQ=255)
                                                                                                                                        aaataGCCTGCCATCTGGCCGGCTATTTAATTAGGGCGGTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAGCAGTTCGCGCAGTTTCGCCAGCATGGTCTGACGGTTGAGACGCAATGCTTTGTCGTCAACCATGACCATCACt                >  1:1626073/1‑150 (MQ=255)
                                                                                                                                         aaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAtggttgatgacaaagaattgcgtatcaaccgtttgaccatgctgg               <  2:390427‑M1/150‑46 (MQ=255)
                                                                                                                                             gCCTGCCATATGGCATGCTATTTAATAACGGCGTTATTATTGCAACATCGAAATATCCGCAACGAGCAGGAACAGTTCGCGCAGTTTCTCCATTATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTTTc           <  1:913310/150‑1 (MQ=255)
                                                                                                                                              ccTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACTGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATca          <  2:982102/150‑2 (MQ=255)
                                                                                                                                                tGCCATCTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCgaa        >  2:633227/1‑150 (MQ=255)
                                                                                                                                                tGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCgaa        >  2:404273/1‑150 (MQ=255)
                                                                                                                                                 gCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCgaag       <  1:807129/150‑1 (MQ=255)
                                                                                                                                                  ccAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGCCCATCACTTTATCgaaga      >  2:1058265/1‑150 (MQ=255)
                                                                                                                                                  ccAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCgaaga      >  2:187917/1‑150 (MQ=255)
                                                                                                                                                  ccAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCgaaga      >  2:882631/1‑150 (MQ=255)
                                                                                                                                                  ccAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTACAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCgaaga      >  1:77342/1‑150 (MQ=255)
                                                                                                                                                      aTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTGGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCGAAGAAAGc  <  1:1120602/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCCGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCGAAGAAAGC  >  NC_000913/3722152‑3722450

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