Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,704,304 A→G 5.7% intergenic (‑241/+31) rnc ← / ← lepB RNase III/leader peptidase (signal peptidase I)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,704,3040AG 5.7% 362.7 / 9.1 123intergenic (‑241/+31)rnc/lepBRNase III/leader peptidase (signal peptidase I)
Reads supporting (aligned to +/- strand):  ref base A (54/62);  new base G (3/4);  total (57/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

Only 100 of 123 total aligned reads displayed.

TACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAA  >  NC_000913/2704156‑2704451
                                                                                                                                                    |                                                                                                                                                   
tACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGACAt                                                                                                                                                    >  2:676319/1‑150 (MQ=255)
   cGGTATATGGTGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGTGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAAc                                                                                                                                                 <  1:1020816/150‑1 (MQ=255)
    ggTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACg                                                                                                                                                >  2:252989/1‑150 (MQ=255)
    ggTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAAACTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTGAGCCCCGGGAGCTTTCTCTCATAAATAATTCACGTTGTCGCCATAACg                                                                                                                                                >  1:263913/1‑150 (MQ=255)
     gTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACgg                                                                                                                                               <  2:221080/150‑1 (MQ=255)
         atGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGAc                                                                                                                                           <  1:1383973/150‑1 (MQ=255)
         atGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGAc                                                                                                                                           >  2:873809/1‑150 (MQ=255)
           ggggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACaa                                                                                                                                         >  1:445955/1‑150 (MQ=255)
             ggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACg                                                                                                                                       <  2:532764/150‑1 (MQ=255)
             ggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACg                                                                                                                                       <  2:792175/150‑1 (MQ=255)
               aCCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACGTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACGTg                                                                                                                                     <  2:1130308/150‑1 (MQ=255)
               aCCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTg                                                                                                                                     <  1:1502977/150‑1 (MQ=255)
                ccAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGCCGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGa                                                                                                                                    >  1:1091655/1‑150 (MQ=255)
                  aaCGCGGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACGTGaat                                                                                                                                  <  2:1575605/150‑2 (MQ=255)
                    cGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGa                                                                                                                                >  2:64607/1‑150 (MQ=255)
                       tGTTTCGTGTGCCGTTGCAACCTGGACGTGCCAACCTGAAACTTCGGGCGAATATTCTATACACAACGACGGGGGAAGTCGTTAGCCACGGGAGAGTTATCTCAAAAAAAATACAAGTATTCCCCAAAACGAGGACAACCTGAACGAAGa                                                                                                                             >  2:116122/1‑150 (MQ=255)
                       tGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAACCTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGa                                                                                                                             <  1:1306503/150‑1 (MQ=255)
                       tGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGa                                                                                                                             >  2:359442/1‑150 (MQ=255)
                             gtgtGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAAAGAAGATGGCTa                                                                                                                       >  1:1220583/1‑150 (MQ=255)
                                tgCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTa                                                                                                                    >  1:1216985/1‑150 (MQ=255)
                                tgCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTa                                                                                                                    >  2:1115182/1‑150 (MQ=255)
                                 gCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTaa                                                                                                                   >  1:207746/1‑150 (MQ=255)
                                     tGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGa                                                                                                               <  1:1607154/150‑1 (MQ=255)
                                     tGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGa                                                                                                               <  2:566497/150‑1 (MQ=255)
                                          aCCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgccg                                                                                                          <  1:133655/150‑1 (MQ=255)
                                           ccTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgccgc                                                                                                         >  2:543199/1‑150 (MQ=255)
                                            cTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGCCGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgccgcc                                                                                                        <  1:1332433/150‑1 (MQ=255)
                                             tGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCAACGGGAGATTTTTCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAAAGTGAACGAAGATGGCTATTATTGGATGCGGCCa                                                                                                       >  2:819520/1‑150 (MQ=255)
                                                aGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATg                                                                                                    <  2:1665914/150‑1 (MQ=255)
                                                    gCCACCCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGAc                                                                                                >  2:1523856/1‑150 (MQ=255)
                                                    gCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGATATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATTCCGCCAATGCGAc                                                                                                >  2:1405060/1‑150 (MQ=255)
                                                     ccAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGTGGATGTCGTTAGCCACGGGAGATTAATCTCATAAATAATTCACCTTGTCTCCATAACGTCGACAACGTGAACGAAGATGTCTATTAATCTATACCGCCAATGCGACt                                                                                               >  1:252942/1‑150 (MQ=255)
                                                     ccAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACt                                                                                               <  2:1626130/150‑1 (MQ=255)
                                                       aaCCTTAAACTTCGGGGGAATATTCTATACACTACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCCCGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTa                                                                                             <  2:874532/150‑1 (MQ=255)
                                                        aCCTTAAACTTCGGGGGAATATTCTATACACAACGACTGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTaa                                                                                            <  1:315931/150‑1 (MQ=255)
                                                          cTTAAACTTCGGGGGATTATTCTCTACACTACGACGGGGGATGTCGTTAGCCCCGGGAGATTTATCTCATAAATAATACACGTTGTCTCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAgc                                                                                          <  1:1125992/150‑1 (MQ=255)
                                                             aaaCTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCa                                                                                       <  1:1115182/150‑1 (MQ=255)
                                                                  tCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGTGCAGACCa                                                                                  <  1:1699350/150‑1 (MQ=255)
                                                                      gggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCg                                                                              >  2:1209823/1‑150 (MQ=255)
                                                                        gAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGc                                                                            <  1:1339705/150‑1 (MQ=255)
                                                                        gAATATTCTAAACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGc                                                                            <  1:743674/150‑1 (MQ=255)
                                                                           tatTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCAt                                                                         <  2:207746/150‑1 (MQ=255)
                                                                            atTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGTCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCAtt                                                                        >  1:901105/1‑150 (MQ=255)
                                                                              tCTATACCCAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCg                                                                      <  1:835607/150‑1 (MQ=255)
                                                                              tCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCg                                                                      <  1:64607/150‑1 (MQ=255)
                                                                               cTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGc                                                                     <  1:873809/150‑1 (MQ=255)
                                                                                tataCACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGcc                                                                    >  2:133686/1‑150 (MQ=255)
                                                                                tataCACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGcc                                                                    <  2:963329/150‑1 (MQ=255)
                                                                                tataCACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGcc                                                                    <  2:1542310/150‑1 (MQ=255)
                                                                                 ataCACAACGACGGTGGCTGGCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATACCGGCGACAACGTGAACGACGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCct                                                                   <  1:198196/150‑1 (MQ=255)
                                                                                   acacAACGACGGGGGCTGGCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttc                                                                 <  1:360335/150‑1 (MQ=255)
                                                                                     acaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttctt                                                               <  1:1516677/150‑1 (MQ=255)
                                                                                     acaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATATCATAAATAATTCACCTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttctt                                                               >  1:1337795/1‑150 (MQ=255)
                                                                                      caACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTg                                                              <  2:872898/150‑1 (MQ=255)
                                                                                      caACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTg                                                              >  2:967567/1‑150 (MQ=255)
                                                                                         cgacgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCtt                                                           <  2:237455/150‑1 (MQ=255)
                                                                                            cgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCACCTTCTTGCTTATc                                                        <  2:1091655/150‑1 (MQ=255)
                                                                                            cgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATAAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATc                                                        <  2:1636432/150‑1 (MQ=255)
                                                                                             gggggATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCg                                                       <  2:1216985/150‑1 (MQ=255)
                                                                                              ggggATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGa                                                      <  1:1440699/150‑1 (MQ=255)
                                                                                                 gATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGc                                                   <  2:1220583/150‑1 (MQ=255)
                                                                                                 gATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGc                                                   >  1:201739/1‑150 (MQ=255)
                                                                                                    gTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACACCGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGTCTTACGCGCAGACCAGTCGGCCATTCGCCTTCTGGCTGATCGAAGCTCa                                                >  1:924701/1‑150 (MQ=255)
                                                                                                    gTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCa                                                >  1:168486/1‑150 (MQ=255)
                                                                                                    gTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCa                                                <  2:201739/150‑1 (MQ=255)
                                                                                                      cGTTAGCCACGGGAGATTTATCTCATAAATAATACACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATc                                              >  2:1592769/1‑150 (MQ=255)
                                                                                                        ttAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCa                                            <  1:252989/150‑1 (MQ=255)
                                                                                                         tAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttctt                                                               <  1:591279/130‑1 (MQ=255)
                                                                                                          aGCCACGGGAGATTTATCTCATTAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGa                                          <  2:1141762/150‑1 (MQ=255)
                                                                                                           gCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGAt                                         <  1:441780/150‑1 (MQ=255)
                                                                                                           gCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGAt                                         >  2:455554/1‑150 (MQ=255)
                                                                                                            ccACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATa                                        >  1:1664210/1‑150 (MQ=255)
                                                                                                            ccACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATa                                        <  2:396180/150‑1 (MQ=255)
                                                                                                              actggAGATTTATCTCATAAATAATTCACGTTGTCGCCATTACGGAGACAACGTGAACTAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGc                                      <  1:260311/147‑1 (MQ=255)
                                                                                                                 ggAGATTTATCTCATAAATTATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGt                                   <  2:1341/150‑1 (MQ=255)
                                                                                                                   ttaTTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGAAAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTg                                 <  1:1264214/148‑1 (MQ=255)
                                                                                                                    gaTTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGc                                <  1:1209823/150‑1 (MQ=255)
                                                                                                                       ttATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCccg                             >  1:1442683/1‑150 (MQ=255)
                                                                                                                        tATCTCATAAATAATTCACGTTGTCGCCATAACGGCGAAAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTAGCCTTATTGCTTATCGAAGATCATACAGATAGCCGTTGCccga                            >  2:1392555/1‑150 (MQ=255)
                                                                                                                         ttctcTTATATTATTCAAGTTGTCGCCATACCGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCGGCCTATGCGACATAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCTTTGCccgac                           <  1:116122/149‑1 (MQ=255)
                                                                                                                         aTCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCccgac                           <  1:359442/150‑1 (MQ=255)
                                                                                                                           ctcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCccgaccg                         <  1:486028/150‑1 (MQ=255)
                                                                                                                            tcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCccgaccga                        >  2:577404/1‑150 (MQ=255)
                                                                                                                             cATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCcgaccgac                       >  1:846726/1‑150 (MQ=255)
                                                                                                                             cATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCcgaccgac                       >  2:742570/1‑150 (MQ=255)
                                                                                                                             cATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCcgaccgac                       >  2:1669620/1‑150 (MQ=255)
                                                                                                                               taaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAg                                       <  1:918342/132‑1 (MQ=255)
                                                                                                                               taaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAg                                       >  2:918342/1‑132 (MQ=255)
                                                                                                                               taaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCa                     >  2:380311/1‑150 (MQ=255)
                                                                                                                               taaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCa                     <  2:948987/150‑1 (MQ=255)
                                                                                                                                 aataatTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGCCCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGa                   >  2:867893/1‑150 (MQ=255)
                                                                                                                                  ataatTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGAt                  <  1:327688/150‑1 (MQ=255)
                                                                                                                                  ataatTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGAt                  <  1:1422668/150‑1 (MQ=255)
                                                                                                                                     atTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACGTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCg               >  1:39692/1‑150 (MQ=255)
                                                                                                                                        cACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCtt            >  2:1310971/1‑150 (MQ=255)
                                                                                                                                          cGTTGTCGCCGTTATGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTcc          >  2:1240269/1‑150 (MQ=255)
                                                                                                                                           gTTGTCGCCGTTATGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCg         >  1:959384/1‑150 (MQ=255)
                                                                                                                                            ttGTCGCCGTTATGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGGCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGCCCGACCAGATTCGCTTCCgg        >  1:485080/1‑150 (MQ=255)
                                                                                                                                               tCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGcac     <  1:967567/150‑1 (MQ=255)
                                                                                                                                                 gCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACaa   >  1:663066/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAA  >  NC_000913/2704156‑2704451

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: