Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 21.4% intergenic (+47/‑149) 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,5810AC21.4% 270.5 / 29.8 186intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (54/89);  new base C (7/36);  total (61/125)
Fisher's exact test for biased strand distribution p-value = 9.23e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 191 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  NC_000913/164432‑164730
                                                                                                                                                     |                                                                                                                                                     
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      >  2:2130887/1‑150 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                    >  2:5932491/1‑150 (MQ=255)
        tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGAACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                                               >  1:2544508/1‑150 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATACGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCAc                                                                                                                                            >  2:5057127/1‑150 (MQ=255)
               agttaTGAAAGGGCTTTATCCCAAACATGTCTGGCCGGACGACCGGTAAAATACTGCACCGACGCGAGGTACTAAAATTTATTCGTATGGTTGTTTTTCTCCCTCTCCCTGTGGTTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTtc                                                                                                                                        <  2:480800/146‑1 (MQ=255)
                 ataTGAAAGGGCTTTATCCCAAACATGTCGGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAAAGTATTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGTAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcac                                                                                                                                      <  2:2887253/148‑1 (MQ=38)
                  gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtggccc                                                                                                                                     >  1:3761199/1‑145 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggccc                                                                                                                                    >  1:1290734/1‑144 (MQ=255)
                    tGAAAGGGCTTTATCCAAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCTAGGGACTAAAAAGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGTTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccct                                                                                                                                   <  2:2386223/150‑1 (MQ=37)
                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTgcccccc                                                                                                                                   >  1:5274309/1‑149 (MQ=255)
                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtggggccc                                                                                                                                   >  2:1843656/1‑143 (MQ=255)
                       aaGGGCGTTATCCCAAACATGTCTGGCGGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTTTTCGTAAGTTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGCCATCAGCGCGCACGTtcaccctcac                                                                                                                                <  2:2964667/150‑1 (MQ=255)
                        cgggCGTTAACCCAACCCGGTCTGGCCGCCCGCCCCGCCAAATACGGCACCGAGGCGACGGACGAAAAAGTATTGGTACGGTTGTTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                               <  2:6237183/149‑1 (MQ=25)
                            cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCACCCaaa                                                                                                                           >  1:270340/1‑150 (MQ=38)
                            cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCaaa                                                                                                                           >  1:270325/1‑150 (MQ=255)
                              ttATCCCAAACTTGTCGGGCGGGACGACCCGGCAAATCCTGCACCGACGCGAGGTACTAAAATGTTTTGTTATGTTTGTTTTTCTCCCTCGCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                         <  2:580982/150‑1 (MQ=25)
                              ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGAACCACACCCaaaac                                                                                                                         >  1:3008147/1‑148 (MQ=255)
                              ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCCAAcc                                                                                                                         >  1:5711868/1‑150 (MQ=255)
                                  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCGCCCGCCCCCGAAccccca                                                                                                                     >  1:1996976/1‑147 (MQ=38)
                                  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCACCCCAAacccca                                                                                                                     >  2:3359900/1‑147 (MQ=37)
                                   ccAAACTTGTCGGGCGGGCCGACCGGGAAAATACTGCCCCGACGCGACGTACGAAAACGTATTGGTATGTTTGGTTTCCCCCCTCTCCGTGGGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                    <  2:5116368/150‑1 (MQ=255)
                                   ccAAACTTGGCTGGCGGTCCCACCGGGAAATTACTGAACCGACGCTTCGTACTAAAATTTTTTCGTATTGTTGGTTTTCTCCCTCTCCGGGTGGGTGTGGGCGGGGGTGTGGGAATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                    <  2:1293006/150‑1 (MQ=255)
                                   ccAAACATGTGTGGCCGGACGCCCGGGCAAATACTGCACCGACGCGAGGTACTAAAAAGTTTTCGTTTGGTTGTTTTTCTCCCTCGCCGTGTGGTTGTGGGCGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                    <  2:3337083/150‑1 (MQ=255)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACACCCGAAaccccac                                                                                                                    >  1:2895463/1‑146 (MQ=37)
                                    caaccctGTCTGGCCGGACGACCCGCCAAATCCGGCCCCGACGCGACGGACGAAACAGTTTTCGTTTGTTGGGTTTTCTCCCTCTCCCTGGGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                   <  2:6931764/144‑1 (MQ=255)
                                    cAAACTTGTGGGGCCGGACGACCCGGAAAATACTGCACCGACGCGCGGGACGAAAATTTATTCTAATGTTTGTTTTTCTCCGTCTCCGTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                   <  2:2290772/150‑1 (MQ=255)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCCCCGACGCGACGGACGAACAAGTATTCGTAAGGTTGTTTTTCTCCCTCACCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                   <  2:3784991/150‑1 (MQ=255)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAAcccacacc                                                                                                                   >  2:7163482/1‑145 (MQ=255)
                                     aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCACCCCAAaccacaacc                                                                                                                  >  2:2662183/1‑144 (MQ=37)
                                     aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAACCCACAccc                                                                                                                  >  2:5869259/1‑150 (MQ=255)
                                     aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCACACCCAAACCCACAccc                                                                                                                  >  1:1685777/1‑150 (MQ=38)
                                          tGTCTGGCCGGCCGACCCGGCAAATACTGCACCGCGGCGACGGACGAAAAAGTTGACGTACGTTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                             <  2:6242503/150‑1 (MQ=38)
                                          tGTCTGGCCGGCCGACCCGCAAAATACCCCACCGCCGCGCCGCCCGAAAAAGTTTTCGAAGGGTGGGTTTCCTCCCTCTCCCTGTGGGTTTGGGCCGGGGTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                             <  2:6132617/150‑1 (MQ=255)
                                          tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTAGTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                             <  2:2308009/150‑1 (MQ=255)
                                          ggTCGGGCGGCACGACCCGGAAATTACTGCACCGACGCGACGGACGAAAAAGTTTTCGTACGTTTGTTTTTCACCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                             <  2:5182982/149‑1 (MQ=37)
                                           gtggggccggccGCCCCGACAAAAACTCCCCCGCCGCGCGGTACTAAATTTTTTTCGTAGGGTTGTTTTTCCCCCTCGCCCTGTGGGTGGGGGGGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                            <  1:6804688/146‑1 (MQ=255)
                                           gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCttgagg                                                                                                            <  2:892051/150‑6 (MQ=255)
                                               ggccggccGCCCCGGAAAATCCGGCACCGCCGCGCCGGCCGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGGGGGTGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                        <  1:6392803/150‑1 (MQ=255)
                                               ggccggccGACCCGGAAAATACTGCACCGCCGCGACGGACGAAATAGTATTCTTAAGGTTGGTTTTCTCCCACTCCCTGGGGTTGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                        <  2:2174896/150‑1 (MQ=255)
                                               gcccGGACCACCCGGCAAAGACTCCCCCGACGCGCGGGCCGAAAATGTTTTCTTATGGTTGGTTTTCTCCCCCTCCCTGGGGGTGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                        <  2:6771636/148‑1 (MQ=255)
                                                gccggccgACCCGGCAAATACTGCACCGACGCGCCGGACGAAAATGTTTTCGTATGGTTGGTTTTCTCCCCCCCCCTGTGGGAGAGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                       <  2:5154517/150‑1 (MQ=255)
                                                gCCGGACGCCCCGGCAAACACGGCCCCGCCGCGACGGCCGAAAAAGTTTTCGTAAGGTTGGTTTTCTCCCTCCCCCTGTGGTCTAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                       <  2:5154713/150‑1 (MQ=37)
                                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACCCCCGAACACCCACCCCCACGgggggggg                                                                                                     >  1:4962556/1‑150 (MQ=25)
                                                   ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCAAACCCAAAACCAAAACCACAACGGAGAgggg                                                                                                    >  2:4105849/1‑150 (MQ=25)
                                                    gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGGAAGGGTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCCCCCAAAACCACACCCCCAACGGCGAGGGGa                                                                                                   >  2:1699595/1‑150 (MQ=25)
                                                     ccgacCCTCCAAATACCCCACCGACGCGCCGGACAAAATAGTTTTCGTATGTTTGGTTTTCCCCCCCTCCCTGGGGTAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggtta                                                                                                  <  1:767108/149‑4 (MQ=255)
                                                     acgacCCGGCAAATACTGCACCGACGCGACGGACTAAAAAGTTTTCGTAAGGTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                                  <  2:4016664/150‑1 (MQ=39)
                                                      ctacCCGGAAAATACTGCCCCGACGGTAGGGACTAAAATTTTTTGTTATGGTTGTTTTGCTCCCTCTGCGTGTGGGTTTGGGCGGGGTTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                                 <  1:4165252/148‑1 (MQ=255)
                                                       gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGGCCCCCACCCAAAACCACACCCACAACGGAGGCGGGGCCg                                                                                                >  1:6054994/1‑150 (MQ=25)
                                                        acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCCAACCCACACCCACAACGGCGCGGGGGccca                                                                                               >  2:5252051/1‑148 (MQ=37)
                                                          ccGGCAAAGACTGCACCGACGCGACGTACGAAAAAGTATGCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                             <  2:6577497/150‑1 (MQ=255)
                                                                      gCACGGACGCGACGGACGAAAATGTTTTCTTTCGTTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAttt                                                                                 <  1:5409374/150‑1 (MQ=37)
                                                                       cACCGCCGCGACGCACGAAAACGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGTGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                <  2:6848919/150‑1 (MQ=38)
                                                                        aCCTACGCTACGGACGAAAATGTATTCTTATGGTTGGTTTTCTCCCGCTCCCTTTGGTTGAGGGCCGGGGTTTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGTACCGATCGAGCACAAATTTTg                                                                               <  1:3943243/150‑1 (MQ=255)
                                                                         cccccgcgccgCACTAAAAAGTATGCTTAAGGTTTTTTTACGCCCTCCCGTGGGGGGTGAGGGGGGGGTTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  1:1302730/146‑1 (MQ=255)
                                                                         ccGAGGCGACTGTCGAAAATGTTTTGGTAGGGTTGGTTTTCGCCCCCGCGGGGGGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  1:4385780/150‑1 (MQ=255)
                                                                         ccGACGCGCCGGACCAAAAAGTATGCGTACGGTTGGTTTTCTCCCTCGCCCTGTGGGAGTGGGCCGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  1:6508618/150‑1 (MQ=255)
                                                                         ccGACGCGACGGACGAAAATGTATTCGTTAGGTTGGTTTTCTCCCTCGCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  1:7072439/150‑1 (MQ=255)
                                                                          cgacgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCCGCGCGCACGTGCACCCCCACCCAAACCCACTCCCACAAGGGGGAGGGGGCCGATCGAGCCCACAATTTgag                                                                             >  2:5947370/1‑150 (MQ=37)
                                                                           gacgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGGGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                            <  2:804064/150‑1 (MQ=255)
                                                                                gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGGCCCCCACCCAAACCCCCTACCCCACAGGCGAGGGGGCCGATCGAGCCCAAAATTTTAGGGaaat                                                                       >  2:1595694/1‑150 (MQ=21)
                                                                                 ccggactAAAATTTTTTCGTTAGGTTGGTTTTCCCCCTCCCCGGGGGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                      <  1:1569227/149‑1 (MQ=255)
                                                                                    gactAAAATTTATTGGTTTGGTTGTTTTTCTCCCTCCCCCTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                                                                   <  1:1157855/150‑1 (MQ=37)
                                                                                        aaaatGTATTCTTAAGTTTGGTTTACTCCAACTCCCTGTGGGAGTGGGCCGGGGTGCGGGAATCAGCGCGAACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                               <  2:713498/150‑1 (MQ=255)
                                                                                           aagaactCGTACGGTGGTATTACTCCCTCGCCGTGTGGGTGAGGGCCGTGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCATGGGTGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                            <  1:5159745/144‑1 (MQ=255)
                                                                                             tttttCGTATGGTTGGTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  2:4989802/147‑1 (MQ=255)
                                                                                             ttATTCGTATGTTTGGTTTCCACCCTCTCCCTGTGGGATGGGGCCGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  1:5877418/149‑1 (MQ=255)
                                                                                             gttttCGTTGGGTGGGTTTTCTCCCTCGCCCTGGGGGCGGGGGGGGGGGTGTGGGCATCCGCGCGCACGTGCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  1:925509/147‑1 (MQ=255)
                                                                                             gTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAAACCACACCCAACAGGGAGCGGGGGCCCAGCGAGCCCCAAAATTGAGAGATATAATATTCTGTGtt                                                          >  2:3472963/1‑150 (MQ=16)
                                                                                              tATTCTTATGGTTGTTTGTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                         <  1:6831894/150‑1 (MQ=39)
                                                                                               aTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCTCACCCGAACCCACACCCTCCAGGGGGAGGGGGCCGATCGAGCACAAAATTTGAGGTCTATCTTATTATGTCTTGt                                                        >  2:4542438/1‑150 (MQ=25)
                                                                                               aTTCGTAAGGTTGGTTTTCTCCCTCTCCCTCTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACAACCAAAACCAATCACGAAAACGGGAGGGGGCCGATAGAGCCCAAAATTTGTGAGATATATTATTCTGTGtttt                                                        >  2:5848903/1‑148 (MQ=12)
                                                                                                 tcttttggttgtttTTCCCCCCCTCCCTGTGGTCTGGGGCGGGGGTGCGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                      <  1:4872354/144‑1 (MQ=255)
                                                                                                 tcttatggttggttTGCTCCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                      <  1:6372632/144‑1 (MQ=255)
                                                                                                 tCGTATGTTTGGTTTCCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                      <  1:5962616/150‑1 (MQ=39)
                                                                                                 tCGTACGGTTGGTTTTCTCCCTCTCCCGGTGGGCGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                      <  1:6617035/150‑1 (MQ=255)
                                                                                                      tggttggttTTCTCCCTCTCCGTGGGGGAGGGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                 <  1:2381360/149‑1 (MQ=255)
                                                                                                         gtggttTTCGCCCTCCCCCTGGGGGATCGGGCGGGGGTGGGGCCCTCAGCGCGCACGTTCACCCTCCCCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                              <  2:1486967/149‑1 (MQ=255)
                                                                                                          tggttTTCCCCCCCTCCCTGGGGGAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                             <  2:5689461/150‑1 (MQ=255)
                                                                                                           gtttttCTCCCTCTCCCTGGGGGTGTGGGCGGGGTTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                            <  1:5634999/148‑1 (MQ=255)
                                                                                                           gttttgCTCCCCCGCCGTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                            <  2:2794412/148‑1 (MQ=255)
                                                                                                             ttttCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                          <  2:1094025/150‑1 (MQ=255)
                                                                                                                   ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGGGCACCCCCCCCCGAAACCTCGCCCGCAAGGGCGGGGGGGCCGCTCGCGCCCACAGTTTGAGAGATATCTTCGTCTGGCGTGTAACAGCAGCACAGAAAAACg                                    >  2:5346661/1‑150 (MQ=16)
                                                                                                                   ccctctccctGGGGGTGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                    <  1:7031549/150‑1 (MQ=255)
                                                                                                                   cccactccctGGGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGGTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                    <  1:1154335/150‑1 (MQ=255)
                                                                                                                    cctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCGCCCACACCCGGACCCTCTCCCTCAAGGGGGTGGGGACCGAGCGGGCGCAAAATTTGGGGGATATCTTGTTATGTGTTGTGAAAGACGTGCAGTAAATAgg                                   >  2:7144613/1‑150 (MQ=16)
                                                                                                                    cccctcccgGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                   <  1:5843845/147‑1 (MQ=255)
                                                                                                                           ctGTGGGTGGGGGCGGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                            <  1:5013843/150‑1 (MQ=255)
                                                                                                                               gggcgcgggCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                        <  1:764362/144‑1 (MQ=255)
                                                                                                                               gggCGAGGGCCGGGGTGAGGGCCTCAGCGCGCCCATGCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAGCGAGCACAACTTTTGAGAGATATCTTCTTCTGTCTAGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                        <  1:2148265/150‑1 (MQ=255)
                                                                                                                                ggtgtgggCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  1:6721122/145‑1 (MQ=255)
                                                                                                                                   gcgggCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATGTCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                    <  1:3269573/148‑1 (MQ=255)
                                                                                                                                        cgggggggCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg               <  2:2544508/148‑1 (MQ=255)
                                                                                                                                         cGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc              <  1:3997103/150‑1 (MQ=255)
                                                                                                                                          ggggTGCGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  1:373276/150‑1 (MQ=255)
                                                                                                                                           gggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCACAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGAGAACTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg            <  1:1803394/150‑1 (MQ=255)
                                                                                                                                                 gggCATCAGCGCGCACGTTCACCCACACACAAACCCTAACCCACCACGGCCAGGGGACCGATCCAGCACCAAATTTTAGAGATATAATCATATGTCTTTTAACAAAAGAAAAGAAAAAAGGGCGTTTGCGCCGGAATATTTCGGAGaaaa      >  2:3461702/1‑150 (MQ=255)
                                                                                                                                                  ggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa     <  1:5070034/150‑1 (MQ=255)
                                                                                                                                                     aTCAGCGCGCACGTTCACCCACAACCAAAACCAAACCCACCAGGGAGAGGGGGCCCATCGAGCACAAAATTTTAGATATAAATAAATATGTCTTGTAAAAAAAAAACACAAAAACGGGCTTTTTCGCCAGAATAATGCGGAGCAAAAACa  >  2:3061351/1‑150 (MQ=255)
                                                                                                                                                     aTCAGCGCGCACGTTCACCCAAACCCAAAACCAAAACCTCAAGGGAGAGGGGACCCATAGATCACAAAATTTGATATATATATTATTAAGTCTTTTAACAGAAGAACAGCAAAACGGGCTTTTTCGCCGGAAAAATTCGGGGAAAAAGCa  >  2:4682484/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  NC_000913/164432‑164730

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