Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,866,488 G→C 7.7% intergenic (+23/+71) 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,4880GC 7.7% 726.8 / 42.5 208intergenic (+23/+71)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 = 1.24e-01

Only 100 of 212 total aligned reads displayed.

TTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCT  >  NC_000913/2866340‑2866635
                                                                                                                                                    |                                                                                                                                                    
 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)
         gtgACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGAcc                                                                                                                                             <  4:1367738/149‑1 (MQ=255)
           gACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCtt                                                                                                                                          <  1:1374495/149‑1 (MQ=255)
            aCCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCttt                                                                                                                                         >  2:1359973/1‑149 (MQ=255)
             ccAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCtttt                                                                                                                                        <  1:816563/149‑1 (MQ=255)
             ccAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCtttt                                                                                                                                        >  2:198299/1‑149 (MQ=255)
             ccAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCtttt                                                                                                                                        >  1:1335057/1‑149 (MQ=255)
              cAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTa                                                                                                                                       >  3:867542/1‑149 (MQ=255)
               aaaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTAt                                                                                                                                      >  1:1249701/1‑149 (MQ=255)
               aaaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTAt                                                                                                                                      >  1:1210494/1‑149 (MQ=255)
                aaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTa                                                                                                                                       >  4:330097/1‑147 (MQ=255)
                aaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTa                                                                                                                                       <  3:330097/147‑1 (MQ=255)
                aaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTAtt                                                                                                                                     <  2:143320/149‑1 (MQ=255)
                 atatCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATtg                                                                                                                                    >  3:438327/1‑149 (MQ=255)
                  tatCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATtgt                                                                                                                                   <  2:120767/149‑1 (MQ=255)
                   atCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGCCCATCGGACCtttttattaa                                                                                                                                  <  2:783086/149‑7 (MQ=255)
                       ggggTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGcaca                                                                                                                              <  4:223008/149‑1 (MQ=255)
                            tGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGaaaa                                                                                                                         <  2:14950/149‑1 (MQ=255)
                             gTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAg                                                                                                                        <  1:665140/149‑1 (MQ=255)
                             gTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAg                                                                                                                        <  3:45876/149‑1 (MQ=255)
                             gTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAg                                                                                                                        >  4:407776/1‑149 (MQ=255)
                              tCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAgg                                                                                                                       <  4:867967/149‑1 (MQ=255)
                                       gACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCt                                                                                                                                           >  1:512610/1‑120 (MQ=255)
                                       gACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCt                                                                                                                                           <  2:512610/120‑1 (MQ=255)
                                         cGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCt                                                                                                            <  3:753076/149‑1 (MQ=255)
                                            tCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGaaaa                                                                                                                         <  3:643628/133‑1 (MQ=255)
                                            tCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTga                                                                                                         <  1:222695/149‑1 (MQ=255)
                                             cTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTgag                                                                                                        >  3:386094/1‑149 (MQ=255)
                                              tGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGAc                                                                                                                                              >  4:470929/1‑111 (MQ=255)
                                              tGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGAc                                                                                                                                              <  3:470929/111‑1 (MQ=255)
                                              tGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTgaga                                                                                                       <  4:683210/149‑1 (MQ=255)
                                              tGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTgaga                                                                                                       >  4:1266358/1‑149 (MQ=255)
                                              tGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTgaga                                                                                                       <  3:960622/149‑1 (MQ=255)
                                                aaaaCCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACt                                                                                                     >  1:101564/1‑149 (MQ=255)
                                                  aaCCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTgg                                                                                                   >  4:1081225/1‑149 (MQ=255)
                                                      tGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGc                                                                                                                      >  2:162518/1‑126 (MQ=255)
                                                       ggACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCttt                                                                                              <  3:927308/149‑1 (MQ=255)
                                                         aCAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCt                                                                                            <  1:466635/149‑1 (MQ=255)
                                                          cAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGCGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGTCTCAAGCGAGGCTGGCCttttct                                                                                           <  1:1096884/149‑4 (MQ=255)
                                                           aGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGTTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGa                                                                                          >  4:785147/1‑149 (MQ=255)
                                                           aGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGa                                                                                          <  4:1272127/149‑1 (MQ=255)
                                                             tgtTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTc                                                                                             <  3:1272536/144‑1 (MQ=255)
                                                             tgtTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACa                                                                                        >  2:394733/1‑149 (MQ=255)
                                                               ttAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGa                                                                                      >  2:946660/1‑149 (MQ=255)
                                                                 aaCGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCtt                                                                                                           >  1:1328245/1‑126 (MQ=255)
                                                                 aaCGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATg                                                                                    <  4:439332/149‑1 (MQ=255)
                                                                      ccATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttac                                                                               >  4:866427/1‑149 (MQ=255)
                                                                       cATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttact                                                                              >  2:892529/1‑149 (MQ=255)
                                                                            tcGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTc                                                                         >  2:659395/1‑149 (MQ=255)
                                                                                  ttACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGaa                                                                   <  1:1172614/149‑1 (MQ=255)
                                                                                    aCCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACa                                                                 >  1:928530/1‑149 (MQ=255)
                                                                                      cGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAgc                                                               <  4:383877/149‑1 (MQ=255)
                                                                                             ttAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGa                                                        >  1:805312/1‑149 (MQ=255)
                                                                                              tAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGTAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGat                                                       <  2:1150968/149‑1 (MQ=255)
                                                                                               aaTTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGata                                                      <  3:868477/149‑1 (MQ=255)
                                                                                                 ttACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATAtt                                                    >  2:1298730/1‑149 (MQ=255)
                                                                                                 ttACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATAtt                                                    >  1:790157/1‑149 (MQ=255)
                                                                                                     cTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGc                                                >  3:786753/1‑149 (MQ=255)
                                                                                                     cTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGc                                                >  4:1424483/1‑149 (MQ=255)
                                                                                                            gTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTga                                                                                                         >  3:561714/1‑85 (MQ=255)
                                                                                                                  gCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCGTCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCa                                   <  1:261745/149‑1 (MQ=255)
                                                                                                                  gCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCa                                   <  1:397257/149‑1 (MQ=255)
                                                                                                                  gCGGTAATGTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCa                                   <  3:405828/149‑1 (MQ=255)
                                                                                                                    ggTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAgg                                 >  3:124533/1‑149 (MQ=255)
                                                                                                                     gTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGa                                >  1:175378/1‑149 (MQ=255)
                                                                                                                          ttGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGa                                                                                          >  4:543679/1‑86 (MQ=255)
                                                                                                                          ttGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTcg                           >  4:821949/1‑149 (MQ=255)
                                                                                                                          ttGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTcg                           <  2:1249701/149‑1 (MQ=255)
                                                                                                                            gACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTcgcg                         >  4:259647/1‑149 (MQ=255)
                                                                                                                              ccTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCa                       >  3:1195399/1‑149 (MQ=255)
                                                                                                                                ttAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGc                                                <  2:832738/122‑1 (MQ=255)
                                                                                                                                ttAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaaa                     >  3:750854/1‑149 (MQ=255)
                                                                                                                                ttAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaaa                     >  1:95189/1‑149 (MQ=255)
                                                                                                                                cacaataaaaGGTCCGATGCCCATCGGACCTTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAg                                  >  1:986435/7‑137 (MQ=38)
                                                                                                                                    taaAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGgc                 >  3:1336153/1‑149 (MQ=255)
                                                                                                                                     aaaaaGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGgcg                >  4:789767/1‑149 (MQ=255)
                                                                                                                                        aaGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg             <  2:674952/149‑1 (MQ=255)
                                                                                                                                         aGGTCCGATGCCCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg            >  4:637491/1‑149 (MQ=255)
                                                                                                                                         aGGTCCGATGCCCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg            >  4:576546/1‑149 (MQ=255)
                                                                                                                                         aGGTCCGATGCCCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTATGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg            >  4:101822/1‑149 (MQ=255)
                                                                                                                                          ggTCCGATGCCCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc           <  3:101822/149‑1 (MQ=255)
                                                                                                                                          ggTCCGATGCCCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc           >  3:209723/1‑149 (MQ=255)
                                                                                                                                          ggTCCGATGCCCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc           >  1:1289435/1‑149 (MQ=255)
                                                                                                                                          ggTCCGATGCCCATCGGACCTTTTATTGTGCACAGAAAAGGCAAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttactta                                                                            >  2:283693/1‑84 (MQ=38)
                                                                                                                                           gTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc          >  1:557367/1‑149 (MQ=255)
                                                                                                                                           gTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc          >  1:698470/1‑149 (MQ=255)
                                                                                                                                              cGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTcg                                                                        <  3:1418274/84‑1 (MQ=255)
                                                                                                                                              cGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTc       <  2:559618/149‑1 (MQ=255)
                                                                                                                                               gATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt         <  4:1209749/146‑1 (MQ=255)
                                                                                                                                               gATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt         >  3:1209749/1‑146 (MQ=255)
                                                                                                                                                aTGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCaa     <  1:324095/149‑1 (MQ=255)
                                                                                                                                                 tGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc    <  3:1266358/149‑1 (MQ=255)
                                                                                                                                                 tGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc    <  2:805312/149‑1 (MQ=255)
                                                                                                                                                  gggCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc   >  1:642604/1‑149 (MQ=255)
                                                                                                                                                   ggCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt  <  3:240698/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCT  >  NC_000913/2866340‑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: