Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,704,304 A→G 12.2% 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,3040AG12.2% 290.3 / 28.8 115intergenic (‑241/+31)rnc/lepBRNase III/leader peptidase (signal peptidase I)
Reads supporting (aligned to +/- strand):  ref base A (48/53);  new base G (5/9);  total (53/62)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-02

Only 100 of 116 total aligned reads displayed.

TTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAG  >  NC_000913/2704155‑2704452
                                                                                                                                                      |                                                                                                                                                    
ttACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCa                                                                                                                                                      >  1:575561/1‑150 (MQ=255)
  aCCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATa                                                                                                                                                    <  2:1658829/150‑1 (MQ=255)
     ggTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACg                                                                                                                                                 >  1:1157609/1‑150 (MQ=255)
     ggTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACg                                                                                                                                                 >  1:234836/1‑150 (MQ=255)
     ggTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACg                                                                                                                                                 >  1:1732562/1‑150 (MQ=255)
     ggTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACg                                                                                                                                                 >  2:376166/1‑150 (MQ=255)
     ggTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCAAGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACg                                                                                                                                                 >  1:311244/1‑150 (MQ=255)
        atatGGTGACCAACGCTGTTTCGTGTGCCGTGGCAACATGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACTGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCg                                                                                                                                              <  1:98957/150‑1 (MQ=255)
            ggggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACaa                                                                                                                                          >  2:583993/1‑150 (MQ=255)
             gggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCGACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAAc                                                                                                                                         >  1:104101/1‑150 (MQ=255)
             gggACCAACGCTGTTTCGTGTGCCGTGGCAACCGGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAAc                                                                                                                                         <  1:1207801/150‑1 (MQ=255)
              ggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACg                                                                                                                                        <  1:612404/150‑1 (MQ=255)
              ggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACg                                                                                                                                        >  2:1712579/1‑150 (MQ=255)
               taCCAACGCAGTTTCGTGTGCCGTGGCAACCTGGACGTGCCAACCTTAAACATCGGGGGAATATTCAATACACAACGACGGGGGATGGCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGt                                                                                                                                       <  1:1563009/149‑1 (MQ=255)
                aCCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACGTg                                                                                                                                      <  1:1391998/150‑1 (MQ=255)
                aCCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACGTg                                                                                                                                      <  2:520179/150‑1 (MQ=255)
                 ccAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGa                                                                                                                                     >  1:1161558/1‑150 (MQ=255)
                   aaCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTTCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACGTGaat                                                                                                                                   <  2:352862/150‑2 (MQ=255)
                   aaCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCTGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAAc                                                                                                                                   >  2:733568/1‑150 (MQ=255)
                   aaCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGaaa                                                                                                                                   >  2:1542717/1‑149 (MQ=255)
                   aaCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAAc                                                                                                                                   >  2:1617071/1‑150 (MQ=255)
                   aaCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAAc                                                                                                                                   <  1:937970/150‑1 (MQ=255)
                   aaCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTAAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAAc                                                                                                                                   <  1:934764/150‑1 (MQ=255)
                    aCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTACACTTCGGGGGAATATTCTATACACCACGACGGGGGATGTCGTTAGCCACGGGAGATTAATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACGTGaatt                                                                                                                                  <  1:58163/150‑3 (MQ=255)
                    aCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACg                                                                                                                                  <  1:268206/150‑1 (MQ=255)
                      gCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGaa                                                                                                                                >  1:579404/1‑150 (MQ=255)
                       cTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAg                                                                                                                               >  1:67061/1‑150 (MQ=255)
                          tttCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATg                                                                                                                            <  2:357555/150‑1 (MQ=255)
                            tCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGc                                                                                                                          <  1:1719697/150‑1 (MQ=255)
                             cGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGTGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCt                                                                                                                         <  1:106083/150‑1 (MQ=255)
                                 tgCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTa                                                                                                                     >  2:1424874/1‑150 (MQ=255)
                                     gTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATgg                                                                                                                 <  2:1732562/150‑1 (MQ=255)
                                     gTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATgg                                                                                                                 >  1:1536795/1‑150 (MQ=255)
                                       ggCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCTACAACGTGAACGAAGATGGCTATTAATGGAt                                                                                                               <  1:1039352/150‑1 (MQ=255)
                                       ggCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGAt                                                                                                               >  1:890037/1‑150 (MQ=255)
                                        gCAAGCTGGACGTGCCAACCTTAAACTTCGGGGGAATATTCTAAACACAACGACGGTGGATGTCGTTAGCCAAGGGAGATTTATCTGATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATg                                                                                                              <  2:121633/150‑1 (MQ=255)
                                            ccTGGAGGTGCCAACCTTAAACTACGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgcccc                                                                                                          >  2:675891/1‑148 (MQ=255)
                                                gAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAAt                                                                                                      <  2:631183/150‑1 (MQ=255)
                                                  ggTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGc                                                                                                    >  2:1304322/1‑150 (MQ=255)
                                                  ggTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAAAGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGc                                                                                                    >  1:339918/1‑150 (MQ=255)
                                                     gCCTACCTTAAACTTCGGGGGAATATTCTATACACAACGTCGGGTGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGTCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGAc                                                                                                 <  2:858378/150‑1 (MQ=255)
                                                      ccAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACt                                                                                                <  1:797481/150‑1 (MQ=255)
                                                        aaCCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTa                                                                                              >  2:1399988/1‑150 (MQ=255)
                                                          ccTTAAACTTCGGGGTAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAg                                                                                            <  1:1243968/150‑1 (MQ=255)
                                                           cTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAgc                                                                                           >  1:1274955/1‑150 (MQ=255)
                                                             tAAACTTCGGTGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATATATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAgcgc                                                                                         <  1:214497/150‑1 (MQ=255)
                                                              aaaCTTCGGGGGAATATTCTATACACAACGACTGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCa                                                                                        <  2:1153537/150‑1 (MQ=255)
                                                               aaCTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAg                                                                                       <  2:247341/150‑1 (MQ=255)
                                                               aaCTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAg                                                                                       <  2:1614997/150‑1 (MQ=255)
                                                               aaCTTCGGCGGAATATTCTATACACCACGACGGGGGAGGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGAAGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAg                                                                                       <  2:193980/150‑1 (MQ=255)
                                                                aCTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGa                                                                                      >  2:574949/1‑150 (MQ=255)
                                                                 cTTCGGGGGAATATTCTATACACAACGACGGTGGATGTCGTTAGCCACGGTAGATTTATATCATAAATAATTCACGTAGTCGCAATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGAc                                                                                     <  1:1561559/150‑1 (MQ=255)
                                                                  ttCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGAcc                                                                                    <  2:824623/150‑1 (MQ=255)
                                                                     gggggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGt                                                                                 <  1:1001315/150‑1 (MQ=255)
                                                                      ggggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTc                                                                                >  2:1291059/1‑150 (MQ=255)
                                                                        ggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCgg                                                                              <  1:702437/150‑1 (MQ=255)
                                                                         gAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGc                                                                             <  1:583993/150‑1 (MQ=255)
                                                                         gAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGc                                                                             <  2:1447241/150‑1 (MQ=255)
                                                                          aaTATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGcc                                                                            <  1:954766/150‑1 (MQ=255)
                                                                            tatTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCAt                                                                          <  2:1210306/150‑1 (MQ=255)
                                                                             atTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCAtt                                                                         >  2:1482115/1‑150 (MQ=255)
                                                                              ttCTATACACAACGACGGGGGATGTCGTTAGCCTCGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATTGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTc                                                                        <  1:423720/150‑1 (MQ=255)
                                                                                 tataCACCACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGcc                                                                     >  1:92666/1‑150 (MQ=255)
                                                                                 tataCACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGCTGGCTATTAATTGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATACGcc                                                                     >  1:279882/1‑150 (MQ=255)
                                                                                  ataCACAACGACGGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGcc                                                                     <  2:736818/150‑1 (MQ=255)
                                                                                  ataCACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCct                                                                    >  1:1691228/1‑150 (MQ=255)
                                                                                       caACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTg                                                               <  2:575561/150‑1 (MQ=255)
                                                                                        aaCGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGc                                                              <  2:136435/150‑1 (MQ=255)
                                                                                         acgacgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTACCCATATCGTCGACAACGTGAACGAAGATGGCTATTAATCGATGCCGCCAATGAGACTTAAGCGCAGACCAGTCGGCCATTCTCCTTCTTGCt                                                             <  2:302648/150‑1 (MQ=255)
                                                                                               gggggATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCg                                                        <  1:134154/150‑1 (MQ=255)
                                                                                                ggggATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGGTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGa                                                       <  1:99697/150‑1 (MQ=255)
                                                                                                ggggAGGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGGTGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGa                                                       <  2:1282616/150‑1 (MQ=255)
                                                                                                 gggATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGTCATTCGCCTTCTTGCTTATCGaa                                                      >  1:531275/1‑150 (MQ=255)
                                                                                                   gATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGc                                                    <  2:189926/150‑1 (MQ=255)
                                                                                                      gTCGTAAGCCACGGGAGATTAATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCa                                                 <  1:276624/150‑1 (MQ=255)
                                                                                                         gTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAAAGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATcc                                              <  1:36499/150‑1 (MQ=255)
                                                                                                            aGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGa                                           <  2:1157609/150‑1 (MQ=255)
                                                                                                               cACGGGTCATTTATCTCTTAAATAATTAACGTTGTCCCCATAACGGCTAAAACTTTAACGAAGATGCCAATTAATGGATGCCGCCAATGCTACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAg                                        <  1:1617071/150‑1 (MQ=255)
                                                                                                               cACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAg                                        <  1:733568/150‑1 (MQ=255)
                                                                                                                aCGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCTGCCAGAAAGc                                       >  2:1143671/1‑150 (MQ=255)
                                                                                                                    gagaTTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGtt                                   >  2:1482275/1‑150 (MQ=255)
                                                                                                                          tATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTAAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCccga                             <  2:451666/150‑1 (MQ=255)
                                                                                                                              tcatcaATAATTCACTGTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCCGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCccgaccga                         <  2:278005/150‑1 (MQ=255)
                                                                                                                              tcATAAATAATTCACGTTGTCGCCATTACGGCGGGCACGTGAACGAAGATGGCTATTATTGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTGTCGATGCTCATCCAGATAGCCGTTGCccgaccga                         >  1:1202005/1‑150 (MQ=255)
                                                                                                                               cATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCcgaccgac                        >  2:1463898/1‑150 (MQ=255)
                                                                                                                                 taaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGATGCCCGACCGACCa                      >  2:1208862/1‑150 (MQ=255)
                                                                                                                                  aaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAg                     <  2:1176130/150‑1 (MQ=255)
                                                                                                                                   aataatTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGa                    >  2:732542/1‑150 (MQ=255)
                                                                                                                                    ataatTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGAt                   <  2:1161558/150‑1 (MQ=255)
                                                                                                                                       gttCACGTTGTCGCCGTTATGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCACATAGCCGTTGCCCGACCGACCAGATTCg                >  1:966439/2‑150 (MQ=255)
                                                                                                                                       atTCACGTTGTCGCCATAACGGCGACAACGTGCACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCGATTCGCCTTCGTGCTTATCGAAGCTCCTCCAGATAGCCGTTGCCCGACCGACCAGAGTCg                >  1:1565337/1‑150 (MQ=255)
                                                                                                                                       atTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGGTGCCCGACCGACCAGATTCg                <  1:147339/150‑1 (MQ=255)
                                                                                                                                           aCGTTGTCGCCGTTATGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTc            >  1:1688180/1‑150 (MQ=255)
                                                                                                                                           aCGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGACGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTc            >  1:1051624/1‑150 (MQ=255)
                                                                                                                                            cGTTGTCGCCGTTATGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGCCCAGATTCGCTTcc           >  2:1429725/1‑150 (MQ=255)
                                                                                                                                            cGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTcc           <  2:1100222/150‑1 (MQ=255)
                                                                                                                                             gTTGTCGCCGTTATGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGCCCAGATTCGCTTCCg          >  2:455353/1‑150 (MQ=255)
                                                                                                                                               tGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGc        <  2:234836/150‑1 (MQ=255)
                                                                                                                                                    ccATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGCCCGACCAGATTCGCTTCCGGCACaaa   >  2:1164923/1‑150 (MQ=255)
                                                                                                                                                     cATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATTCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAg  >  2:827381/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAG  >  NC_000913/2704155‑2704452

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: