Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,866,487 G→C 7.7% intergenic (+22/+72) ygbN → / ← rpoS putative transporter/RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,866,4870GC 7.7% 728.3 / 41.6 208intergenic (+22/+72)ygbN/rpoSputative transporter/RNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base G (97/95);  new base C (11/5);  total (108/100)
Fisher's exact test for biased strand distribution p-value = 1.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

Only 100 of 215 total aligned reads displayed.

GTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCT  >  NC_000913/2866339‑2866635
                                                                                                                                                    |                                                                                                                                                     
gTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATgg                                                                                                                                                       <  1:52149/149‑1 (MQ=255)
 ttCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATggg                                                                                                                                                      >  4:600681/1‑149 (MQ=255)
  tCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGc                                                                                                                                                     <  2:1047090/149‑1 (MQ=255)
    tGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCAt                                                                                                                                                   >  2:466635/1‑149 (MQ=255)
     ggATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATc                                                                                                                                                  <  3:293935/149‑1 (MQ=255)
      gATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCg                                                                                                                                                 <  1:914554/149‑1 (MQ=255)
        ttGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGa                                                                                                                                               <  4:1253349/149‑1 (MQ=255)
        ttGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGa                                                                                                                                               <  4:723270/149‑1 (MQ=255)
             aCCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑ttt                                                                                                                                         >  2:1359973/1‑149 (MQ=255)
              ccAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑tttt                                                                                                                                        >  2:919940/1‑149 (MQ=255)
              ccAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑tttt                                                                                                                                        <  1:816563/149‑1 (MQ=255)
              ccAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑tttt                                                                                                                                        >  4:329610/1‑149 (MQ=255)
              ccAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑tttt                                                                                                                                        >  2:198299/1‑149 (MQ=255)
               cAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTa                                                                                                                                       >  3:867542/1‑149 (MQ=255)
                aaaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTAt                                                                                                                                      >  1:1249701/1‑149 (MQ=255)
                aaaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTAt                                                                                                                                      >  1:1210494/1‑149 (MQ=255)
                 aaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTAtt                                                                                                                                     <  2:143320/149‑1 (MQ=255)
                  atatCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATtg                                                                                                                                    >  3:438327/1‑149 (MQ=255)
                   tatCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATtgt                                                                                                                                   <  2:120767/149‑1 (MQ=255)
                    atCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGCCCATCGGACC‑tttttattaa                                                                                                                                  <  2:783086/149‑7 (MQ=255)
                        ggggTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGcaca                                                                                                                              <  4:223008/149‑1 (MQ=255)
                             tGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGaaaa                                                                                                                         <  2:14950/149‑1 (MQ=255)
                              gTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAg                                                                                                                        >  4:407776/1‑149 (MQ=255)
                              gTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAg                                                                                                                        <  3:45876/149‑1 (MQ=255)
                               tCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAgg                                                                                                                       <  4:867967/149‑1 (MQ=255)
                                cGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGc                                                                                                                      <  2:56062/149‑1 (MQ=255)
                                          cGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCt                                                                                                            <  3:753076/149‑1 (MQ=255)
                                             tCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGaaaa                                                                                                                         <  3:643628/133‑1 (MQ=255)
                                             tCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGaaaa                                                                                                                         >  4:643628/1‑133 (MQ=255)
                                             tCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTga                                                                                                         <  1:222695/149‑1 (MQ=255)
                                              cTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTgag                                                                                                        >  3:386094/1‑149 (MQ=255)
                                               tGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGAc                                                                                                                                              <  3:470929/111‑1 (MQ=255)
                                               tGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGAc                                                                                                                                              >  4:470929/1‑111 (MQ=255)
                                               tGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTgaga                                                                                                       >  4:1266358/1‑149 (MQ=255)
                                               tGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTgaga                                                                                                       <  1:588429/149‑1 (MQ=255)
                                                 aaaaCCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACt                                                                                                     >  1:101564/1‑149 (MQ=255)
                                                 aaaaCCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACt                                                                                                     >  4:253911/1‑149 (MQ=255)
                                                   aaCCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTgg                                                                                                   >  4:1081225/1‑149 (MQ=255)
                                                       tGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGc                                                                                                                      >  2:162518/1‑126 (MQ=255)
                                                       tGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCtt                                                                                               >  1:498358/1‑149 (MQ=255)
                                                        ggACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCttt                                                                                              <  3:927308/149‑1 (MQ=255)
                                                         gACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTc                                                                                             <  3:235232/149‑1 (MQ=255)
                                                           cAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGCGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGTCTCAAGCGAGGCTGGCCttttct                                                                                           <  1:1096884/149‑4 (MQ=255)
                                                            aGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGTTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGa                                                                                          >  4:785147/1‑149 (MQ=255)
                                                              tgtTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCt                                                                                                            <  3:107685/129‑1 (MQ=255)
                                                              tgtTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTc                                                                                             >  4:1272536/1‑144 (MQ=255)
                                                              tgtTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACa                                                                                        <  3:24250/149‑1 (MQ=255)
                                                              tgtTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACa                                                                                        >  2:394733/1‑149 (MQ=255)
                                                                  aaCGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCtt                                                                                                           >  1:1328245/1‑126 (MQ=255)
                                                                  aaCGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATg                                                                                    <  4:439332/149‑1 (MQ=255)
                                                                  aaCGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATg                                                                                    >  3:1399357/1‑149 (MQ=255)
                                                                      aCCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGctta                                                                                <  3:339628/149‑1 (MQ=255)
                                                                       ccATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttac                                                                               >  4:866427/1‑149 (MQ=255)
                                                                        cATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttact                                                                              >  4:177717/1‑149 (MQ=255)
                                                                        cATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttact                                                                              >  3:1308711/1‑149 (MQ=255)
                                                                             tcGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTc                                                                         >  2:659395/1‑149 (MQ=255)
                                                                                   ttACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGaa                                                                   <  1:1172614/149‑1 (MQ=255)
                                                                                     aCCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACa                                                                 >  1:928530/1‑149 (MQ=255)
                                                                                      ccGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAg                                                                <  1:1459404/149‑1 (MQ=255)
                                                                                              ttAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGa                                                        >  1:805312/1‑149 (MQ=255)
                                                                                               tAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGTAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGat                                                       <  2:1150968/149‑1 (MQ=255)
                                                                                                aaTTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGata                                                      <  3:868477/149‑1 (MQ=255)
                                                                                                 aTTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGatat                                                     >  1:951022/1‑149 (MQ=255)
                                                                                                  ttACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATAtt                                                    >  1:790157/1‑149 (MQ=255)
                                                                                                  ttACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATAtt                                                    >  2:1298730/1‑149 (MQ=255)
                                                                                                     ccTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAg                                                 <  4:359115/149‑1 (MQ=255)
                                                                                                      cTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGc                                                >  3:786753/1‑149 (MQ=255)
                                                                                                             gTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTga                                                                                                         >  3:561714/1‑85 (MQ=255)
                                                                                                                tGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCCGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGttt                                      <  1:954444/149‑1 (MQ=255)
                                                                                                                   gCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCGTCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCa                                   <  1:261745/149‑1 (MQ=255)
                                                                                                                   gCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCa                                   <  1:397257/149‑1 (MQ=255)
                                                                                                                   gCGGTAATGTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCa                                   <  3:405828/149‑1 (MQ=255)
                                                                                                                     ggTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAgg                                 >  3:124533/1‑149 (MQ=255)
                                                                                                                      gTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGa                                >  4:1000772/1‑149 (MQ=255)
                                                                                                                      gTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGa                                >  1:175378/1‑149 (MQ=255)
                                                                                                                           ttGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGa                                                                                          >  4:543679/1‑86 (MQ=255)
                                                                                                                           ttGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTcg                           <  2:1249701/149‑1 (MQ=255)
                                                                                                                           ttGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTcg                           >  4:821949/1‑149 (MQ=255)
                                                                                                                               ccTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCa                       >  3:1195399/1‑149 (MQ=255)
                                                                                                                                 ttAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGc                                                <  2:832738/122‑1 (MQ=255)
                                                                                                                                 ttAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaaa                     >  3:750854/1‑149 (MQ=255)
                                                                                                                                 cacaataaaaGGTCCGATGCCCATCGGACCTTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAg                                  <  2:986435/131‑1 (MQ=38)
                                                                                                                                     taaAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGgc                 >  1:674952/1‑149 (MQ=255)
                                                                                                                                     taaAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGgc                 <  1:660449/149‑1 (MQ=255)
                                                                                                                                      aaaaaGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGgcg                >  4:789767/1‑149 (MQ=255)
                                                                                                                                         aaGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg             <  2:674952/149‑1 (MQ=255)
                                                                                                                                          aGGTCCGATGCCCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg            >  4:576546/1‑149 (MQ=255)
                                                                                                                                          aGGTCCGATGCCCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTATGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg            >  4:101822/1‑149 (MQ=255)
                                                                                                                                          aGGTCCGACGCCCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg            >  1:1339144/1‑149 (MQ=255)
                                                                                                                                           ggTCCGATGCCCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc           >  1:1289435/1‑149 (MQ=255)
                                                                                                                                           ggTCCGATGCCCATCGGACC‑TTTTATTGTGCACAGAAAAGGCAAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttactta                                                                            >  2:283693/1‑84 (MQ=38)
                                                                                                                                            gTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc          >  1:557367/1‑149 (MQ=255)
                                                                                                                                               cGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTcg                                                                        >  4:1418274/1‑84 (MQ=255)
                                                                                                                                               cGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCtt        >  2:553601/1‑148 (MQ=255)
                                                                                                                                                gATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt         >  3:1209749/1‑146 (MQ=255)
                                                                                                                                                gATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt         <  4:1209749/146‑1 (MQ=255)
                                                                                                                                                 aTGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCaa     <  1:324095/149‑1 (MQ=255)
                                                                                                                                                 aTGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCaa     <  3:177717/149‑1 (MQ=255)
                                                                                                                                                  tGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc    <  3:1266358/149‑1 (MQ=255)
                                                                                                                                                  tGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc    <  2:805312/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCT  >  NC_000913/2866339‑2866635

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