Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,621 A→G 5.4% intergenic (+87/‑109) 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,6210AG 5.4% 873.1 / 24.5 283intergenic (+87/‑109)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (94/162);  new base G (7/8);  total (113/170)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

Only 100 of 297 total aligned reads displayed.

CATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAA  >  NC_000913/164472‑164770
                                                                                                                                                     |                                                                                                                                                     
 ttGTCTGGCGGCACGACCCGCAAAATACTCCACCGACGCGACGCACTAAAATTTATTCGTATGTTTGTTTTTCTCCCTCCCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                                                      <  2:5030470/149‑1 (MQ=25)
 aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                                                      <  2:2684807/150‑1 (MQ=255)
  tgtggggccggccGACCCGCAAAATACTCCCCCGCCGCGCGGGACTAAAAATTTTTCTTTTGGTTGGTTTTCGCCCCCGCCGTGGGGTTGGGGGGGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                                                     <  1:6223955/145‑1 (MQ=255)
  tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGAAAGTTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                                                     <  2:9084564/150‑1 (MQ=255)
   gtggggCCGGACGCCCGGCCAAATAGGGCCCCGACGCGCGGTCCTAAAATTTTTTGGTAGGGTTGTTTTTCCCCGTCTCCCGGGGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                                                                    <  2:3045612/146‑1 (MQ=255)
   gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAACCTCTCCCTCAAgg                                                                                                                                                    >  2:1744447/1‑150 (MQ=255)
       cccgattttctgttcttctgttacaagacagaagaagatatctctcaaaatttgtgctcgatcggtcccctctcccttgagggagagggttagggtgagggtgaacgtgcgcgctgatggccgcacccccccccactcccccacggcgag                                                                                                                                                >  1:4144460‑M2/141‑146 (MQ=255)
         ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGCGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                                                              <  1:9879184/150‑1 (MQ=38)
         ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAACCCACACCCGCCAgggggcgg                                                                                                                                              >  1:5469430/1‑145 (MQ=37)
             acgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTGGGCCCGCACCCCAACCCACACCCACAAGGGGGAgggggg                                                                                                                                          >  1:10818156/1‑148 (MQ=37)
             acgacCCGGAAAATACTGCCCCGCGGCGACGGACTAAAATTTATTCGTATGTTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                                                                          <  2:5997348/150‑1 (MQ=255)
               gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCGCACCCAAAACCACAACCCCAAGCGGGAGGGGCCCg                                                                                                                                        >  1:8476966/1‑150 (MQ=25)
               gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCCCCCGAAACCCCAACCCCAAGGGGGTGGGGGCCg                                                                                                                                        >  1:1865563/1‑150 (MQ=25)
                  ccgtaaaaTCCTGCACCGCCGCGCGGGACCAAAACTTTTTGGTAGGGTTGGTTTGCCCCCCCCCCGTGGGGGTGCGGGCGGGGGGGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                                                                     <  2:3860258/145‑1 (MQ=255)
                     caaaaTCCGGCCCCGCCGCGCGGGACCAAAAATTTTTCGTAGGGTTGGTTTCCCCCCACCCCCGGGGGGTGTGGGCGGGGGTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                                                                  <  1:6061820/148‑1 (MQ=255)
                              gCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCACACCCACACCCGCTCCCGCCCGGGAGAGGGGGCCGAGCGCGCCCAAAatt                                                                                                                         >  2:10851448/1‑150 (MQ=25)
                               cacccccgcgacgGACGAAAATTTATCCTCATGGTTTTTTAACCCCCTAGCCCTGTGTTGGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                                                        <  1:4045932/144‑1 (MQ=255)
                               cACCGACGCGTCGGCCGAAAATTTATTCGTATGGTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                                                        <  1:6568061/150‑1 (MQ=38)
                                 ccgccgCGCGGGACTAAAATTTTTTCGTATGGTGGGTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                      <  1:7252915/150‑1 (MQ=255)
                                 ccgccgCGACGGACTAAAATGTTTTCGTATGTTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                      <  2:8476966/150‑1 (MQ=38)
                                 cccacgcgacgCACCAAAAAGTATTCGAAAGGTTGGTTTACACCCACCCCCTGTGGCAGAGGGCCGGGGTGAGGCAATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                      <  2:8796251/147‑1 (MQ=255)
                                 ccGACGCGCCGGACGAAAAAGTATTCGTACGTTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                      <  2:5403685/150‑1 (MQ=255)
                                 ccGACGCGACGTACGAAAAATTTTTCTTATGTTTGTTTTTCTCCCTCTCCGTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                      <  1:2683272/150‑1 (MQ=38)
                                 ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                      <  2:2789610/150‑1 (MQ=255)
                                   gcggcgacgTACTAAAAAGTTTTCTTATGGTTGTTTTCCGCCCTCCCCCGGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                                                                    <  1:6852933/147‑1 (MQ=255)
                                        gagggaccAAAATTTATTCTTACGTTTGGTTTTCCCCCACCCCGGGTGGGTGGGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                                                                                                               <  2:1972041/147‑1 (MQ=255)
                                             acgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCCCACCCAAAAACACACCCACAAGGGGGCGGGGCCCGATCGAGCCCAAATTTTTAGATATAtcatat                                                                                                          >  1:9941741/1‑146 (MQ=21)
                                                aaaacGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                                                                       <  2:10046109/150‑1 (MQ=255)
                                                   attttttCGTTGGGTTGTTTTTCTCCCTCGCCGTGGGGGAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                                                                    <  1:637321/145‑1 (MQ=255)
                                                   atgTATTCGTATGGTTGGTTTTCTCCCTCGCCCGGTGGGTGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                                                                    <  2:7958154/148‑1 (MQ=255)
                                                   aaGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                                                                    <  1:7958573/150‑1 (MQ=255)
                                                     gTATTCGTATGGTGGGTTTCCCCCCTCTCCCTGTGGGTGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                                  <  1:8834309/150‑1 (MQ=255)
                                                     gTATTCGTAAGGTTGTTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                                  <  2:3618661/150‑1 (MQ=255)
                                                      tATTCTTATGGTTGTTTTTCTCCCTCTCCCGGGGGGAGTGGGCGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAGTTTGAGAGAGATCTTCTTCTGTCTTg                                                                                                 <  2:4908414/150‑1 (MQ=255)
                                                         tcttatggttggttTTCTCCCTCTCCCTGTGGGTGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                                                              <  2:3501642/144‑1 (MQ=39)
                                                         tCGTATGGTTGGTTTCCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCACGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGATCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                                                              <  1:3112173/150‑1 (MQ=255)
                                                         tCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                                                              <  2:4639352/150‑1 (MQ=255)
                                                         tCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCAACCGAACCCACTCCCACAAGGGAGAGGGGACCGATCGAGCACAAAATTTGAGTGATATCTTATTCTGTGTTGTaa                                                                                              >  1:1001509/1‑150 (MQ=37)
                                                           gTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                                                            <  1:10740861/150‑1 (MQ=255)
                                                             acggttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                                                          <  1:8218235/148‑1 (MQ=255)
                                                              cggttggttTTCCCCCTCTCCCTGTTGGAGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                                                         <  1:2910042/149‑1 (MQ=255)
                                                              cggttggttTACTCCCTCCCCCTGGGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGTTATCTTCTTCTGTCTTGTAACagaa                                                                                         <  1:921224/149‑1 (MQ=255)
                                                                gttgtttTTCCCCCGCGCCCGGTGGGCGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                                                                       <  2:5155676/150‑1 (MQ=255)
                                                                gttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCCAACCCGCCCCCCCAAGGGAGAGGGGGCCGCGCGAGCCCCAAATTTTAGTGATATAATCATCGGTCGTGTAAAagaaga                                                                                       >  2:2573804/1‑150 (MQ=18)
                                                                   ggttTACGCCCTCCCCGGGTGGGTGCGGGCCGGGGTGCGGCCCACAGCGCGCACGGTCACCCTCCCCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAACTTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                                                                    <  1:9463696/150‑1 (MQ=255)
                                                                     ttttCTCCCTCTCCGTGGGGGTGTGGGGGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                                                                  <  1:7402179/150‑1 (MQ=255)
                                                                     ttttCTCCCTCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                                                                  <  2:5898068/150‑1 (MQ=255)
                                                                     ttgtcgccctctccctGTGGGAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                                                                  <  1:1213530/144‑1 (MQ=255)
                                                                           ccccccccctGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGAGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                                                            <  1:1031213/144‑1 (MQ=255)
                                                                                  cctGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                                                                     <  1:10034739/150‑1 (MQ=255)
                                                                                      gGGGAGTGGGCGGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                                                                 >  2:2304225‑M1/2‑150 (MQ=255)
                                                                                         gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACAACCAAAACCACACCCACAACGGGGAGGGGGCCGCAAGAGCACAAAATTTGAGAGATTTATTAGTATGTCTTGTAACAAAAGAACACAAAAAAGGGCTTTTTCGccag                                                              >  2:7087920/1‑148 (MQ=255)
                                                                                           GAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCCAACCCACTCCCACAAGGGGGAGGGGACCCATCGAGCACAAAATTTGAGTGATAAATTATTATGTCTTGTAACAGAAGAACAGAAAATAGGGCTTTTGCGCCGGaa                                                            <  2:8203458‑M1/150‑1 (MQ=255)
                                                                                              ggttagggtgagggtgaacgtgcgcgctgatgccctcaccccggccctctcccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat                                                         >  2:6887016‑M2/54‑150 (MQ=255)
                                                                                              ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACCCCCAAAACCACACCCGCAATGGGGAGGGGGCCGATCGAGCCCAAAATTTGAGATATATATTATTATGTCTTTTAAAAGAAGAAAAGAAAATACGGCTTTTTCGCCAGAaaat                                                         >  2:6401717/1‑150 (MQ=255)
                                                                                                 tagggtgagggtgaacgtgcgcgctgatgccctcaccccggccctctcccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                                                                       <  1:1535925‑M2/51‑1 (MQ=255)
                                                                                                 gggggTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                                      <  2:7858328/149‑1 (MQ=255)
                                                                                                 gggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                                      <  1:2324145/149‑1 (MQ=255)
                                                                                                  ggggTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATACCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                     <  1:1308967/150‑1 (MQ=255)
                                                                                                  ggggTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                     <  1:6841043/150‑1 (MQ=255)
                                                                                                  ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                     <  2:8591153/150‑1 (MQ=255)
                                                                                                  ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                     <  1:10570686/150‑1 (MQ=255)
                                                                                                   gggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg                                                    <  1:10426756/150‑1 (MQ=255)
                                                                                                    ggTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga                                                   <  2:1070668/150‑1 (MQ=255)
                                                                                                      ggcgggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCGGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgaga                                                 <  2:6543634/147‑1 (MQ=255)
                                                                                                          ggCCTCAGCGCGCACGTTCACCCTCACCCGACCCCTCTCCCTCAAGGGAGAGGGGACCGAGCGAGCACAAATTTTGAGAGATATCTTCTTCGGTCGTGTAACCGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa                                             <  2:10539784/150‑1 (MQ=255)
                                                                                                                   gcgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCg                                    <  2:9695559/150‑1 (MQ=255)
                                                                                                                    cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg                                   <  1:5165425/150‑1 (MQ=255)
                                                                                                                       aCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGaa                                <  2:898655/150‑1 (MQ=255)
                                                                                                                       aCGTTCACCCACACCCTAACCCTCTCCCACAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATATTCTTATGTCTTGTAACAGAAGAACAGAAAATAGGGCTTTTGCGCCTGAAAATTTCGGAGAAAAAGCATGGCCGGGaa                                >  2:6301448/1‑150 (MQ=255)
                                                                                                                        cGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAAt                               <  1:3057059/150‑1 (MQ=255)
                                                                                                                          ttCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGa                             <  2:10147577/150‑1 (MQ=255)
                                                                                                                            caccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGAcc                           <  1:9367442/150‑1 (MQ=255)
                                                                                                                                ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAGGCATGGCCGGGAATGACCGCGa                       >  1:5993787/1‑150 (MQ=255)
                                                                                                                                ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                       >  1:3331941/1‑150 (MQ=255)
                                                                                                                                ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                       >  1:213027/1‑150 (MQ=255)
                                                                                                                                     cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                  >  1:1822975/1‑150 (MQ=255)
                                                                                                                                     cctAACCCACTCCCTCAAGGGAGAGGGGACCGATCGAGCCCAAATTTTGAGAGATATCTTCTTCTGTCTTGTTACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAAGATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCac                  >  1:9064554/1‑149 (MQ=255)
                                                                                                                                      ctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAt                 >  1:7552014/1‑150 (MQ=255)
                                                                                                                                       tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt                >  2:5495536/1‑150 (MQ=255)
                                                                                                                                        aaCCCTCTCCCTCACGGGAGAGGGGACCGATCGAGCCCAAATTTTGAGAGATATCTTCTTCTGTCGTGGAACAGAAGAACAGAAAATCGGGCTGTTGCGCCTGAAGATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGAGccccttg               >  2:982554/1‑145 (MQ=255)
                                                                                                                                        aaCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg               >  1:4971037/1‑150 (MQ=255)
                                                                                                                                          ccctctcccTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt                <  2:7449421‑M2/138‑1 (MQ=255)
                                                                                                                                           cctctcccTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt                >  2:5692872‑M2/9‑146 (MQ=255)
                                                                                                                                           cctctcccTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt                <  1:5692872‑M2/138‑1 (MQ=255)
                                                                                                                                            ctctcccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg                                   <  2:8704012‑M2/119‑1 (MQ=255)
                                                                                                                                            ctctcccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg           <  2:9297903‑M2/143‑1 (MQ=255)
                                                                                                                                               tcccTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                                                              <  1:4403211‑M2/44‑1 (MQ=255)
                                                                                                                                               tcccTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                                                              >  2:4403211‑M2/5‑48 (MQ=255)
                                                                                                                                                 ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGAGCCCATTGGACGCAAAg      >  1:10270633‑M2/3‑150 (MQ=255)
                                                                                                                                                  cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg     >  1:5551394‑M2/2‑150 (MQ=255)
                                                                                                                                                  cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg     <  2:1323307‑M2/149‑1 (MQ=255)
                                                                                                                                                  cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg     <  1:5772187‑M2/149‑1 (MQ=255)
                                                                                                                                                   tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                                      <  1:4536757/100‑1 (MQ=255)
                                                                                                                                                   tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAggg    <  2:1498320/150‑1 (MQ=255)
                                                                                                                                                   tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAggg    >  1:4661653/1‑150 (MQ=255)
                                                                                                                                                     aaGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGaa  <  2:3107729/150‑1 (MQ=255)
                                                                                                                                                     aaGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGaa  >  1:4815372/1‑150 (MQ=255)
                                                                                                                                                     aaGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGaa  <  2:9826758/150‑1 (MQ=255)
                                                                                                                                                     aaGGGAGAGGGGACAGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGaa  <  2:2930916/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAA  >  NC_000913/164472‑164770

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: