Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,621 A→G 7.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 7.4% 543.6 / 25.3 181intergenic (+87/‑109)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (48/115);  new base G (9/4);  total (61/120)
Fisher's exact test for biased strand distribution p-value = 5.34e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.35e-01

Only 100 of 185 total aligned reads displayed.

ATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAA  >  NC_000913/164473‑164770
                                                                                                                                                     |                                                                                                                                                     
tttagtGGCCGGCCCCCCGGGCAAATCCTGCCCCGACGCCCCGTACGAAAAAGTTTTCTTATGTTGGGTTTTCCCCCCCGCCCGGGGGGTGTGGGCGGGGTTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                                                      <  2:5050563/145‑1 (MQ=255)
aTGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                                                      <  2:4887262/150‑1 (MQ=255)
 tgtgtgGCGGGACGCCCCGGAAAATCCTGAACCGCCGCGCGGGACTAAAATGTATTCGTATGGTTGGTTTTCTCCCTCCCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                                                     <  2:1710973/150‑1 (MQ=25)
 tgtggggccggccGCCCCGGAAAATACTGCCCCGCGGCGCCGTACTAAAATTTTTTGGTTTGGTTGTTTTTCACCCACTCCGTGGGGGTGTGGGGGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                                                     <  1:5851939/145‑1 (MQ=255)
 tGTCTGGCGGGACGACCCGGCAAATCCGGCACCGCCGCGACGGACTAAAATTTTTTCGTTTGGTTGGTTTTCCCCCTCTCCGTGTGGGTTTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                                                     <  2:1711122/150‑1 (MQ=255)
     tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACAACCAAACCCACTCCCACAAggggg                                                                                                                                                 >  2:4505283/1‑148 (MQ=38)
      ggCGGGACGCCCCGCCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGTTTGGTTTTCGCCCTCACCGGGGGGGCGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                                                                <  2:4156635/150‑1 (MQ=255)
      ggCCGGACGACCCGGCAATTACTGCACGGACGCGACGGACGAAAAAGTTTTGGTAAGGTTGGTTTCCTCCCTCGCCCGGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                                                                <  2:4814256/150‑1 (MQ=255)
       gcggtacgccCCTGAAAATACGGCACGGCCGCGCCGGACGAAAATTTTTTCTTTTGTTTGTTTTTCACCCTCTCCGTGTGGTTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                                                               <  1:3684220/145‑1 (MQ=255)
       gCCGGACGACCCGGCAAATACTGCACCGACGCGACGTACGAAATAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                                               <  2:4864509/150‑1 (MQ=255)
       gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGCCCCGCACCCCAAACCACGCCCGCAAggggggg                                                                                                                                               >  1:1478369/1‑146 (MQ=37)
       cccGGACCCCCCGGAAAATACTGCACCGCCGCGACGCACGAAAAATTATTGTTACGGTTGTTTTTCGCCCACCCCCGGTGGGCGAGGGCCGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                                                               <  2:2216801/149‑1 (MQ=255)
        ccGGACCACCGGGAAAATACTGCACCCACGCGAGGGACCAAAATGTTTGCGTACGGTGGGTTTTCACCCACCCCCTGTGGGAGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                                                              <  2:1597878/150‑1 (MQ=255)
         cgtacgacCCCTAAAATATTCCACCTACGCGCCTCACTAAAATTTTTTCTTATGTTTGTTTTTGCCCCTCTCCTTGTGGTTGTGGGCGGGGGTTTGGGCATCAGCGCGCACGTTCACACTCACCCTAACACTCTCACTCAAGGGAGAggg                                                                                                                                             <  2:5054183/147‑1 (MQ=255)
         cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCACACCCAAAACCACACCCCCCCGGGGGAggg                                                                                                                                             >  2:1340272/1‑150 (MQ=25)
           gacgccCCGGAAAATACGGCCCCGACGCGCCGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCGCGCCCTGTGGTAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                                                           <  2:729266/150‑1 (MQ=38)
           gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCGAACCCACACCCGCACGGGGGAggggg                                                                                                                                           >  2:5470541/1‑149 (MQ=37)
             ccacCCGGAAATACCTGCACCGCCGCGACTCACTAAAATTTATTCATAGGGTTGTTTTTCCCCCTCGCCGTGTGGGAGTGGGCCGGGGTGAGGGACTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                                                                         <  2:24796/148‑1 (MQ=255)
                  cGGCAAATACTGCACCGACGCGACGGACGACAAAGTATTCGTCAGGGTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCACCCGCACCCCAACCCGCGCCCCCCCGGGCGCGGGGAccccgcg                                                                                                                                    >  2:2491752/1‑145 (MQ=21)
                  cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCGCACCCAAACCCACACCCACAACGGAGCGGGGGCCGAACg                                                                                                                                    >  1:5740558/1‑150 (MQ=25)
                  cGGAAAATACTGCACCGCCGCGCCGGAATAAAATTTTTTCTTATGTTTGTTTTTCTCCCTCCCCCTGTGGGAGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCg                                                                                                                                    <  1:5887380/150‑1 (MQ=255)
                    gCAAATACTGCCCCGACGCGACGGACCAAAAAGTTTTCTTAAGGTGGGTTTTCACCCTCTCCCTGGGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                                                                  <  1:4356040/150‑1 (MQ=255)
                    gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                                                                  <  2:4951167/150‑1 (MQ=255)
                     aaaaTCCTCCACCGACGCGACGGACTAAAATTTATTCTTCTGTTTGGTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGc                                                                                                                                 <  1:2872994/149‑1 (MQ=255)
                      aaaTAGTGCACCGCCGCTACGCACGAAAATTTTTGCGTCTGTTTGTTTTTCGCCCTCGCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                                                                <  2:200292/150‑1 (MQ=255)
                              gCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCACACCCAAAACCAAAACCACAAAGGGGAGGGGGCCCCACGAGcacacaaaa                                                                                                                         >  2:4702842/1‑145 (MQ=21)
                               cACCGCCGCGACGGACGAAAATGTTTTCTTATGGTTGGTTTTCGCCCTCTCCGGGTGGGAGTGGGCGGGGGGGTGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                                                        <  1:1362674/150‑1 (MQ=255)
                               cACCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCCCCCTCTCCCTGTGGGAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                                                        <  1:853925/150‑1 (MQ=255)
                                acggccgcgacgGACGAAAATTTATAGGTATGGTGGGTTTCCCCCCGCTCCGGGGGGGATCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                                                                       <  1:3263290/145‑1 (MQ=255)
                                 ccGACGCGACGGACGAAAAATTATTCTTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGTGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                      <  2:1459364/150‑1 (MQ=38)
                                 ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCGGGGGGCGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                      <  2:3654599/150‑1 (MQ=255)
                                           ggacgAAAATGTATTCGTATGTTTGGTTTCCGCCCTCTCCCTGTGGTAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGTGACCGATCGAGCACAAATTTTGAGAGATAtctt                                                                                                            <  2:4235759/150‑1 (MQ=38)
                                           ggacgAAAAAGTATGCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtctt                                                                                                            <  2:3395320/150‑1 (MQ=255)
                                               gAACAAGTATTCGTACGGGTGGTTTTCTCCCTCGCCCTGGGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                                                                        <  2:1967817/150‑1 (MQ=255)
                                                aaaacGTATTGGTAAGGTTGTTTTTCCCCCTCTCCGTGTGGGAGAGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                                                                       <  2:2298686/150‑1 (MQ=255)
                                                    cgtcttCGTACGGGTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACTAATTTTGAGCGATATCTTCTTCTGTCt                                                                                                   <  2:5422791/146‑1 (MQ=38)
                                                     gTATTGGTTGGGTTGGTGTTCTCCCTCGCCGGGGGGGTGTGGGCGGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                                  <  2:609880/150‑1 (MQ=255)
                                                         tggTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCTCCCTCACGCTAACCCTCTCCCTGACGGGAGAGGGGCCCGATGGAGCACAACTTTTGAGAGATCTCGTCTTCTGTGTTGTaa                                                                                              <  2:3976378/148‑1 (MQ=25)
                                                         tCGTATGGTTGGTTTTCTCCCTCTCCCTGGGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                                                              <  1:5381416/150‑1 (MQ=255)
                                                                gttggttTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                                                                       <  1:4505439/150‑1 (MQ=255)
                                                                gttggttTTCTCCCTCTCCCTGTGGCTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                                                                       <  1:4504978/150‑1 (MQ=255)
                                                                   ggttTTCGCCATCTCCGGGTGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                                                                    <  1:2760826/150‑1 (MQ=255)
                                                                          cccctccccttGTGGTTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                                                             <  1:512413/149‑1 (MQ=255)
                                                                           ccctctccctGTGGGAGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTATTGTAACAGAAGAACAGAAAATCg                                                                            <  2:3707100/150‑1 (MQ=255)
                                                                            cctctccgtGTGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                                                           <  1:3277824/150‑1 (MQ=255)
                                                                             ctcgccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                                                                          <  2:3592649/150‑1 (MQ=255)
                                                                             ctccccctGTGGGAGAGGCCCGGGGTGAGGGCATCAGCGCGCACGTACACCCTCACCCTAACCCACTCCCTCAAGGGAGAGGGGACCGAACGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                                                                          <  2:5482463/150‑1 (MQ=255)
                                                                                 ccctGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                                                                      <  1:2847954/150‑1 (MQ=255)
                                                                                  cctGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                                                                     <  2:5817471/150‑1 (MQ=255)
                                                                                    gcggggAGCGGGCGGGGGTGCGGCCCTCAGCGCGCACGTGCACCCGCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGCTCGAGCCCAAATTGTGAGAGAGATCTACTTCTGTCTGGTAACAGAAGAACCGAAACGCGGGCTTTTgc                                                                   <  2:4108696/147‑1 (MQ=255)
                                                                                      tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAACCCACACCCGCAAGGGGGTGGGGGCCCATCGAGCACAAAATTTTAGATATATATTAATATGTCTTTTAACAGAAGAAAAGAAAAAAGGGCGTTTtcgc                                                                 >  2:3007351/1‑150 (MQ=255)
                                                                                      ggggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                                                                 <  1:1386496/149‑1 (MQ=255)
                                                                                      ggggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                                                                 <  2:2894026/149‑1 (MQ=255)
                                                                                        ggAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                                               <  2:2375443/150‑1 (MQ=255)
                                                                                        ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCTCACCCAAAACCAAGACCAAAAGGGGGGGGGGGCCCGGCGCGGACAAAATTTGGGAGAAATGTTGTTGTGTGTTTTTACAGAAGAACAGAAAATCGGGCgttttggccg                                                               >  2:3368834/1‑144 (MQ=255)
                                                                                          agagGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCAAACCCACGCCCGCAAGGGAGAGGGGACCCATCGAGCACAAAATTTGAGAGATATCTTCTTATGTCTTGTAACAGAAGAACACAAAAAAGGGCGTTTGCGCCGGa                                                             >  2:3785982/1‑150 (MQ=255)
                                                                                             gggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAata                                                          <  2:4376385/150‑1 (MQ=255)
                                                                                                 cGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                                      <  2:1478369/150‑1 (MQ=255)
                                                                                                  ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                     <  1:1950152/150‑1 (MQ=255)
                                                                                                        aGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCACTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaa                                               <  2:4799987/150‑1 (MQ=255)
                                                                                                         gggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCGGAATATTGCGGAGaaaa                                              <  2:2431823/150‑1 (MQ=255)
                                                                                                             aTCAGCGCGCACGTTCACCCACACCCAAAACCACAACCAAAAGGGAGAGGGGACCCATCGAGCACAAAATTTGAGAGAAATAATATTATGTCTTGTAAAAGAAGAACAGAAAAAAGGGCGTTTTCGCCGGAATAATTCGGAGAAAAAACa                                          >  2:3498292/1‑150 (MQ=255)
                                                                                                              tCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAACGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCAt                                         <  2:5250969/150‑1 (MQ=255)
                                                                                                                   gcgcACGTTCACCCTCACCCTAACCCTCTCCCTCTAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCg                                    <  1:3694863/150‑1 (MQ=255)
                                                                                                                    cgcACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg                                   <  1:4132333/150‑1 (MQ=255)
                                                                                                                     gcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCggg                                  <  2:3295111/150‑1 (MQ=255)
                                                                                                                             accctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcg                          <  1:3502273/150‑1 (MQ=255)
                                                                                                                                ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                       >  1:1811995/1‑150 (MQ=255)
                                                                                                                                ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                       >  2:317989/1‑150 (MQ=255)
                                                                                                                                 tcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAg                      >  2:1751472/1‑150 (MQ=255)
                                                                                                                                  caccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGc                     <  1:3544596/150‑1 (MQ=255)
                                                                                                                                  caccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGc                     <  2:5702485/150‑1 (MQ=255)
                                                                                                                                  caccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGc                     <  1:3363177/150‑1 (MQ=255)
                                                                                                                                  caccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGc                     <  1:5665603/150‑1 (MQ=255)
                                                                                                                                  caccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGc                     <  2:5702380/150‑1 (MQ=255)
                                                                                                                                     cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                  >  1:1778293/1‑150 (MQ=255)
                                                                                                                                     cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                  >  1:2200576/1‑150 (MQ=255)
                                                                                                                                     cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGAGCCTGAATATTGCGGAGAAAAAGAATGGACGGGAATTACCGCGAGCCaa                  >  2:3562515/1‑150 (MQ=255)
                                                                                                                                     cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAGTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGAGccca                  >  2:4560853/1‑148 (MQ=255)
                                                                                                                                       tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt                <  1:4845125/150‑1 (MQ=255)
                                                                                                                                       tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt                <  1:4090540/150‑1 (MQ=255)
                                                                                                                                       tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt                <  1:25518/150‑1 (MQ=255)
                                                                                                                                        tccccTCTCCCTTGAGGGAGAGGGTACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg               >  1:180958/3‑150 (MQ=255)
                                                                                                                                        tccccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGAGCCAATTg               >  2:854584/3‑150 (MQ=255)
                                                                                                                                        tccccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg               >  2:445923/3‑150 (MQ=255)
                                                                                                                                          cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa             <  2:4445202/150‑1 (MQ=255)
                                                                                                                                          cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa             <  2:2574520/150‑1 (MQ=255)
                                                                                                                                          cccTCTCCCGTGAGGGAGAGGGGACCGATCGAGCCCAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGACAAAGCATGGCCGGGAATGACCGCGAGCCACTTGGa             >  2:3572010/1‑150 (MQ=255)
                                                                                                                                           ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc            >  1:3902196/1‑150 (MQ=255)
                                                                                                                                           ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc            <  1:4741143/150‑1 (MQ=255)
                                                                                                                                           ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAGGTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATAGCGGAGAAAAAGCATGGCCGGGAATGACCGCTAGCCAATTGGAc            >  2:5568521/1‑150 (MQ=255)
                                                                                                                                            ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg           <  1:1337007/150‑1 (MQ=255)
                                                                                                                                            ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg           <  1:1337003/150‑1 (MQ=255)
                                                                                                                                             tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc          >  1:421002/1‑150 (MQ=255)
                                                                                                                                               acccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTAGAGAGATATCTTCTTCTGTCTTGAAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCaa        <  2:710523/149‑1 (MQ=255)
                                                                                                                                                  cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg     <  2:4982623/150‑1 (MQ=255)
                                                                                                                                                  cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg     >  2:2313494/1‑150 (MQ=255)
                                                                                                                                                   tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACACAAAggg    >  1:189001/1‑150 (MQ=255)
                                                                                                                                                    cAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGa   <  2:3452718/150‑1 (MQ=255)
                                                                                                                                                    cAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGa   >  1:656724/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAA  >  NC_000913/164473‑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: