Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,329,709 C→G 8.1% G30A (GGT→GCT) ‡ yfaQ ← tandem DUF2300 domain‑containing protein YfaQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,329,7090CG 8.1% 321.5 / 13.6 97G30A (GGT→GCT) ‡yfaQtandem DUF2300 domain‑containing protein YfaQ
Reads supporting (aligned to +/- strand):  ref base C (48/41);  new base G (5/3);  total (53/44)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

GCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCA  >  NC_000913/2329563‑2329858
                                                                                                                                                  |                                                                                                                                                     
gctgctCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACcgc                                                                                                                                                    <  1:8691/150‑1 (MQ=255)
gctgctCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACcgc                                                                                                                                                    <  1:1215871/150‑1 (MQ=255)
   gctCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCgca                                                                                                                                                 <  1:86105/150‑1 (MQ=255)
   gctCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCgca                                                                                                                                                 <  1:933074/150‑1 (MQ=255)
     tCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCgcagc                                                                                                                                               <  1:735068/150‑1 (MQ=255)
      ccGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCgcagca                                                                                                                                              <  1:1081611/150‑1 (MQ=255)
      ccGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCgcagca                                                                                                                                              <  2:637497/150‑1 (MQ=255)
        ggTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCAcc                                                                                                                                            >  1:454138/1‑150 (MQ=255)
        ggTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCAcc                                                                                                                                            >  2:1493933/1‑150 (MQ=255)
         gTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCGGTGCGCAGCAcga                                                                                                                                           <  2:1108315/150‑3 (MQ=255)
            gATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCGGTGCGCAGCAcgatca                                                                                                                                        <  2:1165214/150‑6 (MQ=255)
            gATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGtt                                                                                                                                        >  1:449815/1‑150 (MQ=255)
             aTGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTg                                                                                                                                       <  1:1666271/150‑1 (MQ=255)
               gggTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGt                                                                                                                                         >  2:888923/1‑146 (MQ=255)
               gggTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGt                                                                                                                                         <  1:888923/146‑1 (MQ=255)
               gggTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTa                                                                                                                                     >  2:1290406/1‑150 (MQ=255)
                ggTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTaa                                                                                                                                    >  2:1316743/1‑150 (MQ=255)
                ggTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTaa                                                                                                                                    >  1:270369/1‑150 (MQ=255)
                 gTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAAc                                                                                                                                   <  1:1622527/150‑1 (MQ=255)
                     cctccAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCg                                                                                                                               >  2:1104290/1‑150 (MQ=255)
                      ctccAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGt                                                                                                                              >  2:547470/1‑150 (MQ=255)
                           gccagGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTctcct                                                                                                                         <  2:1026784/150‑1 (MQ=255)
                                gCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCa                                                                                                                    <  1:954479/150‑1 (MQ=255)
                                 cGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATAAAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCaa                                                                                                                   >  2:1035033/1‑150 (MQ=255)
                                  gTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAg                                                                                                                  <  1:516082/150‑1 (MQ=255)
                                     aaCATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAgtgt                                                                                                               >  2:411447/1‑150 (MQ=255)
                                       cATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTg                                                                                                             >  1:489746/1‑150 (MQ=255)
                                       cATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTg                                                                                                             >  1:729951/1‑150 (MQ=255)
                                         tACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGc                                                                                                           >  2:1663088/1‑150 (MQ=255)
                                         tACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGc                                                                                                           <  2:5469/150‑1 (MQ=255)
                                          aCAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCa                                                                                                          >  2:1528353/1‑150 (MQ=255)
                                           cAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCGGTGCGCAGCAcgatcaa                                                                                                                                       >  2:712436/1‑114 (MQ=255)
                                           cAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCGGTGCGCAGCAcgatcaa                                                                                                                                       <  1:712436/120‑7 (MQ=255)
                                             gTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCagag                                                                                                       <  1:1550669/150‑1 (MQ=255)
                                              tttCCACAAACTCCCTAACGGCGTGGTCAGCTAGTCAGTTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGt                                                                                                      <  2:1391890/150‑1 (MQ=255)
                                                 ccACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAAc                                                                                                   <  2:31833/150‑1 (MQ=255)
                                                  cacaAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAAcc                                                                                                  >  2:740376/1‑150 (MQ=255)
                                                  cacaAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAAcc                                                                                                  >  1:421968/1‑150 (MQ=255)
                                                    caAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAg                                                                                                                  >  2:1473685/1‑132 (MQ=255)
                                                    caAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAg                                                                                                                  <  1:1473685/132‑1 (MQ=255)
                                                      aaCTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACg                                                                                              >  1:873501/1‑150 (MQ=255)
                                                       aCTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGc                                                                                             >  1:161770/1‑150 (MQ=255)
                                                        cTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGcc                                                                                            >  1:12642/1‑150 (MQ=255)
                                                        cTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGcc                                                                                            <  2:785709/150‑1 (MQ=255)
                                                           ccTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAAt                                                                                         <  1:547470/150‑1 (MQ=255)
                                                           ccTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGATGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAAt                                                                                         >  2:1486674/1‑150 (MQ=255)
                                                             tAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAg                                                                                       <  1:672413/150‑1 (MQ=255)
                                                             tAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAg                                                                                       <  2:161770/150‑1 (MQ=255)
                                                                 ggCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAg                                                                                   >  2:956346/1‑150 (MQ=255)
                                                                    gtggtgAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACg                                                                                >  2:488535/1‑150 (MQ=255)
                                                                       gtgAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGAtt                                                                             >  2:39342/1‑150 (MQ=255)
                                                                        tgAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTc                                                                            >  1:522849/1‑150 (MQ=255)
                                                                           gCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTgc                                                                         >  1:188067/1‑150 (MQ=255)
                                                                            cGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTgcg                                                                        >  2:267001/1‑150 (MQ=255)
                                                                              aGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGcc                                                                      <  1:1316743/150‑1 (MQ=255)
                                                                                  aGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGtt                                                                  >  2:1274503/1‑150 (MQ=255)
                                                                                      aaTGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACt                                                              <  1:1528353/150‑1 (MQ=255)
                                                                                      aaTGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGAAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACt                                                              <  2:1154504/150‑1 (MQ=255)
                                                                                       aTGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACtt                                                             <  1:1104290/150‑1 (MQ=255)
                                                                                       aTGCTGAGACTTTGGTCACGCCTGAAGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTCTTGGCAGAGAAACCAACGCCAATAGCCAGACCATTCTCCGCCAGTTCAcat                                                             >  1:1529506/1‑148 (MQ=255)
                                                                                         gCTGAGACTTTGGTCACGACTGATGAAGACAGCTGATACAGTTGATCGTGCTGCGCACCGAGCAGAACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCAAATAGCCAGAAGATTCTGCGCCAGTTAACTTAt                                                           >  2:573322/1‑150 (MQ=255)
                                                                                          cTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTAtt                                                          >  2:14573/1‑150 (MQ=255)
                                                                                           tGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTAttt                                                         <  2:873501/150‑1 (MQ=255)
                                                                                            gagaCTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTc                                                        >  2:1314080/1‑150 (MQ=255)
                                                                                             agaCTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCa                                                       <  2:12642/150‑1 (MQ=255)
                                                                                             agaCTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCa                                                       <  1:1486674/150‑1 (MQ=255)
                                                                                               aCTTTGGTCACGCCTGATGACGACAGCTGATAAAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCAcc                                                     <  1:1035033/150‑1 (MQ=255)
                                                                                                       cACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATccc                                             >  1:1622833/1‑150 (MQ=255)
                                                                                                             tgatgaCGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCa                                       >  1:1153975/1‑150 (MQ=255)
                                                                                                             tgatgaCGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCa                                       >  1:61031/1‑150 (MQ=255)
                                                                                                             tgatgaCGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCa                                       >  1:647225/1‑150 (MQ=255)
                                                                                                                tgaCGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTc                                    <  1:411447/150‑1 (MQ=255)
                                                                                                                tgaCGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTc                                    >  2:1211190/1‑150 (MQ=255)
                                                                                                                      cAGCTGATACAGTTGATCGTGCTGCGCACCGCGCCGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGCCGATTCTGCGCCAGTTCACTTATTTCACCAGCATCCCGGTCCATTCACTccc                              >  2:921378/1‑150 (MQ=255)
                                                                                                                       aGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCt                             >  1:1166091/1‑150 (MQ=255)
                                                                                                                          tGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGcc                          >  2:731158/1‑150 (MQ=255)
                                                                                                                          tGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGcc                          >  2:442867/1‑150 (MQ=255)
                                                                                                                           gATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCg                         <  1:507790/150‑1 (MQ=255)
                                                                                                                            aTACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGAGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGc                        <  2:52750/150‑1 (MQ=255)
                                                                                                                             tACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCa                       <  2:449815/150‑1 (MQ=255)
                                                                                                                              aCAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCaa                      >  1:961845/1‑150 (MQ=255)
                                                                                                                                aGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAaca                    >  1:367861/1‑150 (MQ=255)
                                                                                                                                 gTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAacac                   >  1:781128/1‑150 (MQ=255)
                                                                                                                                 gTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAacac                   <  1:1213916/150‑1 (MQ=255)
                                                                                                                                    gATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACttt                >  1:1531565/1‑150 (MQ=255)
                                                                                                                                     ctggTGCTGCGCTGTGCGCAGCACCAGTTCTAACGGCCTCTCCTCCGCAAGAGTTGGCAGAGTGACCAACGCCAATAGCCACACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACATtg               >  2:927106/4‑150 (MQ=255)
                                                                                                                                     ctggTGCTGCGCGGTGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTtg               >  1:1095115/4‑150 (MQ=255)
                                                                                                                                     aTCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACTATTCTGCACCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTtg               >  2:117735/1‑150 (MQ=255)
                                                                                                                                      tCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTtgc              <  1:1314080/150‑1 (MQ=255)
                                                                                                                                       ggTGCTGCGCGGTGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTtgct             >  1:1249740/2‑150 (MQ=255)
                                                                                                                                        gTGCTGCGCGGTGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTtgctg            >  1:333642/1‑150 (MQ=255)
                                                                                                                                          gctgcGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTgcg          >  1:1501550/1‑150 (MQ=255)
                                                                                                                                          gctgcGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTgcg          <  2:188067/150‑1 (MQ=255)
                                                                                                                                            tgcGCACCGCGCAGCACCAGTTGTAACGGAGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGTTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCa        <  1:738720/150‑1 (MQ=255)
                                                                                                                                                cACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTg    <  1:213794/150‑1 (MQ=255)
                                                                                                                                                 aCCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGc   <  1:1290406/150‑1 (MQ=255)
                                                                                                                                                 aCCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGc   <  1:956346/150‑1 (MQ=255)
                                                                                                                                                  ccGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTc                                         <  1:774814/111‑1 (MQ=255)
                                                                                                                                                  ccGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTc                                         >  2:774814/1‑111 (MQ=255)
                                                                                                                                                  ccGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCa  <  1:1026287/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCA  >  NC_000913/2329563‑2329858

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: