Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 857,580 G→A 5.0% intergenic (+25/+216) 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,5800GA 5.0% 577.1 / 10.6 172intergenic (+25/+216)ybiT/ybiUABC‑F family putative regulatory ATPase/DUF1479 family protein
Reads supporting (aligned to +/- strand):  ref base G (87/74);  new base A (4/5);  total (93/79)
Fisher's exact test for biased strand distribution p-value = 7.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 173 total aligned reads displayed.

GTTTGGTAAGTTAATGATGCAGAAGCCGAACATTCTGATCATGGACGAACCGACCAACCACCTGGATATGGAATCCATTGAGTCGCTGAACATGGCACTGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAG  >  NC_000913/857306‑857837
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                 
                     gAAGCCGAACATTCTGATCAGGGAAGAACCGACCAAACACCTGGATGTGGAATCCATTGAGTCGCTGAACATGGAACTGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtgat                                                                                                                                                                                                                                               <  2:300783/274‑1 (MQ=255)
                                            aCGAACCGACCAACCACCTGGATATGGAATCCATTGAGTCGCTGAACATGGCACTGGAACTGCATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGCTTACCTGAGTAGGAACGGGATCGAGTCAGTGATCATCCGTCACCCCGATACTTTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCATGATgc                                                                                                                                                                                                                       >  1:49530/1‑275 (MQ=255)
                                                        acagccTGGATATGGGATCCATTGGGTCGCTGAACATGGCCCTGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTAt                                                                                                                                                                                                            <  4:78341/270‑1 (MQ=255)
                                                                                                  tGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGAGCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTTACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGAAATAATTTAACTCTGCAGGCCtt                                                                                                                                                                 >  4:312034/1‑275 (MQ=255)
                                                                                                                  cACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATgc                                                                                                                                                                                                                       >  2:165036/1‑205 (MQ=255)
                                                                                                                    cGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTTAGGTTTGCAGGCCTAAGCAGTACCTGCAATATACTGAATCTGGAGGACTTTGTAGGCCGGATGAGGCg                                                                                                                                               >  2:215238/1‑275 (MQ=255)
                                                                                                                    cGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACCCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCg                                                                                                                                               >  4:286872/1‑275 (MQ=255)
                                                                                                                      cTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTa                                                                                                                                                                                                                   <  4:165036/205‑1 (MQ=255)
                                                                                                                        gATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCTTACGCAGTACCTGCa                                                                                                                                                                                       <  2:299222/231‑1 (MQ=255)
                                                                                                                               gTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGt                                                                                                                                                                                                                                                                  <  5:1223702/149‑1 (MQ=255)
                                                                                                                               gTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGt                                                                                                                                                                                                                                                                  <  6:1163071/149‑1 (MQ=255)
                                                                                                                                tttCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTc                                                                                                                                                                                                                                                                 <  5:595284/149‑1 (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)
                                                                                                                                      acGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTg                                                                                                                                                                                                                                                           >  7:43995/1‑149 (MQ=255)
                                                                                                                                        gACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAc                                                                                                                                                                                                                                                         >  8:486730/1‑149 (MQ=255)
                                                                                                                                        gACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAc                                                                                                                                                                                                                                                         <  7:477796/149‑1 (MQ=255)
                                                                                                                                          ccGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAcgg                                                                                                                                                                                                                                                       >  7:898589/1‑149 (MQ=255)
                                                                                                                                          ccGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAcgg                                                                                                                                                                                                                                                       >  8:831957/1‑149 (MQ=255)
                                                                                                                                           cGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAcggc                                                                                                                                                                                                                                                      >  5:1137041/1‑149 (MQ=255)
                                                                                                                                            gTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTAGCCGTCACCCCGACAGTATCGGGGTGAcggat                                                                                                                                                                                                                                                     >  5:238072/1‑147 (MQ=255)
                                                                                                                                             tGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAcggcgg                                                                                                                                                                                                                                                    >  8:623127/1‑149 (MQ=255)
                                                                                                                                             tGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGAATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAcggcgg                                                                                                                                                                                                                                                    >  8:180200/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)
                                                                                                                                                        gCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACaaa                                                                                                                                                                                                                                         <  8:41540/149‑1 (MQ=255)
                                                                                                                                                        gCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGa                                                                                                                                                                              <  4:226947/208‑1 (MQ=255)
                                                                                                                                                          tCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATg                                                                                                                                                                                                                                       >  5:742203/1‑149 (MQ=255)
                                                                                                                                                             cTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtgat                                                                                                                                                                                                                                               >  8:944468/1‑138 (MQ=255)
                                                                                                                                                             cTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCaa                                                                                                                                                                                                                                    >  7:407663/1‑149 (MQ=255)
                                                                                                                                                               ggCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGGGCGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAAt                                                                                                    <  2:21388/274‑1 (MQ=255)
                                                                                                                                                                  gACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCggg                                                                                                                                                                                                                                                              <  6:610534/118‑1 (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)
                                                                                                                                                                                   tCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATCCAAAATAGCCTGATGCGCTACGCtt                                                                                                                                                                                                              >  7:954285/1‑149 (MQ=255)
                                                                                                                                                                                    cACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTa                                                                                                                                                                                                             >  7:1384809/1‑149 (MQ=255)
                                                                                                                                                                                     aCCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTAt                                                                                                                                                                                                            >  5:274396/1‑149 (MQ=255)
                                                                                                                                                                                     aCCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTAt                                                                                                                                                                                                            >  6:410292/1‑149 (MQ=255)
                                                                                                                                                                                       cccGGAACGCGTGATCGACTTTAGCGGTAATTATGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACCATTATCAGGCCTACGCAGTCCCTGCAATATATCGAATCTGCAGGACTTTGTAGTCCGGATAAGGCTTTAACGCCGCATCCGGCATTAACAAACCCAATTTTTCAACAATTTGAGTAACGTCGAACGGATTTTc                                                                            >  4:287072/1‑275 (MQ=255)
                                                                                                                                                                                        ccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAg                                                                                                                                                                                                         <  8:144203/149‑1 (MQ=255)
                                                                                                                                                                                            aaCGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCt                                                                                                                                                                                                     >  5:1344539/1‑149 (MQ=255)
                                                                                                                                                                                                 gTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGCGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCa                                                                                                                                                                                                <  8:1253033/149‑1 (MQ=255)
                                                                                                                                                                                                    aTCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTa                                                                                                                                                                                             <  5:502211/149‑1 (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                                                                                                                                                                                >  7:131348/1‑141 (MQ=255)
                                                                                                                                                                                                                         gAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATAtt                                                                                                                                                                                <  8:131348/141‑1 (MQ=255)
                                                                                                                                                                                                                         gAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGc                                                                                                                                                                        >  5:807880/1‑149 (MQ=255)
                                                                                                                                                                                                                         gAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGc                                                                                                                                                                        >  8:315305/1‑149 (MQ=255)
                                                                                                                                                                                                                         gAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTGAcggat                                                                                                                                                                                                                                                     <  8:237959/72‑3 (MQ=25)
                                                                                                                                                                                                                         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)
                                                                                                                                                                                                                                                 aTCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGc                                                                                                                                                <  5:411991/149‑1 (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)
                                                                                                                                                                                                                                                          gTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAg                                                                                                                                                                                               >  5:1353885/1‑93 (MQ=255)
                                                                                                                                                                                                                                                          gTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAg                                                                                                                                                                                               <  6:1353885/93‑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)
                                                                                                                                                                                                                                                                     tCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATgcgc                                                                                                                                                                                                                     >  4:171109/1‑60 (MQ=255)
                                                                                                                                                                                                                                                                       aCCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatatat                                                                                                                                                                                 <  5:990085/94‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aCCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTa                                                                                                                          >  7:810746/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       aCCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTa                                                                                                                          >  8:646371/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         ctcGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAAGAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAAc                                                                                                                        >  5:1079347/3‑149 (MQ=255)
                                                                                                                                                                                                                                                                         cccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAAc                                                                                                                        >  6:1070/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         cccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAAc                                                                                                                        >  5:524183/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          ccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACa                                                                                                                       >  6:732010/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                           cGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACaa                                                                                                                      >  6:354140/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                 tGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCa                                                                                                                >  5:831812/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                 
GTTTGGTAAGTTAATGATGCAGAAGCCGAACATTCTGATCATGGACGAACCGACCAACCACCTGGATATGGAATCCATTGAGTCGCTGAACATGGCACTGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAG  >  NC_000913/857306‑857837

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