Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 857,578 C→G 5.2% intergenic (+23/+218) ybiT → / ← ybiU ABC‑F family putative regulatory ATPase/DUF1479 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913857,5780CG 5.2% 608.9 / 17.0 172intergenic (+23/+218)ybiT/ybiUABC‑F family putative regulatory ATPase/DUF1479 family protein
Reads supporting (aligned to +/- strand):  ref base C (89/74);  new base G (4/5);  total (93/79)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

Only 100 of 174 total aligned reads displayed.

GTTTGGTAAGTTAATGATGCAGAAGCCGAACATTCTGATCATGGACGAACCGACCAACCACCTGGATATGGAATCCATTGAGTCGCTGAACATGGCACTGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAG  >  NC_000913/857306‑857837
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                   
                     gAAGCCGAACATTCTGATCAGGGAAGAACCGACCAAACACCTGGATGTGGAATCCATTGAGTCGCTGAACATGGAACTGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtgat                                                                                                                                                                                                                                               <  2:300783/274‑1 (MQ=255)
                                            aCGAACCGACCAACCACCTGGATATGGAATCCATTGAGTCGCTGAACATGGCACTGGAACTGCATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGCTTACCTGAGTAGGAACGGGATCGAGTCAGTGATCATCCGTCACCCCGATACTTTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCATGATgc                                                                                                                                                                                                                       >  1:49530/1‑275 (MQ=255)
                                                                                                  tGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGAGCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTTACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGAAATAATTTAACTCTGCAGGCCtt                                                                                                                                                                 >  4:312034/1‑275 (MQ=255)
                                                                                                                  cACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATgc                                                                                                                                                                                                                       >  2:165036/1‑205 (MQ=255)
                                                                                                                    cGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTTAGGTTTGCAGGCCTAAGCAGTACCTGCAATATACTGAATCTGGAGGACTTTGTAGGCCGGATGAGGCg                                                                                                                                               >  2:215238/1‑275 (MQ=255)
                                                                                                                    cGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAcc                                                                                                                                                                                           >  4:299222/1‑231 (MQ=255)
                                                                                                                    cGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACCCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCg                                                                                                                                               >  4:286872/1‑275 (MQ=255)
                                                                                                                      cTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTa                                                                                                                                                                                                                   <  4:165036/205‑1 (MQ=255)
                                                                                                                        gATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCTTACGCAGTACCTGCa                                                                                                                                                                                       <  2:299222/231‑1 (MQ=255)
                                                                                                                          tCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTTTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTCCGCTTATCAGGCCTACGCAGTACCTGCAATCTATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcg                                                                                                                                         >  4:185409/1‑275 (MQ=255)
                                                                                                                            tttGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAc                                                                                                                                                                                                                                                                     >  8:189641/1‑149 (MQ=255)
                                                                                                                               gTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGt                                                                                                                                                                                                                                                                  <  6:1163071/149‑1 (MQ=255)
                                                                                                                               gTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGt                                                                                                                                                                                                                                                                  <  5:1223702/149‑1 (MQ=255)
                                                                                                                                tttCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTc                                                                                                                                                                                                                                                                 >  6:1279639/1‑149 (MQ=255)
                                                                                                                                   cacaCGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCggg                                                                                                                                                                                                                                                              <  8:460230/149‑1 (MQ=255)
                                                                                                                                   cacaCGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCggg                                                                                                                                                                                                                                                              <  6:1364189/149‑1 (MQ=255)
                                                                                                                                    acacGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCgggg                                                                                                                                                                                                                                                             <  6:973705/149‑1 (MQ=255)
                                                                                                                                     cacGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGt                                                                                                                                                                                                                                                            >  8:1253840/1‑149 (MQ=255)
                                                                                                                                     cacGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGt                                                                                                                                                                                                                                                            <  7:501167/149‑1 (MQ=255)
                                                                                                                                     cacGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGt                                                                                                                                                                                                                                                            <  5:472208/149‑1 (MQ=255)
                                                                                                                                      acGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTg                                                                                                                                                                                                                                                           >  7:43995/1‑149 (MQ=255)
                                                                                                                                        gACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAc                                                                                                                                                                                                                                                         <  7:477796/149‑1 (MQ=255)
                                                                                                                                        gACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTGAc                                                                                                                                                                                                                                                         <  5:1147368/149‑1 (MQ=255)
                                                                                                                                          ccGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAcgg                                                                                                                                                                                                                                                       >  8:831957/1‑149 (MQ=255)
                                                                                                                                          ccGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAcgg                                                                                                                                                                                                                                                       >  7:898589/1‑149 (MQ=255)
                                                                                                                                            gTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTAGCCGTCACCCCGACAGTATCGGGGTGAcggat                                                                                                                                                                                                                                                     >  5:238072/1‑147 (MQ=255)
                                                                                                                                             tGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAcggcgg                                                                                                                                                                                                                                                    >  8:623127/1‑149 (MQ=255)
                                                                                                                                                gTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTGAcggataatc                                                                                                                                                                                                                                                 <  6:848083/149‑7 (MQ=39)
                                                                                                                                                 ttCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtga                                                                                                                                                                                                                                                <  5:560281/149‑1 (MQ=255)
                                                                                                                                                    gTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtgatga                                                                                                                                                                                                                                             <  7:1158715/149‑1 (MQ=255)
                                                                                                                                                    gTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATCTTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCAt                                                                                                               >  2:322997/1‑275 (MQ=255)
                                                                                                                                                        gCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGa                                                                                                                                                                                                                                                          <  6:941910/132‑1 (MQ=255)
                                                                                                                                                        gCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGa                                                                                                                                                                                                                                                          >  5:941910/1‑132 (MQ=255)
                                                                                                                                                        gCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGa                                                                                                                                                                              <  4:226947/208‑1 (MQ=255)
                                                                                                                                                         cTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAAt                                                                                                                                                                                                                                        <  8:205141/149‑1 (MQ=255)
                                                                                                                                                          tCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATg                                                                                                                                                                                                                                       >  5:742203/1‑149 (MQ=255)
                                                                                                                                                             cTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtgat                                                                                                                                                                                                                                               >  8:944468/1‑138 (MQ=255)
                                                                                                                                                             cTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCaa                                                                                                                                                                                                                                    >  7:407663/1‑149 (MQ=255)
                                                                                                                                                             cTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCaa                                                                                                                                                                                                                                    >  6:502211/1‑149 (MQ=255)
                                                                                                                                                               ggCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGGGCGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAAt                                                                                                    <  2:21388/274‑1 (MQ=255)
                                                                                                                                                                  gACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCggg                                                                                                                                                                                                                                                              >  5:610534/1‑118 (MQ=255)
                                                                                                                                                                   aCCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGc                                                                                                                                                                                                                              <  7:613897/149‑1 (MQ=255)
                                                                                                                                                                     ccGCATTCTGGAAATCACCCCGGAACGCCTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCt                                                                                                                                                                                                                            >  7:1253033/1‑149 (MQ=255)
                                                                                                                                                                        cATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGAt                                                                                                                                                                                                                         <  7:1253840/149‑1 (MQ=255)
                                                                                                                                                                           tCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATgcg                                                                                                                                                                                                                      <  7:1442627/149‑1 (MQ=255)
                                                                                                                                                                              ggAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTa                                                                                                                                                                                                                   >  5:628638/1‑149 (MQ=255)
                                                                                                                                                                              ggAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTa                                                                                                                                                                                                                   >  8:506957/1‑149 (MQ=255)
                                                                                                                                                                                   tCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATCCAAAATAGCCTGATGCGCTACGCtt                                                                                                                                                                                                              >  7:954285/1‑149 (MQ=255)
                                                                                                                                                                                    cACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTa                                                                                                                                                                                                             >  7:1384809/1‑149 (MQ=255)
                                                                                                                                                                                     aCCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTAt                                                                                                                                                                                                            >  6:410292/1‑149 (MQ=255)
                                                                                                                                                                                     aCCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTAt                                                                                                                                                                                                            >  5:274396/1‑149 (MQ=255)
                                                                                                                                                                                        ccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCaaaa                                                                                                                                                                                                                                  >  6:725004/1‑124 (MQ=255)
                                                                                                                                                                                        ccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAg                                                                                                                                                                                                         <  8:144203/149‑1 (MQ=255)
                                                                                                                                                                                          ggAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGc                                                                                                                                                                                                       <  5:1478283/149‑1 (MQ=255)
                                                                                                                                                                                            aaCGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCt                                                                                                                                                                                                     >  5:1344539/1‑149 (MQ=255)
                                                                                                                                                                                                 gTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGCGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCa                                                                                                                                                                                                <  8:1253033/149‑1 (MQ=255)
                                                                                                                                                                                                   gATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGt                                                                                                                                                                                              >  8:241336/1‑149 (MQ=255)
                                                                                                                                                                                                    aTCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTa                                                                                                                                                                                             <  5:502211/149‑1 (MQ=255)
                                                                                                                                                                                                     tCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAc                                                                                                                                                                                            >  8:485746/1‑149 (MQ=255)
                                                                                                                                                                                                       gACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCt                                                                                                                                                                                          >  5:16755/1‑149 (MQ=255)
                                                                                                                                                                                                        aCTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTg                                                                                                                                                                                         <  6:628638/149‑1 (MQ=255)
                                                                                                                                                                                                               cGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatata                                                                                                                                                                                  <  5:851027/149‑1 (MQ=255)
                                                                                                                                                                                                                 gTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATAtt                                                                                                                                                                                <  6:615728/149‑1 (MQ=255)
                                                                                                                                                                                                                    aTTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCCTCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGCATCGCTTTTCTTAATCCCACATCCATTCACTTCCATTaa                                               >  4:186783/1‑275 (MQ=255)
                                                                                                                                                                                                                       aCGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCt                                                                                                                                                                          <  7:883618/149‑1 (MQ=255)
                                                                                                                                                                                                                         gAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCGCCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAgg                                                                                                                                                                                                        <  7:643998/117‑1 (MQ=255)
                                                                                                                                                                                                                         gAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATAtt                                                                                                                                                                                <  8:131348/141‑1 (MQ=255)
                                                                                                                                                                                                                         gAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATAtt                                                                                                                                                                                >  7:131348/1‑141 (MQ=255)
                                                                                                                                                                                                                         gAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTGAcggat                                                                                                                                                                                                                                                     >  7:237959/1‑70 (MQ=25)
                                                                                                                                                                                                                          aaGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCa                                                                                                                                                                       >  6:247917/1‑149 (MQ=255)
                                                                                                                                                                                                                            gATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatata                                                                                                                                                                                  >  4:59176/1‑136 (MQ=255)
                                                                                                                                                                                                                             aTTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGa                                                                                                                                                                    >  5:1081791/1‑149 (MQ=255)
                                                                                                                                                                                                                                aCCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGa                                                                                                                                                                              <  2:59176/136‑1 (MQ=255)
                                                                                                                                                                                                                                  cTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTg                                                                                                                                                               <  6:7149/149‑1 (MQ=255)
                                                                                                                                                                                                                                   tGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGACGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGt                                                                                                                                                              <  7:1029535/149‑1 (MQ=255)
                                                                                                                                                                                                                                    gCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTa                                                                                                                                                             >  6:497854/1‑149 (MQ=255)
                                                                                                                                                                                                                                     cGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACTCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGACTAACGTCGAATTGATTTTCTTAATCCCCCATCAATTCACTTCCATTAATTCTCAAACATCCATg                               >  2:243108/1‑274 (MQ=255)
                                                                                                                                                                                                                                       tagtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTTNATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCAGCATCAGGCATTAACAAAGCGCATTTGTCANCAATCTGAGTCAAGTTGAATCGATTTTTTTAATCCAACAAAAATTNANTTCAACTAATTATAAAACAACCAttacc                            >  3:35069/1‑271 (MQ=37)
                                                                                                                                                                                                                                          taAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCgg                                                                                                                                                       >  8:710866/1‑149 (MQ=255)
                                                                                                                                                                                                                                           aaaGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGa                                                                                                                                                      >  5:1155360/1‑149 (MQ=255)
                                                                                                                                                                                                                                            aaGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGAt                                                                                                                                                     >  8:320851/1‑149 (MQ=255)
                                                                                                                                                                                                                                               ggATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAg                                                                                                                                                  >  5:1044309/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    gAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTa                                                                         >  2:53193/1‑217 (MQ=255)
                                                                                                                                                                                                                                                    gAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTAc               >  2:42683/1‑275 (MQ=255)
                                                                                                                                                                                                                                                       tAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAc                                                                                                                                          >  7:554816/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        aaGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATcc                                                                     <  4:53193/217‑1 (MQ=255)
                                                                                                                                                                                                                                                         aGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgc                                                                                                                                        >  8:108712/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         aGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGt            <  2:286872/273‑1 (MQ=255)
                                                                                                                                                                                                                                                           tGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgccg                                                                                                                                      >  7:1446338/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            gATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgccgc                                                                                                                                     >  6:828880/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 tCCGTCGCCCCTATTCTGTCGGGGGGCCGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGATAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAg  <  1:35069/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATAtt                                                                                                                                                                                <  6:1277605/101‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAg                                                                                                                                <  6:412220/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aCCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTa                                                                                                                          >  7:810746/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       aCCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTa                                                                                                                          >  8:646371/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        ccccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTaa                                                                                                                         <  5:828880/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ctcGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAAGAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAAc                                                                                                                        >  5:1079347/3‑149 (MQ=255)
                                                                                                                                                                                                                                                                         cccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAAc                                                                                                                        >  6:1070/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         cccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAAc                                                                                                                        >  5:524183/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                           cGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACaa                                                                                                                      >  6:354140/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                   
GTTTGGTAAGTTAATGATGCAGAAGCCGAACATTCTGATCATGGACGAACCGACCAACCACCTGGATATGGAATCCATTGAGTCGCTGAACATGGCACTGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAG  >  NC_000913/857306‑857837

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