Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,049,986 G→T 8.2% noncoding (88/109 nt) spf → small regulatory RNA Spot 42

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,049,9860GT 8.2% 601.5 / 24.3 180noncoding (88/109 nt)spfsmall regulatory RNA Spot 42
Reads supporting (aligned to +/- strand):  ref base G (85/80);  new base T (8/7);  total (93/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

Only 100 of 186 total aligned reads displayed.

CTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTA  >  NC_000913/4049839‑4050135
                                                                                                                                                   |                                                                                                                                                     
cTGAAAAAAAATTACAAAATGTGCTTTCTGCACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcagga                                                                                                                                                     >  2:2402574/1‑148 (MQ=255)
cTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTa                                                                                                                                                     >  1:765958/1‑150 (MQ=255)
cTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTa                                                                                                                                                     >  2:535178/1‑150 (MQ=255)
  gAAAAAAAATTACAAAAAGTCCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAAt                                                                                                                                                   <  1:1556984/150‑1 (MQ=255)
   aaaaaaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCGGTAATg                                                                                                                                                  >  2:2360705/1‑150 (MQ=255)
   aaaaaaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATg                                                                                                                                                  <  1:2476057/150‑1 (MQ=255)
    aaaaaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGa                                                                                                                                                 >  2:1785795/1‑150 (MQ=255)
     aaaaaaTTACAAAAAGTGCTTACTGAACTGAACAAAAAAGAGTAAAGTTAGTCCCGTAGGCTACAGAGGTAAGATCTTCTATCTTTCAGACCTTTTACTGAACGTAAGCCGAATTGGCTGAATATTTTAGCCGCCCCAGGCAGTAATGAc                                                                                                                                                >  2:196823/1‑150 (MQ=255)
      aaaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACt                                                                                                                                               >  1:349600/1‑150 (MQ=255)
          aTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTgggg                                                                                                                                           <  2:1748161/150‑1 (MQ=255)
            tACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCg                                                                                                                                         <  2:852079/150‑1 (MQ=255)
                aaaaGTGCTTTCTTAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttt                                                                                                                                     <  1:1129897/150‑1 (MQ=255)
                  aaGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGGCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCACTCAGTAATGACTGGGGCGtttttt                                                                                                                                   <  1:680272/150‑1 (MQ=255)
                   aGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTa                                                                                                                                  <  1:2326264/150‑1 (MQ=255)
                   aGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTa                                                                                                                                  <  2:1928050/150‑1 (MQ=255)
                      gCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTg                                                                                                                               <  1:2218365/150‑1 (MQ=255)
                       cTTTCTGAACTGAACAAAAAAGAGAAAAGTTAGTCGCGTAGGGTAAAGAGGTAAGATGTTATATATTTAACACCTTTTACTAAACGTAAAAGTATTTGGATGAATATTTTACCCGCCCCATACAGTAATTAATGGGGCGTTTTTTATTgg                                                                                                                              >  2:530710/1‑150 (MQ=255)
                        tttCTGAACTTAACAAAAAAGACTAAAGTTAGTCGCGTAGGTTACAGAGGGAAGCTGTGCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTggc                                                                                                                             <  1:648818/150‑2 (MQ=255)
                        tttCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGGACAGAGGGAAGATGTGCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTggg                                                                                                                             <  1:2433437/150‑1 (MQ=255)
                            tgaactgaacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaa                                                                                                                         <  2:1403982/150‑1 (MQ=255)
                                ctcaacAAAAAAGAGTAAAGTTAGTCGCCTAGGGTACAGAGCTAAGATGTTGTATCTTACAGACCATAAAATTCACGTTATAGGATTTGGCTGAATTTTTTAGACGCCCCAGACAGTAATGACTGGGGCGTTTTTTTATGGGCgatagaa                                                                                                                     >  2:262387/4‑150 (MQ=255)
                                 tgaacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaagaaa                                                                                                                    <  1:2021819/150‑1 (MQ=255)
                                  gaacAATAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaa                                                                                                                   >  1:796399/1‑150 (MQ=255)
                                  gaacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaa                                                                                                                   >  2:1435031/1‑150 (MQ=255)
                                   aacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaag                                                                                                                  <  1:364885/150‑1 (MQ=255)
                                   aacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaag                                                                                                                  <  1:1814611/150‑1 (MQ=255)
                                    acAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaaga                                                                                                                 >  1:1257570/1‑150 (MQ=255)
                                     cAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGAt                                                                                                                <  2:1315300/150‑1 (MQ=255)
                                     cAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGAt                                                                                                                >  1:1757462/1‑150 (MQ=255)
                                       aaaaaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATcc                                                                                                              <  1:615693/150‑1 (MQ=255)
                                       aaaaaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATcc                                                                                                              <  2:1496033/150‑1 (MQ=255)
                                        aaaaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCg                                                                                                             <  1:1832286/150‑1 (MQ=255)
                                         aaaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGt                                                                                                            <  1:2154812/150‑1 (MQ=255)
                                              gTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATATTTCACACCTTTTCCTTGCCGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACAGGGGCGTTTTTTATTGGGCGATAGAAAAGCTCCGTAATGc                                                                                                       >  2:1322557/1‑150 (MQ=255)
                                               tAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGcc                                                                                                      <  1:1400472/150‑1 (MQ=255)
                                                aaaGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCt                                                                                                     <  2:226688/150‑1 (MQ=255)
                                                   gTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGAt                                                                                                  >  2:1211891/1‑150 (MQ=255)
                                                       gTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcgc                                                                                              >  2:2286857/1‑150 (MQ=255)
                                                       gTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcgc                                                                                              <  1:1897570/150‑1 (MQ=255)
                                                               gggTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTa                                                                                      >  1:1972265/1‑150 (MQ=255)
                                                               gggTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTa                                                                                      <  1:1193373/150‑1 (MQ=255)
                                                                 gTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATc                                                                                    <  2:2450132/150‑1 (MQ=255)
                                                                    cAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAgg                                                                                 <  1:535178/150‑1 (MQ=255)
                                                                      gagGTAAGATGTTCTATCTTTCGGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGcc                                                                               <  2:894533/150‑1 (MQ=255)
                                                                      gagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGcc                                                                               >  2:593170/1‑150 (MQ=255)
                                                                      gagGAAAGATGTTCTATCTATCAGACAATTTACTTCACGTAATCGGATTTCGCAGAATATTTTAACCCCCACAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGGGCTATGTTTATCAGGcc                                                                               >  2:1723241/1‑150 (MQ=255)
                                                                        ggTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCaa                                                                             >  2:1008135/1‑150 (MQ=255)
                                                                            aGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCGGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGt                                                                         <  1:2360705/150‑1 (MQ=255)
                                                                             gATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTa                                                                        <  1:18475/150‑1 (MQ=255)
                                                                             gATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGACGCCACAGTCATTAATGACTGATGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGACTGATGGGGTATGTTTATAAGGACAACGGTa                                                                        >  1:2061839/1‑150 (MQ=255)
                                                                              aTGTTCTATCTTTAAGACCTTTTACTTCACGTAATCGTATTTGTCTTAATATTTTATCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCTAAAGAAAATATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACTGTAg                                                                       <  2:870708/150‑1 (MQ=255)
                                                                               tGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGa                                                                      >  2:487531/1‑150 (MQ=255)
                                                                                gTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGaa                                                                     >  2:2249105/1‑150 (MQ=255)
                                                                                  tCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAAtt                                                                   >  1:2040781/1‑150 (MQ=255)
                                                                                  tCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAAtt                                                                   <  2:349600/150‑1 (MQ=255)
                                                                                   cTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTg                                                                  >  1:704854/1‑150 (MQ=255)
                                                                                   cTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGGCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATCCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTg                                                                  <  2:2214225/150‑1 (MQ=255)
                                                                                     aTCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTa                                                                >  1:929802/1‑150 (MQ=255)
                                                                                        tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc                                                             >  1:977603/1‑150 (MQ=255)
                                                                                          tCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTa                                                           >  1:1378898/1‑150 (MQ=255)
                                                                                          tCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTa                                                           >  2:160679/1‑150 (MQ=255)
                                                                                               ccTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGaa                                                      <  1:1589058/150‑1 (MQ=255)
                                                                                               ccTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGaa                                                      >  1:2137506/1‑150 (MQ=255)
                                                                                                 ttttACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAAtt                                                    >  1:1909926/1‑150 (MQ=255)
                                                                                                   ttACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTa                                                  <  2:66072/150‑1 (MQ=255)
                                                                                                     aCTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCAAACGGTAGAATTGTAATCTATTGAATTTACg                                                >  1:1792696/1‑150 (MQ=255)
                                                                                                        tCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGc                                             <  1:1008135/150‑1 (MQ=255)
                                                                                                         cACGTAATCGGATTTGACAGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCAAACGGTAGAATTGTAATCTATTGAATTTACGGGcc                                            >  2:73573/1‑150 (MQ=255)
                                                                                                          aagTAATCGGATTTGGATGAATATTTTAGCCGCCCCATTCATTAATTAATTTGGCTTTTTTTATTTGTCGAAAGAAAAGATCCTTAATGACTGATGCGCTATGTTTATAAGTCCAACGTTATAATTGTAATCTATTGAATTTACgtgact                                           <  2:2136633/148‑6 (MQ=255)
                                                                                                               aTCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATAc                                      <  1:1435031/150‑1 (MQ=255)
                                                                                                                  ggATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGcc                                   <  2:1571002/150‑1 (MQ=255)
                                                                                                                  ggATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGcc                                   <  1:1350779/150‑1 (MQ=255)
                                                                                                                  ggATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGcc                                   <  1:1503216/150‑1 (MQ=255)
                                                                                                                        ggCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATcc                             >  2:1353495/1‑150 (MQ=255)
                                                                                                                         gCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCg                            >  2:1528747/1‑150 (MQ=255)
                                                                                                                         gCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCg                            >  2:1040571/1‑150 (MQ=255)
                                                                                                                           ttaaTATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATTCTCTATGTTTATCAGTCCAACTGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGc                          <  2:765958/148‑1 (MQ=255)
                                                                                                                           tGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGc                          <  1:1364392/150‑1 (MQ=255)
                                                                                                                            gAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGca                         >  1:1984152/1‑150 (MQ=255)
                                                                                                                              atatTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGcaca                       <  2:458170/150‑1 (MQ=255)
                                                                                                                                 ttttccccgccCCAGTCAGTCAGGCCGGGGGCGTTTTTTATTGGGCGAACGAAAACATCCGTAATGCCTGAGGCGCTAGGGTTCTCAGGCCAACGGTAGAATTTTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGc                    <  2:1102953/144‑1 (MQ=255)
                                                                                                                                  tttAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCa                   <  1:720601/150‑1 (MQ=255)
                                                                                                                                    tAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCAtt                 >  2:1659416/1‑150 (MQ=255)
                                                                                                                                     aGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTa                >  2:1245926/1‑150 (MQ=255)
                                                                                                                                      gccgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGTGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTaa               >  1:1666141/1‑150 (MQ=255)
                                                                                                                                      gccgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTaa               <  2:1294333/150‑1 (MQ=255)
                                                                                                                                      gccgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTaa               >  2:725388/1‑150 (MQ=255)
                                                                                                                                       ccgccCCATTCAGTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATCCGCCACATGCGGCACAAGCATTAAg              <  2:2031499/150‑1 (MQ=255)
                                                                                                                                        cgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAgg             >  2:550888/1‑150 (MQ=255)
                                                                                                                                        cgccACTGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAgg             >  2:960457/1‑150 (MQ=255)
                                                                                                                                         gccCCAGTCATTACTGCCTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCCGCACAAGCATTAAGGc            >  2:568776/1‑150 (MQ=255)
                                                                                                                                         gccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGc            >  2:1778123/1‑150 (MQ=255)
                                                                                                                                          ccccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACCAGCATTAAGGCa           >  1:44788/1‑150 (MQ=255)
                                                                                                                                          ccccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCa           >  1:40013/1‑150 (MQ=255)
                                                                                                                                          ccccAGGCAGTAATGACTGGGGCGTTTTTTATTGGGCGACAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCa           <  1:1700981/150‑1 (MQ=255)
                                                                                                                                           cccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCaa          >  2:1057197/1‑150 (MQ=255)
                                                                                                                                           cccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCaa          >  1:1006410/1‑150 (MQ=255)
                                                                                                                                                 cagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAat    <  2:497469/150‑1 (MQ=255)
                                                                                                                                                 cagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAat    <  1:1040571/150‑1 (MQ=255)
                                                                                                                                                   gtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAatta  <  1:487531/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTA  >  NC_000913/4049839‑4050135

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