Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,376,265 G→C 6.7% intergenic (+29/+14) dcuD → / ← sspB putative transporter/ClpXP protease specificity enhancing factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,376,2650GC 6.7% 217.6 / 6.0 83intergenic (+29/+14)dcuD/sspBputative transporter/ClpXP protease specificity enhancing factor
Reads supporting (aligned to +/- strand):  ref base G (42/35);  new base C (3/2);  total (46/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

TGGTCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGC  >  NC_000913/3376118‑3376413
                                                                                                                                                   |                                                                                                                                                    
tggtCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTgg                                                                                                                                                     <  1:970892/149‑1 (MQ=255)
tggtCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGGGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTgg                                                                                                                                                     <  1:668724/149‑1 (MQ=255)
  gtCGTTACGTCGGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGc                                                                                                                                                   >  2:1002959/1‑149 (MQ=255)
  gtCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCTCCTGTGc                                                                                                                                                   >  1:990432/1‑149 (MQ=255)
   tCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGcc                                                                                                                                                  >  1:864607/1‑149 (MQ=255)
     gTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTg                                                                                                                                                >  2:797331/1‑149 (MQ=255)
      ttACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACCATCACGCTATTTTATTAAGTCATTAAAAAGACTAAACAGGCCGCCTGGGCCTGt                                                                                                                                               >  1:851863/1‑149 (MQ=255)
        aCGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGttt                                                                                                                                             <  2:614124/149‑1 (MQ=255)
        aCGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGttt                                                                                                                                             >  1:301999/1‑149 (MQ=255)
         cGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCCAGGCGGCCTGtttt                                                                                                                                            <  2:177255/149‑1 (MQ=255)
         cGTCAGGAATTGCAGGCGTTTCGCCTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTg                                                                                                                                           <  2:121258/149‑1 (MQ=255)
           tCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGt                                                                                                                                          <  2:947082/149‑1 (MQ=255)
           tCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGGCGGAGTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTCAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGt                                                                                                                                          <  1:113692/149‑1 (MQ=255)
            cAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGCAGCGGACAGCGATCCCCATGGCAGTAGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAATAAGACAAAACAGGCCGCCTGGGCCTGTTTTGtt                                                                                                                                         >  1:656036/1‑148 (MQ=255)
            cAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCCAGGCGGCCTGTTTTGtc                                                                                                                                         <  1:1397/149‑2 (MQ=255)
            cAGGAATTGCAGGCGGTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGCTCCCCATGGCAGTCGGTTTCGTGGTGTTTATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGACGCCTGGGCCTGTTTTGtt                                                                                                                                         >  1:788462/1‑148 (MQ=255)
             acgAATTGCACGCGTTTCGCCTTTTGCGGTTGTGAAGCGGACATCGATCCCCATGGCAGTCGTTTTCTTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACACGCCGCCTGGGCCTGTTTTGTAt                                                                                                                                        <  1:1067623/147‑1 (MQ=255)
              ggAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTAtt                                                                                                                                       >  2:264355/1‑149 (MQ=255)
              ggAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTAtt                                                                                                                                       >  1:470012/1‑149 (MQ=255)
              ggAATTGCAGGCGTGTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTAtt                                                                                                                                       <  1:1032981/149‑1 (MQ=255)
               gAATCGCAGGCGCTTCGCCTTTTGCGGTGGTCAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGGGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTa                                                                                                                                      <  2:774080/149‑1 (MQ=255)
                   tGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTc                                                                                                                                  >  1:173434/1‑149 (MQ=255)
                      aGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTAGTGGTTAATATGATTTCCACAATCACGCTATTTAATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTcaca                                                                                                                               <  2:374724/149‑1 (MQ=255)
                           tttCGCCTTTTGCGGTGGTGAAGCGGACACCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTTTTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAAcgcg                                                                                                                          <  2:716672/149‑1 (MQ=255)
                             tCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTa                                                                                                                        >  2:574225/1‑149 (MQ=255)
                                 cTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTGCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGc                                                                                                                    <  2:13767/149‑1 (MQ=255)
                                  ttttGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGcc                                                                                                                   >  2:837549/1‑149 (MQ=255)
                                      gCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCTTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTATAAAGACAAAACAGGCCGCCTAGGCCTATATTGTATTACTTCAAAACGCGAAATCCCGTTc                                                                                                               >  1:818675/1‑149 (MQ=255)
                                      gCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAAAGCGTAATGCCGGTc                                                                                                               >  1:730998/1‑149 (MQ=255)
                                      gCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGGGGTGAATATGATTGCCACAATCACGCTATTTTATTCAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTc                                                                                                               <  2:1054137/149‑1 (MQ=255)
                                        ggtggtGAAGCGGACAGAGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAAAGCCGGTCgc                                                                                                             >  2:186213/1‑148 (MQ=255)
                                          tggtgAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGacc                                                                                                           <  2:505170/149‑1 (MQ=255)
                                             tgAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGaccacc                                                                                                        >  1:760097/1‑149 (MQ=255)
                                                    gACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCTTAATGCAGGTCGACACCCGCgtggt                                                                                                 >  2:1059850/1‑149 (MQ=255)
                                                       aGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGc                                                                                              >  2:227645/1‑149 (MQ=255)
                                                           aTCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGACGGTCGACCACCCCGTGGTGGCTGCg                                                                                          >  1:110707/1‑149 (MQ=255)
                                                            tCCCCATGGCCGTCGGTTTCGTGGTGAATATGATTGCCACATTCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCgg                                                                                         <  1:168606/149‑1 (MQ=255)
                                                              cccATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAg                                                                                       <  1:264355/149‑1 (MQ=255)
                                                              cccAGGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAACAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAg                                                                                       <  2:517502/149‑1 (MQ=255)
                                                               ccATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAgg                                                                                      >  2:533374/1‑149 (MQ=255)
                                                               ccATGGCAGTCGGTTTCGTGGGTAATATGATTGCCACAATCACGCTATTTTATTAAGGCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAgg                                                                                      <  2:1141823/149‑1 (MQ=255)
                                                                cATGGCAGTCGGTTTCGTGGTGAATTTGATTGCCACTATGACGCTATTTTATTAAGTCTTTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTAGTTCACAACGCGTATTTCCGGTCGACCACCGAGTTGTGGATGCGGAGGt                                                                                     >  2:816206/1‑149 (MQ=255)
                                                                   ggCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCCCCCCCCGGTTGTTGCATGCGGAGGTTCa                                                                                  >  1:287327/1‑149 (MQ=255)
                                                                     cAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAACCAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATc                                                                                <  2:128537/149‑1 (MQ=255)
                                                                     cAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATc                                                                                <  2:409368/149‑1 (MQ=255)
                                                                     cAGTCGGATTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATc                                                                                <  2:124263/149‑1 (MQ=255)
                                                                           gTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGCCCTGTTATGTATTACTTCAAAACCCGTAATGCAGGTCGACCCCCGCGTGGTGGCTGCGGAGGTTCATCGTCAgg                                                                          >  1:755977/1‑149 (MQ=255)
                                                                            tatCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGa                                                                         <  1:823926/147‑1 (MQ=255)
                                                                                  ggtgAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAgtg                                                                   >  1:409627/1‑149 (MQ=255)
                                                                                   gttaaTATGATTGCCACAATCACGCTATATTATTAAGTCATTAAATAGACAAAAGAGGGCGCCTGGGCCTGTTTTGTATTACTTCACACCGCGTAATGACGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAgtgt                                                                  <  1:475014/146‑1 (MQ=255)
                                                                                    tgAATTTGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACATAACAGGGCGCCTGGGGCTGTTTTGTATTACTTCACATCGCGTATTGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGtc                                                                 <  1:597786/149‑1 (MQ=255)
                                                                                     gAATATGATTGCCAAAATCACGCTAAATTATTAAGTCATAAAAAAGACAAAACAGGCCGCCTGGGCCTGATTTGTATAACTTCAAAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCAACCTCAGGATGAATGtca                                                                >  2:747487/1‑149 (MQ=255)
                                                                                        tatGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGTGGAGTGtctgca                                                             >  1:798109/1‑145 (MQ=255)
                                                                                        tatGATTGCCACAATAACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGtcatca                                                             <  2:851863/149‑1 (MQ=255)
                                                                                           gATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCg                                                          >  2:264437/1‑149 (MQ=255)
                                                                                                ccACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATc                                                     >  1:454824/1‑149 (MQ=255)
                                                                                                 cacaATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCt                                                    <  2:280853/149‑1 (MQ=255)
                                                                                                  acaATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTg                                                   >  2:1026085/1‑149 (MQ=255)
                                                                                                   caATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTgg                                                  <  1:797331/149‑1 (MQ=255)
                                                                                                          gCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGTTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCg                                           <  2:1064227/149‑1 (MQ=255)
                                                                                                           cTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCAGGGCCAGTTTTGTATTACTTCACAACTCGTAATGCCGGTCGACCACCCCGTGGTGGCTGCGGAGCTTCATCGTCAGGATGAGTGTAATCATCCTGATCTGGCTTGACGc                                          >  1:594329/1‑149 (MQ=255)
                                                                                                               tttATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATc                                      >  1:1012487/1‑149 (MQ=255)
                                                                                                                ttattaAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCa                                     >  2:1004202/1‑149 (MQ=255)
                                                                                                                 tattaAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCaa                                    <  1:1002959/149‑1 (MQ=255)
                                                                                                                 tattaAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCaa                                    >  2:1153547/1‑149 (MQ=255)
                                                                                                                  attaAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAAt                                   <  2:301999/149‑1 (MQ=255)
                                                                                                                  attaAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAAt                                   >  2:456168/1‑149 (MQ=255)
                                                                                                                       gTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCg                              >  1:461575/1‑149 (MQ=255)
                                                                                                                        tCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGa                             <  2:1074753/149‑1 (MQ=255)
                                                                                                                           ttAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGGGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACAt                          >  2:451977/1‑149 (MQ=255)
                                                                                                                              gtaataCAAAACAGGCCCAGGCGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAAc                       >  2:354194/6‑149 (MQ=255)
                                                                                                                               aaaGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCTGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACg                      >  1:454034/1‑149 (MQ=255)
                                                                                                                                 ataCAAAACAGGCCCAGTCGGCCTGTTTTGTATTACTTCACAACGCTTAATGCCGGACGACCACCGCGTGGTGGCTGCGTAGGATCATCTACATGATGAGTGACATAATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGt                    >  2:680120/3‑149 (MQ=255)
                                                                                                                                  taCAAAACAGGCCCAGGCGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGtt                   >  1:345934/2‑149 (MQ=255)
                                                                                                                                   aCAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAAGCGACATAACGGgtt                  >  1:860574/1‑149 (MQ=255)
                                                                                                                                     aaaaCAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGCCATAACGGTTTCg                >  1:988174/1‑149 (MQ=255)
                                                                                                                                     aaaaCAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCg                >  2:817232/1‑149 (MQ=255)
                                                                                                                                      aaaCAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGt               <  2:749121/149‑1 (MQ=255)
                                                                                                                                       aaCAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGtt              >  2:1130239/1‑149 (MQ=255)
                                                                                                                                         cAGGCCGCCTGGGCCTGTTTTGTATTACTTCACCACGCGTAATGCCGGGCGTCCCCCGCGGGGTGGCTGCGGAGGTTCATCGTCAGGAGGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCATTAACCGACATAACGGTTTCGTTGt            <  2:906493/149‑1 (MQ=255)
                                                                                                                                         cAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGt            <  1:533374/149‑1 (MQ=255)
                                                                                                                                          aGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTc           >  2:56399/1‑149 (MQ=255)
                                                                                                                                             ccgccTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGc        <  2:470012/149‑1 (MQ=255)
                                                                                                                                                 cTGTGCCTGTTTTGTATTACTTCACAACGCTTAATGCGTGTCGACCACCGCGGGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGAt    <  1:994638/149‑1 (MQ=255)
                                                                                                                                                  tGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATg   <  2:903817/149‑1 (MQ=255)
                                                                                                                                                   gggCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGc  <  1:574225/149‑1 (MQ=255)
                                                                                                                                                   gggCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGc  <  2:760097/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGGTCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGC  >  NC_000913/3376118‑3376413

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: