Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 17.7% 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,5810AC17.7% 437.9 / 35.3 252intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (90/113);  new base C (9/40);  total (99/153)
Fisher's exact test for biased strand distribution p-value = 9.69e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 260 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  NC_000913/164432‑164729
                                                                                                                                                     |                                                                                                                                                    
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                     >  2:675544/1‑150 (MQ=255)
       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAgcgc                                                                                                                                               >  1:7764284/1‑150 (MQ=255)
         gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCgcgcgc                                                                                                                                             >  1:5380533/1‑150 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACGCGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGGAc                                                                                                                                           >  1:2024667/1‑150 (MQ=255)
             aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                                         >  1:2964139/1‑150 (MQ=255)
              aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGGGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtg                                                                                                                                        >  2:743108/1‑149 (MQ=255)
                 agaTTAAAGGGGGTTACCCCAAACTGGTGGGGCCGGACGACCCGCCAAATACTCCCCGGACGCGACGGACTAAAATGTATTCGTATGGTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTtcac                                                                                                                                     <  2:476279/150‑1 (MQ=255)
                  tattaatgggCTTTAACCAAAACCTGTCGGGCGGTACGACCCGGAAATTACTGCCCCGACGCGAGTTAATAAAATTTTTTCTTATGGTTGTTTTTCCCCCTCCCCCTGTGGTTTGGGGCGGGGTTGCGGGCATCAGCGCGCACGTtcacc                                                                                                                                    <  2:3142937/143‑1 (MQ=255)
                  gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggcc                                                                                                                                    >  1:9568471/1‑145 (MQ=255)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTtgaccccc                                                                                                                                 >  1:261030/1‑148 (MQ=255)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctc                                                                                                                                 >  2:7958573/1‑150 (MQ=255)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccgc                                                                                                                                 >  1:6242219/1‑148 (MQ=255)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGtcacccac                                                                                                                                 >  2:6443624/1‑148 (MQ=255)
                       aaGGGCGTTATCCCAACCTGGTCTGGCCGGACGACCCGTCAAATCCTGCACCGACGCGACGGACGAAAATTTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                               <  2:617059/150‑1 (MQ=37)
                         gggCGTTATCCCCAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAACAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGtgccccccaccc                                                                                                                             >  2:10199823/1‑150 (MQ=38)
                          ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCAGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgggcccccacccc                                                                                                                            >  1:3150405/1‑149 (MQ=38)
                              ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCCCAACCaaaac                                                                                                                        >  2:240274/1‑148 (MQ=255)
                              ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTGTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccgcccccacacc                                                                                                                        >  1:10745336/1‑145 (MQ=38)
                              ttATCCCAAACATGTCGGGCCGGACGACCCGGAAAATACGGCACCGACGCGACGGACGAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                        <  2:5380533/150‑1 (MQ=38)
                              ttATCCCAAACATGTCGGGCCGCACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                        <  2:1695908/150‑1 (MQ=255)
                              ttATCCCAAACAGGTCTGGCCGGACGACCCGGCAAATACGGCCCCGAGGCTAGTGACTAAAATTTTTTCTTATGTTTGTTTTTCTCCCTACCCCGGTGGGTGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                        <  2:293729/150‑1 (MQ=255)
                              ttATCCAAACCTTCTCGGGCCGCACGACCCGGCAAATACTGCACCGCCGCGCCGGACTAAAAAGTTTTCGTATGGTTGGTTTTCCCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                        <  2:5179103/150‑1 (MQ=25)
                                aTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCGCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCct                                                                                                                      <  2:9004316/150‑1 (MQ=255)
                                    caaccttGTCGGGCGGGACGCCCGGGCAAATCCGGCCCCGCCGCGTGGGCCGAAAATGTTTTCGTTGGGTGGGTTTTCTCCCTCGCCGTGGGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCtctcn                                                                                                                  <  1:10373179/144‑2 (MQ=255)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCCAAcccacacc                                                                                                                  >  2:8364988/1‑145 (MQ=38)
                                      aaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCCTCAGCGCGCACGTTCCCCCACCCCCAAACCCACAccca                                                                                                                >  1:1707669/1‑149 (MQ=38)
                                        cATGTCGGGCCGTCCGCCCCGCAAAAAACGCCCCCGCCGCGACGGACTAAAATTTATTCGTATGTTTGGTTTACCCCCACCCCCGGGGGTAGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                                                              <  2:2053471/150‑1 (MQ=255)
                                         aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:2684807/150‑1 (MQ=255)
                                          tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGAAAGTTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                            <  2:9084564/150‑1 (MQ=255)
                                          ggtggggCGGGACGACCCGGAAAATACTGCCCCGACGCTCCGGACTAAAATTTTTTGTTATGGTTGTTTTTCTCCCTCTCCGTGGGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                            <  2:9255213/145‑1 (MQ=255)
                                           gtggggCCGGACGCCCGGCCAAATAGGGCCCCGACGCGCGGTCCTAAAATTTTTTGGTAGGGTTGTTTTTCCCCGTCTCCCGGGGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                           <  2:3045612/146‑1 (MQ=255)
                                           gTCTGTCCGGCCGACCCGGCAAATACTGCACCGCCGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                           <  2:7296340/150‑1 (MQ=38)
                                            tCTGGCCGGACGACCCGGCAAATACTGCACCGACGAGACGGACGAAAAAGTATTCGTAAGGTTGGATTTCTCCCTCTCCCTGAGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACTCTCACCCTAACCCTCTCCCTCAAggg                                                                                                          >  2:515053/1‑150 (MQ=255)
                                              tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGCCCCACCCCCAAACCCCCACCCGCaccgggg                                                                                                        >  1:9826758/1‑144 (MQ=25)
                                              gggccggccGCCCCGGCAAATACGGCCCCGCCGCGCCGGCCGAACAAGTATTCGTACGGTGGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCCTCCGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTTGAGGgag                                                                                                        <  2:10037604/149‑1 (MQ=25)
                                               ggccggccTCCCCGCAAAATACTCCCCCGACGCGACTCACTAAAATTTATTGTTTTGGTGGTTTTACCCCCTCGCCGTGTGGGTGCGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGgaga                                                                                                       <  1:2137710/150‑1 (MQ=255)
                                               ggCCGGACCACCCGGCAAATACTGCACCTACGCGACGGACTAAAATGTTTTCGTATGGTTGTTTTTCTCCAACGCCCTGTGGGTGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                       <  2:573528/150‑1 (MQ=255)
                                                gCCGGACGCCCCGCAAAATCCGCCCCCTACGCTCCGGACTAAATATTATGCGTTTGGTTGTTTTTCTCCCGCTCCGTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                      <  1:2387064/150‑1 (MQ=255)
                                                 ccGTACGACCCGTAAAAAACTGCACCGACGCGCCGGACTAAAATGTATTCGTATGGTTGTTTTACCCCCTCCCCGTGTGGGTGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                     <  2:5006582/150‑1 (MQ=255)
                                                 ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGCGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                     <  1:9879184/150‑1 (MQ=38)
                                                 ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAACCCACACCCGCCAgggggcgg                                                                                                     >  1:5469430/1‑145 (MQ=37)
                                                    gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCTAACCCACTCCCACAAGGGGGAggggg                                                                                                  >  2:3662965/1‑149 (MQ=255)
                                                    gacgacCCGGCAAACACTGCACCCACGCGACGCCCGACAAAGTATTCGTAAGGTTGTTTTACCCCCCCACCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                  <  2:8486357/150‑1 (MQ=255)
                                                     acgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTGGGCCCGCACCCCAACCCACACCCACAAGGGGGAgggggg                                                                                                 >  1:10818156/1‑148 (MQ=37)
                                                     acgacCCGGAAAATACTGCCCCGCGGCGACGGACTAAAATTTATTCGTATGTTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                                 <  2:5997348/150‑1 (MQ=255)
                                                       gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCACCCCAAAACCACAACCACCAGGGGGAGGGGGCCg                                                                                               >  1:1865546/1‑150 (MQ=25)
                                                          ccGGCAAATACGGCACCGACGCGACGGACTAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                            <  2:8660648/150‑1 (MQ=38)
                                                             gCAAATCCGGCCCCGACGCGCCGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCCCCCCCGGGGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                         <  2:5811129/150‑1 (MQ=255)
                                                             gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGGGGCCCGCCCCCAAAACCACACCCCCCAGCGAGCGGGGCCCGAACGAg                                                                                         >  1:683384/1‑150 (MQ=21)
                                                             caaaaTCCGGCCCCGCCGCGCGGGACCAAAAATTTTTCGTAGGGTTGGTTTCCCCCCACCCCCGGGGGGTGTGGGCGGGGGTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                         <  1:6061820/148‑1 (MQ=255)
                                                               aaaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAACCCACACCCACCAGGGCGAGGGGGCCGATCGAGcc                                                                                       >  1:4639582/1‑149 (MQ=37)
                                                                       cacccccgcgacgGACGAAAATTTATCCTCATGGTTTTTTAACCCCCTAGCCCTGTGTTGGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                               <  1:4045932/144‑1 (MQ=255)
                                                                       cACCGACGCGTCGGCCGAAAATTTATTCGTATGGTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                               <  1:6568061/150‑1 (MQ=38)
                                                                       cACCGACGCGAGGGACGAAAAAGTTTTCGTATGTTGGGTTTTCTCCCTCCCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                               <  2:7875394/150‑1 (MQ=38)
                                                                       aagacagaagaagatatctctcaaaatttgtgctcgatcggtcccctctcccttgaGGGAGAGGGTTAGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                                    >  1:10262255‑M1/57‑97 (MQ=255)
                                                                         ccgccgCTAGGGACTAAAATGTATTCGTATTGTTGGTTTCCTCCCTCTCCGTGGGGGTGTGGGCGGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                             <  1:8513688/150‑1 (MQ=255)
                                                                         ccgccgCGCGGGACTAAAATTTTTTCGTATGGTGGGTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                             <  1:7252915/150‑1 (MQ=255)
                                                                         cccacgcgacgCACCAAAAAGTATTCGAAAGGTTGGTTTACACCCACCCCCTGTGGCAGAGGGCCGGGGTGAGGCAATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                             <  2:8796251/147‑1 (MQ=255)
                                                                           gcggcgacgTACTAAAAAGTTTTCTTATGGTTGTTTTCCGCCCTCCCCCGGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                           <  1:6852933/147‑1 (MQ=255)
                                                                            agaagaagatatctctcaaaatttgtgctcgatcggtcccctctcccttgaGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTa                                                                                                                           >  2:1535925‑M1/52‑101 (MQ=255)
                                                                            acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACAACCACAAGGGAGAGTGGACCCAACGAGCACAAAATTTtagag                                                                          >  2:2789576/1‑150 (MQ=25)
                                                                                   ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCACACCCCAACCCACAACCACAAGGGACAGGGGCCCGAACGACCCCAAAATTTTGGAGATAtctt                                                                   >  2:10836548/1‑150 (MQ=21)
                                                                                      cgAAAAAGTATTCTTACGGTTGGTTTCCCCCCTCTCCCTGTGGCAGAGGCCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcta                                                                <  1:5734634/150‑2 (MQ=255)
                                                                                        aaaaaGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                              <  2:10046251/150‑1 (MQ=255)
                                                                                        aaaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                              <  1:2686761/150‑1 (MQ=255)
                                                                                           attttttCGTTGGGTTGTTTTTCTCCCTCGCCGTGGGGGAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                           <  1:637321/145‑1 (MQ=255)
                                                                                           atgTATTCGTATGGTTTGTTTTCTCCCTCTCCCTGTGGGAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                           <  1:10659855/148‑1 (MQ=255)
                                                                                           aaGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                           <  1:7958573/150‑1 (MQ=255)
                                                                                             gTATTCGTAAGGTTGTTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                         <  2:3618661/150‑1 (MQ=255)
                                                                                               aTTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                       <  1:7361234/150‑1 (MQ=255)
                                                                                                ttCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGTAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa                                                      <  1:5805344/150‑1 (MQ=255)
                                                                                                 tCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCAACCGAACCCACTCCCACAAGGGAGAGGGGACCGATCGAGCACAAAATTTGAGTGATATCTTATTCTGTGTTGTaa                                                     >  1:1001509/1‑150 (MQ=37)
                                                                                                      cggttggttTACTCCCTCCCCCTGGGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGTTATCTTCTTCTGTCTTGTAACagaa                                                <  1:921224/149‑1 (MQ=255)
                                                                                                      aGGTTGTTTTTCGCCCTCTCCCTGTGGGTTTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                <  2:1281411/150‑1 (MQ=255)
                                                                                                       gtttggttTTCGCCCGCTCCGTGGGGGTGTGGGCGGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                               <  2:187998/148‑1 (MQ=255)
                                                                                                        gttgtttTTCCCCCCCCCCCTGTGGGAGCGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                              <  2:243285/150‑1 (MQ=255)
                                                                                                        gttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCCAACCCGCCCCCCCAAGGGAGAGGGGGCCGCGCGAGCCCCAAATTTTAGTGATATAATCATCGGTCGTGTAAAagaaga                                              >  2:2573804/1‑150 (MQ=18)
                                                                                                        gttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCTCACCCGAACCCACCCCCTCAAGGGAGAGGGGACCGATCGCGCCCAAAGTTTGTGAGATATCTTCTTCTGTCTTGTGACagaaga                                              >  1:8167551/1‑150 (MQ=37)
                                                                                                           ggttTACGCCCTCCCCGGGTGGGTGCGGGCCGGGGTGCGGCCCACAGCGCGCACGGTCACCCTCCCCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAACTTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                           <  1:9463696/150‑1 (MQ=255)
                                                                                                               ttCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACTACCGCAACGGGGAGGGGCCCCAACGAGCCCAAAATTTGTGAGATATATTATTATGTGTTGTAAAAGAAGAACAGaaa                                       >  2:6396835/1‑150 (MQ=21)
                                                                                                                 ctccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAACCCACACCCGCAAGGGGGAGGGGCCCCCTAGCGCACAAATTTTGGGAGATATAATATTCAGTCTTGGTACAGAAGAACAGtaaaa                                     >  2:4560660/1‑149 (MQ=19)
                                                                                                                     ctcgccgtggGGGAGTGGGGGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                                 >  1:4144460‑M1/11‑150 (MQ=255)
                                                                                                                            tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGGGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                          >  2:7051134/1‑150 (MQ=255)
                                                                                                                              ggggtgtgggCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                        <  1:10086408/143‑1 (MQ=255)
                                                                                                                              gGGGAGTGGGCGGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                        >  2:2304225‑M1/2‑150 (MQ=255)
                                                                                                                                gtcgtGGGCGGGGGGGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                      >  1:10563457‑M1/6‑150 (MQ=255)
                                                                                                                                 gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACAACCAAAACCACACCCACAACGGGGAGGGGGCCGCAAGAGCACAAAATTTGAGAGATTTATTAGTATGTCTTGTAACAAAAGAACACAAAAAAGGGCTTTTTCGccag                     >  2:7087920/1‑148 (MQ=255)
                                                                                                                                   GAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCCAACCCACTCCCACAAGGGGGAGGGGACCCATCGAGCACAAAATTTGAGTGATAAATTATTATGTCTTGTAACAGAAGAACAGAAAATAGGGCTTTTGCGCCGGaa                   <  2:8203458‑M1/150‑1 (MQ=255)
                                                                                                                                       gCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAACCCACACCCACCACGGGGAGGGGACCCAAACGGCACAAAATTTGGGTTGTATAATATTAAGTCGTGTAAAAGAAGCAAAGAAAAAACTGCGTTTGCGCCCGAAAAtt               >  2:9812629/1‑150 (MQ=255)
                                                                                                                                         cGGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc             <  2:5182419/150‑1 (MQ=255)
                                                                                                                                          ggggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  2:3382118/150‑1 (MQ=255)
                                                                                                                                          ggggTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  1:486328/150‑1 (MQ=255)
                                                                                                                                          ggggTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATACCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  1:1308967/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  2:3534802/150‑1 (MQ=255)
                                                                                                                                           gggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg           <  1:10426756/150‑1 (MQ=255)
                                                                                                                                            ggTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga          <  1:355232/150‑1 (MQ=255)
                                                                                                                                              ggcgggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCGGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgaga        <  2:6543634/147‑1 (MQ=255)
                                                                                                                                                  ggCCTCAGCGCGCACGTTCACCCTCACCCGACCCCTCTCCCTCAAGGGAGAGGGGACCGAGCGAGCACAAATTTTGAGAGATATCTTCTTCGGTCGTGTAACCGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa    <  2:10539784/150‑1 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGTTCACTCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGATATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa    <  1:515053/150‑1 (MQ=255)
                                                                                                                                                    cATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc  <  2:5600345/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  NC_000913/164432‑164729

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: