Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,594 T→G 5.2% intergenic (+60/‑136) hrpB → / → mrcB putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,5940TG 5.2% 675.7 / 2.1 259intergenic (+60/‑136)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (61/178);  new base G (12/8);  total (73/186)
Fisher's exact test for biased strand distribution p-value = 2.94e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 265 total aligned reads displayed.

AAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGA  >  NC_000913/164445‑164743
                                                                                                                                                     |                                                                                                                                                       
aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                                                         >  2:9828039/1‑150 (MQ=255)
  agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgg                                                                                                                                                       >  1:10383707/1‑148 (MQ=255)
       tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggccca                                                                                                                                                  >  2:5771373/1‑143 (MQ=255)
         aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctca                                                                                                                                                <  2:10063433/150‑1 (MQ=255)
         aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGAGAGGGCATCCGCGCGCACGTtcacccacc                                                                                                                                                >  2:11005171/1‑147 (MQ=255)
             ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTgggcccccacccc                                                                                                                                            >  1:8406921/1‑149 (MQ=38)
                tttATCCCAAACTGGGGGGGCGGGACGACCGGGAAAATACTGCACCGACGCGAGGGACTAAAAAGTATTGTTATGGTTGGTTTTCTCCCTCCCCCTGTGGTTGAGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                                         <  2:9790605/149‑1 (MQ=255)
                gTTATCCCAACCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGATGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCCCGTTCACCCTCACCCTAAc                                                                                                                                         <  2:1452506/150‑1 (MQ=38)
                 ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCGCACCCAAAcc                                                                                                                                        >  1:1764662/1‑150 (MQ=255)
                 ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCACACCCCAAcc                                                                                                                                        >  2:7942323/1‑150 (MQ=255)
                 taaacccAAACATGACTCCCCCTACGACCCCCCAAAAACCCCACCCACCCCACGCACCAAAAACTATTCGAAAGGTTGGATTACACCCACACCCTGTGGCAGAGGCCCGGGGTGAGGGCAACACCGCCCACGTTCACCCTCACCCTAAcc                                                                                                                                        <  2:7240832/146‑1 (MQ=255)
                        aaaCATGTCTGGCAGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTAACCCACAACCAAAaccacacac                                                                                                                                 >  1:6623919/1‑144 (MQ=37)
                            aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTCTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGCACCCGCACCCtcaccctcacccgccc                                                                                                                             >  2:11122598/1‑146 (MQ=38)
                             tGTCTGGCCGGACGACCCGGCAAATACTGCCCCGACGCGACGGACGAAAAAGTATTCTTATGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCTCGTTCACCCTCACGCTAACCCTCTCCCTCAAg                                                                                                                            <  2:8216097/150‑1 (MQ=38)
                               tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTCGTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTGCACCCGCACCCGAACCCCCACCCACACggg                                                                                                                          >  2:4793137/1‑150 (MQ=37)
                               tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGTTCGCCCACACCCACACCCACACCCGCAAggg                                                                                                                          >  2:7901970/1‑150 (MQ=37)
                                  ggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAATTTTTTGGTATGGTTGTTTTTCTCCCTCGCGCTGTGGTCGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                                       <  1:9157925/150‑1 (MQ=255)
                                  ggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCACCCTCTCCCGGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                                       <  2:821401/150‑1 (MQ=255)
                                   gcgggccgacCCGGCAAATCCTGCACCGCCGCGACGGACGAAAAAGTTTTCGTATGGTTGGTTTTCTCCCTCGCCCGGTGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgagat                                                                                                                      <  2:5423386/144‑2 (MQ=255)
                                    ccGGACGCCCCGCAAAATCCTCCACCGCCGCGCCGCACTAAAATTTTTGCTTTTGGTTGGTTTCCCCCCTCTCCGTGTGGGAGTGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                                     <  1:6398284/150‑1 (MQ=255)
                                    ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCACCCAAAACCACACCCCCCAgggggggg                                                                                                                     >  1:1202921/1‑145 (MQ=25)
                                    ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCACCCACACCCTAACCCACACCCCCCAgggggcgg                                                                                                                     >  2:4507393/1‑145 (MQ=37)
                                     cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGGGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                                    <  2:5843079/150‑1 (MQ=255)
                                         cgacCCGGCAAATACTGCACCGACGCTCCGGACTACAAATTATTCGTATGGTTGGTTTTCTCCCTATCCGGGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                                                <  2:5692465/150‑1 (MQ=255)
                                          gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCACAAACCCGAAACCACAACCAAAAGGGAGATGGGAcag                                                                                                               >  1:5630358/1‑148 (MQ=25)
                                             ccGGCAAAACCACCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                                            <  2:9692227/150‑1 (MQ=38)
                                             ccGGCAAAAACTGCACCGACGCGAGGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                                            <  2:6451638/150‑1 (MQ=255)
                                              cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCCAAACCACAACCCCAACGGCGAGGGGGcccaaag                                                                                                           >  2:4322513/1‑145 (MQ=25)
                                              cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCAAACCCTCTCCCCCAACGGGGGGGGGGCCgcgcg                                                                                                           >  2:7364037/1‑146 (MQ=37)
                                                gaaaaTCCTGCCCCGACGCGCCGGACTAAAAATTATTCGTAAGGTTGGTTTACTCCCTCGCCCTGTGGTAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                                         <  1:6783149/148‑1 (MQ=255)
                                                gCAAATACTGCACCGACGCGACGGACGAAAAAGTGTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCACCCACAACCCCACCCCCACCGGGGAGGGGGCCgggaggg                                                                                                         >  1:2340931/1‑144 (MQ=21)
                                                  aaaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCCCACCCAAAACCACACCCGCACGGGGGGTGGGGCCCATCGAGcc                                                                                                       >  1:4633865/1‑149 (MQ=25)
                                                  aaaTACTGCACCGACGCGACGCACTAAATTTTATTCGTATGGTTGGTTTTCCCCCACCCCCTGTGGGTGTGGGCGGGGTTGAGCGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                                       <  1:10829495/150‑1 (MQ=255)
                                                       cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCGAACCCTCGACCCCCAGGGCGAGGGGGCCGATCGAGCAC‑‑ACAt                                                                                                >  2:7283855/1‑150 (MQ=37)
                                                          ccccgccgCTACGCACTAAAATTTTTTCGTCCGGTTTTTTTTCCCCCCCCCCGTGTGGGCGGGGGGGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAAtttt                                                                                             <  1:7439/148‑1 (MQ=255)
                                                          cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCGCACCCGAACCCACTCCCACAAGGGCGAGGGGGCCGAGCGAGcac‑‑acaattg                                                                                             >  2:7299431/1‑144 (MQ=37)
                                                            ccgccgCGACGGCCGAAAAATTATTGGTAGGGTTGGTTTTCTCCCTCCCCGGGTGGGTGCGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTga                                                                                           <  2:5231854/150‑1 (MQ=255)
                                                            ccGACGCGCCGGACTAAAAAGTTTTCGTCCGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTga                                                                                           <  1:9889806/150‑1 (MQ=255)
                                                            ccGACGCGACGTACTAAAATTTATTCTTATGGTGGTTTTCCGCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTga                                                                                           <  1:6938040/150‑1 (MQ=37)
                                                            ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCAAACCCAAAACCACACACACAAGGGGGAGGGGGCCGAAAGAGCCC‑‑AAAatttta                                                                                           >  1:9098699/1‑148 (MQ=21)
                                                                 gcgacgTACTAAAATGTTTTCTTTTGGTTGTTTGTCTCCCTCTCCGTGTGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGat                                                                                      <  1:4287811/150‑1 (MQ=255)
                                                                 gcgacgGACTAAAAAGTATTCGTAAGGTTGGATTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGat                                                                                      <  1:9973545/150‑1 (MQ=255)
                                                                 gcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGat                                                                                      <  2:10056570/150‑1 (MQ=255)
                                                                    acggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATAtc                                                                                   <  2:10905445/150‑1 (MQ=255)
                                                                     gggactAAAATTTATTCTTATGGTTGTTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATAtct                                                                                  <  2:8706149/149‑1 (MQ=255)
                                                                      ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCGCCCCCAAAACCACGACCGCCAGGGGGGGGGGGGCGGGCGGGCACAAAAATTTTGAGatatatt                                                                                   >  2:10525351/1‑149 (MQ=12)
                                                                       gacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATAtcttc                                                                                <  2:7068672/150‑1 (MQ=255)
                                                                       gacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCCAACCCAAACCCACCAGGGAGAGGGGTCCCATACAGCAC‑‑AAAATTTTAGAGATAtattc                                                                                >  2:9303307/1‑150 (MQ=25)
                                                                       gaccAAAATGTTTTCTTACGGTGGGTTCGCCCCCTCGCCGGGGGGGTGCGGGCGGGGGGTCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATAtcttc                                                                                <  1:3016227/150‑1 (MQ=255)
                                                                          gAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGAATCAGCGCGCACGTTGACCCACACCCAAAACCAAAACCCCAAGGGAGGGGGGCCCGAGCGACCCC‑‑AAAGTGTTTGAGATAacttattc                                                                             >  2:3098820/1‑150 (MQ=16)
                                                                                gTATTCTTACGGTTGGTGTTCTCCCGCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCtt                                                                       <  1:8127041/150‑1 (MQ=39)
                                                                                  tttCTTACGGTGGGTTTACTCCCTCTCCCGGTGGGTGTGGGCGGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                                     <  2:9405387/149‑1 (MQ=255)
                                                                                  atgcGTACGGTGGGTTTTCTCCCTCCCCCTGTGGGAGTGGGCGGGGGGGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                                     <  2:1071062/147‑1 (MQ=255)
                                                                                    tCGTATGGGTGTTTTACTCCCACTCCCGGTGGTAGAGGGCCGGGGTGCGGGACTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                                   <  2:373888/150‑1 (MQ=255)
                                                                                     cGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                                  <  1:10083126/150‑1 (MQ=255)
                                                                                      ttttggttgtttTTCTCCCTCGCCGGGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                                 <  2:8384/146‑1 (MQ=255)
                                                                                      gtttggtgggttTGCTCCCTCGCCGGGGGGGTGGGGGCGGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                                 <  1:9501670/146‑1 (MQ=255)
                                                                                      gTACGGTTGGTTTTCTCCCTCGCCCTGTGGGAGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                                 <  2:21645/150‑1 (MQ=255)
                                                                                      gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCCAACCCACTCCCCCAAGGGCGAGGGGACCCATAGAGCAC‑‑AAAAATTGATAGATATCTTCTTCTGTCTTGTAAcc                                                                 >  2:11013837/1‑149 (MQ=37)
                                                                                        ttgtttgtttTCCTCCCTCTCCGTGTGGTTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                               <  1:4030417/146‑1 (MQ=255)
                                                                                        atggttgtttTCCCCCCTCGCCGTGGGGTTGGGGGGGGGGGTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                               <  1:532092/148‑1 (MQ=255)
                                                                                          ggttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                                             <  2:6455074/150‑1 (MQ=255)
                                                                                           gttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCTCACCCTAACCCTCTCCCTCAAGGGGGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTTACagaaga                                                            >  1:711465/1‑150 (MQ=255)
                                                                                            tgggttTTCTCCCTCTCCCTGTGGGCGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                                           <  1:3423140/148‑1 (MQ=255)
                                                                                             tggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                                          <  2:2793564/150‑1 (MQ=255)
                                                                                               gtgttCGCCCTCCCCCGGGGGGAGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAg                                                        <  2:3451114/147‑1 (MQ=255)
                                                                                                     gccctctcggtGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                                  <  2:7645810/149‑1 (MQ=255)
                                                                                                      ccctctccctGTGGGAGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAACTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                                 <  2:6974794/150‑1 (MQ=255)
                                                                                                        ctctccctGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                                               <  1:1817178/150‑1 (MQ=255)
                                                                                                         tctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCTCCCGCACCCACACCCTCTCCCACAAGGGGGTGGGGCCCGGTCGAGCAC‑‑AAATTTTTAGAGATATAATATTAGGTCTTGTAACAGAAGGACAGTAAATCGGGc                                              >  2:2680202/1‑150 (MQ=21)
                                                                                                             cctGGGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                                          <  1:10720002/150‑1 (MQ=255)
                                                                                                                 tgggtgggGGCGGGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                                      <  2:6835720/150‑1 (MQ=255)
                                                                                                                   ggtgtgggCGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                    <  2:3217462/145‑1 (MQ=255)
                                                                                                                   ggAGAGGGCCGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCACTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                    <  1:1051648/150‑1 (MQ=255)
                                                                                                                     agagGGCCGGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                                  <  2:1743938/150‑1 (MQ=255)
                                                                                                                       cgggCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAat                                <  1:7580389/149‑1 (MQ=255)
                                                                                                                        gggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAata                               <  1:11131818/150‑1 (MQ=255)
                                                                                                                         ggCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTGGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat                              <  2:8884760/150‑1 (MQ=255)
                                                                                                                         ggCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCCCACCCAAAACCAAAACCCCAACGGAGAGGGGACCGAACGAGCCC‑‑CAAGTTTGTGAGATATATTGATATGTGTTGTTACAGAACAACACTCAAACGGGCGTTTGCGccggcagat                              >  2:2572252/1‑143 (MQ=255)
                                                                                                                         ggCCGGGGTGAGGGCATCAGCGCGCACGTGAACCCACAACCGAAACCACACCCACAAGGGGGAGGGGGCCCATCGCGCAC‑‑AAAAATTGAGTGATATAATATTAGGTGATTTAAAACAAGAAAATAAAAAAGGGCGATTGCGCCTGGatat                              >  2:3598108/1‑150 (MQ=255)
                                                                                                                         ggCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCCCACCCGAACCCGCTCCCGCCAGGGCGAGCGGCCCCCTCCCGCCC‑‑CAATGGTGAGAGATAAATTCTTCTGGCGGGGAACAGCAGAACCGCACAGCGGGCGGTTGCGCCGGacgag                              >  2:10122746/1‑146 (MQ=255)
                                                                                                                          gCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt                             <  1:10189294/150‑1 (MQ=255)
                                                                                                                            cGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                           <  1:7653905/150‑1 (MQ=255)
                                                                                                                            cGGGGTGAGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                           <  2:8997430/150‑1 (MQ=255)
                                                                                                                             ggggTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGGAGAACTGAAAATCGGGCTTTTGCGCCTGAATATTGCg                          <  2:5110124/150‑1 (MQ=255)
                                                                                                                             ggggTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                          <  1:11005171/150‑1 (MQ=255)
                                                                                                                             ggggTGCGGGCCTCAGCGCGCACGTTCACCCTCACGCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                          <  2:7557684/150‑1 (MQ=255)
                                                                                                                             ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                          <  1:5923829/150‑1 (MQ=255)
                                                                                                                             ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                          <  1:8331718/150‑1 (MQ=255)
                                                                                                                               ggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga                        <  2:8106545/150‑1 (MQ=255)
                                                                                                                               ggTGAGGGCCTCAGCGCGCACGTGCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCCGAATATTGCGga                        <  2:1202587/150‑1 (MQ=255)
                                                                                                                                 tGAGGGCATCAGCGCGCACGTGCACCCACCCCCCAACCCACGCCCCCCACGGAGCGGGGCCCCCGCGAGCAC‑‑AAAATTTGGGGGATATATTCATCGGTCTTGTAACACAAGAACAGCAAATACGGCGGTTGCGCCCGCATATTGCGggga                      >  2:2613292/1‑150 (MQ=255)
                                                                                                                                               cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg        <  2:2285871/150‑1 (MQ=255)
                                                                                                                                               cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAACGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGTAGAAAAAGCATGGCCgg        <  1:1485868/150‑1 (MQ=255)
                                                                                                                                                 cACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGa      <  1:588296/150‑1 (MQ=255)
                                                                                                                                                  aCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGaa     <  1:4112752/150‑1 (MQ=255)
                                                                                                                                                   cGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAAt    <  1:9051238/150‑1 (MQ=255)
                                                                                                                                                    gTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATGTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATg   <  2:8856373/150‑1 (MQ=255)
                                                                                                                                                     ttCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGa  <  1:6204488/150‑1 (MQ=255)
                                                                                                                                                     ttCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGa  <  1:10006682/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
AAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGA  >  NC_000913/164445‑164743

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