Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,063,034 C→A 5.7% G248W (GGG→TGG)  cbpA ← DnaK co‑chaperone; curved DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,063,0340CA 5.7% 332.6 / 3.8 97G248W (GGG→TGG) cbpADnaK co‑chaperone; curved DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (39/52);  new base A (3/3);  total (42/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGACTGGGCGTCTGCCAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTCCAGATCCTGGCCGACAATATCAAACAGCGGATGTGGCG  >  NC_000913/1062889‑1063180
                                                                                                                                                 |                                                                                                                                                  
gactccTGGGCGTAGGCCAGTTGCTGCCACAGCGCGGCAGTTTTTTCATCCGGTTTCGGCGTCATCACTATTTTCAGTACCGAATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGc                                                                                                                                                 <  3:77638/144‑1 (MQ=255)
      tGGGCGTCTGCCAGTTGCTGCCACAGAGCGGCTGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCAGTTTTTTGCTAAACAGACCTTTGCCTTTAACGCGCAATCATTGCCCGGCCTGGCt                                                                                                                                           >  1:352956/1‑149 (MQ=255)
        ggCGTCTGCCAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGc                                                                                                                                         >  3:134051/1‑149 (MQ=255)
         gCGTCTGCCAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGcc                                                                                                                                        >  3:176178/1‑149 (MQ=255)
           gTCTGCCCGTTGCTGCCTGAGCGCGGCAGTGTTTTCTTCGGGGTTCGGCGGCATGACGATTTTGAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCTGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTg                                                                                                                                      <  3:254319/149‑1 (MQ=255)
               gCCAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCgg                                                                                                                                  >  8:229794/1‑149 (MQ=255)
                ccAGTTGCTGCCACAGCGCGGCTGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCGTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCggg                                                                                                                                 >  2:215822/1‑149 (MQ=255)
                ccAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGc                                                                                                                                    <  3:227486/146‑1 (MQ=255)
                ccAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCggg                                                                                                                                 <  2:361649/149‑1 (MQ=255)
                ccAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCgg                                                                                                                                  >  1:15228/1‑148 (MQ=255)
                ccAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGGCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCgg                                                                                                                                  <  5:237802/148‑1 (MQ=255)
                 cAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGa                                                                                                                                <  4:121570/149‑1 (MQ=255)
                       ctgcCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGt                                                                                                                            >  4:177635/1‑147 (MQ=255)
                        tgcCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAg                                                                                                                         >  3:339070/1‑149 (MQ=255)
                          ccACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAg                                                                                                                         <  7:104396/147‑1 (MQ=255)
                          ccACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGc                                                                                                                        >  2:408440/1‑148 (MQ=255)
                          ccACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGc                                                                                                                        >  8:44826/1‑148 (MQ=255)
                           cacaGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCa                                                                                                                       >  1:401936/1‑148 (MQ=255)
                             caGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCaaaa                                                                                                                    <  7:44826/149‑1 (MQ=255)
                                cgcgGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGc                                                                                                                 <  4:348248/149‑1 (MQ=255)
                                     cAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGctttct                                                                                                            <  5:243771/149‑1 (MQ=255)
                                      aGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGctttctt                                                                                                           <  2:356975/149‑1 (MQ=255)
                                      aGTGGTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCGGCCTGGCGGGATAGTCAGCAAAATGctttctt                                                                                                           <  1:325696/149‑1 (MQ=255)
                                       gtgtTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGctttcttt                                                                                                          <  3:177635/149‑1 (MQ=255)
                                          ttttCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCa                                                                                                        <  2:286488/148‑1 (MQ=255)
                                               aTCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTg                                                                                                  <  4:176178/149‑1 (MQ=255)
                                               aTCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTg                                                                                                  <  7:229794/149‑1 (MQ=255)
                                                  cGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGaa                                                                                               >  1:376967/1‑149 (MQ=255)
                                                       tCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGt                                                                                           <  5:33757/148‑1 (MQ=255)
                                                       tCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTAAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGt                                                                                           <  1:7639/148‑1 (MQ=255)
                                                             gcATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTa                                                                                    <  2:9003/149‑1 (MQ=255)
                                                                tCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATGGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTa                                                                                    <  7:20331/146‑1 (MQ=255)
                                                                tCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATGGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTa                                                                                    >  8:20331/1‑146 (MQ=255)
                                                                tCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCa                                                                                 >  1:191467/1‑149 (MQ=255)
                                                                  aCGATTTTCAGTTCCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCACAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGt                                                                                           <  2:85748/137‑1 (MQ=255)
                                                                  aCGATTTTCAGTGCCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCACAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGt                                                                                           >  1:85748/1‑137 (MQ=255)
                                                                  aCGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACttt                                                                                     <  6:214194/143‑1 (MQ=255)
                                                                  aCGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACttt                                                                                     >  5:214194/1‑143 (MQ=255)
                                                                      ttttCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGt                                                                                           >  6:231826/1‑133 (MQ=255)
                                                                      ttttCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGt                                                                                           <  5:231826/133‑1 (MQ=255)
                                                                         tCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAgcgc                                                                         <  1:408440/148‑1 (MQ=255)
                                                                           aGTACCGCATACAGATCGCCGGGCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGcgccg                                                                       <  5:340433/148‑1 (MQ=255)
                                                                            gTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCt                                                                     >  1:148557/1‑149 (MQ=255)
                                                                               ccGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTcc                                                                  >  3:281879/1‑149 (MQ=255)
                                                                                 gCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCa                                                                <  2:191467/149‑1 (MQ=255)
                                                                                  cATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAgt                                                                                                      <  7:235267/110‑1 (MQ=255)
                                                                                  cATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAgt                                                                                                      >  8:235267/1‑110 (MQ=255)
                                                                                    tACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACgg                                                             <  2:376967/149‑1 (MQ=255)
                                                                                     aCAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACgg                                                             <  4:134051/148‑1 (MQ=255)
                                                                                      cAGTTCGCCGGTCTTTATTATGCTCACCAGACCTTTGCCTTTAACGCGCAATCATTGCCCGGCCTGGCTGCCTGGCGGGATTGTCAGCAATATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACggg                                                            <  2:352956/148‑1 (MQ=255)
                                                                                      cAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACggg                                                            <  2:401936/148‑1 (MQ=255)
                                                                                       aGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCt                                                          <  5:379643/149‑1 (MQ=255)
                                                                                         aTCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACg                                                              >  6:201084/1‑143 (MQ=255)
                                                                                         aTCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACg                                                              <  5:201084/143‑1 (MQ=255)
                                                                                         aTCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTaa                                                        >  7:261453/1‑149 (MQ=255)
                                                                                              cGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAgt                                                                                                      >  1:383635/1‑98 (MQ=255)
                                                                                              cGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAgt                                                                                                      <  2:383635/98‑1 (MQ=255)
                                                                                              cGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCgg                                                    >  2:301387/1‑148 (MQ=255)
                                                                                              cGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCgg                                                    >  5:371752/1‑148 (MQ=255)
                                                                                              cGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGc                                                   <  1:301387/149‑1 (MQ=255)
                                                                                                     ttttttGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACaa                                            >  1:278246/1‑149 (MQ=255)
                                                                                                     accgttccaacactgaaagaaagcattttgctgactatcccGCCAGGCAGCCAGGCCGGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCt                                                                     >  1:385432‑M1/42‑124 (MQ=255)
                                                                                                     accgttccaacactgaaagaaagcattttgctgactatcccGCCAGGCAGCCAGGCCGGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCt                                                                     <  2:385432‑M1/83‑1 (MQ=255)
                                                                                                       ttttGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAAt                                           <  4:281879/148‑1 (MQ=255)
                                                                                                          tGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTc                                        >  8:198047/1‑148 (MQ=255)
                                                                                                           gCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCCGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTcc                                       <  1:215822/148‑1 (MQ=255)
                                                                                                            cTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTCCa                                      <  6:371752/148‑1 (MQ=255)
                                                                                                             aacactgaaagaaagcattttgctgactatcccGCCAGGCAGCCAGGCCGGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTc                                        <  7:35729‑M1/112‑1 (MQ=255)
                                                                                                               aCCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAg                                                                                   <  4:343455/100‑1 (MQ=255)
                                                                                                               aCCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTaa                                                        <  4:125716/127‑1 (MQ=255)
                                                                                                               aCCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTaa                                                        >  3:125716/1‑127 (MQ=255)
                                                                                                                 cAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTCCAGATCCt                                >  6:39869/1‑149 (MQ=255)
                                                                                                                  gatgtgaaagcattttgctgactatcccGCCAGGCAGCCAGGCCGGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGGTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTCCAGATCCTg                               >  3:253420‑M1/29‑149 (MQ=255)
                                                                                                                   gACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTCCAGATCCTgg                              <  7:198047/149‑1 (MQ=255)
                                                                                                                       tttGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGcac                                                 <  7:92801/126‑1 (MQ=255)
                                                                                                                       tttGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGcac                                                 >  8:92801/1‑126 (MQ=255)
                                                                                                                        aaagcattttgctgactatcccGCCAGGCAGCCAGGCCGGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTCCAGATCCTg                               <  4:253420‑M1/121‑1 (MQ=255)
                                                                                                                         tGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTc                                        >  3:191022/1‑133 (MQ=255)
                                                                                                                         tGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTc                                        <  4:191022/133‑1 (MQ=255)
                                                                                                                           ccTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTCCAGATCCTGGCCGACaa                       <  5:39869/148‑1 (MQ=255)
                                                                                                                            cattttgctgactatcccGCCAGGCAGCCAGGCCGGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTCCAGATCCCGGCCGACAata                     >  8:35729‑M1/19‑149 (MQ=255)
                                                                                                                             tttAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTccc                                                                 >  7:253530/1‑104 (MQ=255)
                                                                                                                             tttAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTccc                                                                 <  8:253530/104‑1 (MQ=255)
                                                                                                                               tAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTCCAGATCCTGGCCGACAATATc                   <  7:249032/148‑1 (MQ=255)
                                                                                                                                 aCGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTCCAGATCCTGGCCGACAATATCaaa                <  4:339070/149‑1 (MQ=255)
                                                                                                                                 aCGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTCCAGATCCTGGCCGACAATATCaaa                >  6:121358/1‑149 (MQ=255)
                                                                                                                                  cgcgCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTCCAGATCCTGGCCGACAATATCAAAc               <  2:148557/149‑1 (MQ=255)
                                                                                                                                     gcAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTCCAGATCCTGGCCGACAATATCAAACAGc            >  1:9233/1‑149 (MQ=255)
                                                                                                                                         tCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAAtt                                          >  8:213396/1‑115 (MQ=255)
                                                                                                                                         tCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAAtt                                          <  7:213396/115‑1 (MQ=255)
                                                                                                                                           gTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTCCa                                      <  5:344890/117‑1 (MQ=255)
                                                                                                                                           gTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTCCa                                      >  6:344890/1‑117 (MQ=255)
                                                                                                                                              GCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTCCAGATCCTGGCCGACAATATCAAACAGCGGATGTgg    >  2:278246‑M1/1‑148 (MQ=255)
                                                                                                                                              GCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGAACCCAGCGCCGCTTCCCACGGGCTAACCGTCACCACAATTTCCAGATCCTGGCCGACAATATCAAACAGCGGATTTGgc   <  2:165864‑M1/149‑1 (MQ=255)
                                                                                                                                               CCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTCCAGATCCTGGCCGACAATATCAAACAGCGGATGTGgcg  >  8:261453‑M1/1‑149 (MQ=255)
                                                                                                                                                CCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGc                                                   >  7:35624‑M1/1‑99 (MQ=255)
                                                                                                                                                CCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGc                                                   <  8:35624‑M1/99‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GACGACTGGGCGTCTGCCAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGATCGCCGGTCTGTTTTTTGCTCACCAGACCTTTGCCTTTAACGCGCAATCGTTGCCCGGCCTGGCTGCCTGGCGGGATAGTCAGCAAAATGCTTTCTTTCAGTGTTGGAACGGTGACTTTAGCACCCAGCGCCGCTTCCCACGGGCTAACCGGCACCACAATTTCCAGATCCTGGCCGACAATATCAAACAGCGGATGTGGCG  >  NC_000913/1062889‑1063180

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