Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,620 C→T 6.6% intergenic (+86/‑110) 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,6200CT 6.6% 870.7 / 32.9 273intergenic (+86/‑110)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (83/170);  new base T (8/10);  total (92/181)
Fisher's exact test for biased strand distribution p-value = 3.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 275 total aligned reads displayed.

ATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGA  >  NC_000913/164473‑164769
                                                                                                                                                   |                                                                                                                                                     
aTGTCTGGCCGGACCACCCGGAAAATACGGCCCCGACGCGACGGACTAAAATTTTTTGGTATGGTTGTTTTTCTCCCACCCCGTGGGGGATTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                                                     <  2:3175311/150‑1 (MQ=255)
 ggtggggCCGGACGACCGGGCACATACGGCCCCGCCGCGCCGGCCGAAAAAGTATTGGTATGGTGGGTTTTCTCCCTCGCCCGGGGGGAGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                                                    <  2:293317/145‑1 (MQ=255)
  gTCTGGCCGGACGACCCGGCAAAAACTGCACCGACGCGACGTACCAAAAATTATTCTTAAGGTTGGTTTTCCCCCTCTCCCTGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                                                                   <  2:1789872/150‑1 (MQ=38)
   tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGGGCACCCGCACCCAAACCCACACCCGCCAggg                                                                                                                                                  >  2:8171921/1‑150 (MQ=37)
      ggCCGCCCGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATGCGTACGGTTGGTTTTCTCCCCCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGgaga                                                                                                                                               <  1:6103928/150‑1 (MQ=38)
       gCCGGACGACCCGCCAAATCCGGCACCGACGCGACGGACGAAAAAGTATGCGTATGGTGGGTTTTCTCCCTCACCCTGTGGGTGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                                              <  1:8389784/150‑1 (MQ=255)
        ccGGACGCCCCGGCAAATACTGCACCGACGCGCCGGACGAAAACGTATTCGTAAGGTTGGTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                                                             <  2:5463411/150‑1 (MQ=38)
          ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTGAACCACACCCAAACCCACACCCCCCACGGGGAgggg                                                                                                                                           >  1:1365179/1‑150 (MQ=37)
           gacgacCCGGAAAATACGGCACCGCGGCGCCGGACTAAAATTTTTTCTTATGGTGGTTTTTCCCCCCCGCCCGGTGGGGGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                                                          <  1:1957795/150‑1 (MQ=255)
              gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAACCCGCTCCCACCAGGGAGAGGGGCCCg                                                                                                                                       >  2:4881150/1‑150 (MQ=38)
                cccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCCCACCCAAACCCCCACCCCCCAGGGAGGGGGGGCCgcg                                                                                                                                     >  1:6356193/1‑148 (MQ=25)
                  cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCACACCCCAACCCACTCCCACAAGGGCGGGGGGGCCCATCg                                                                                                                                   >  2:7342754/1‑150 (MQ=25)
                  cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCCCCCAAAACCCCAACCACAAGGGGGGGGGGGccccgag                                                                                                                                   >  2:1641224/1‑145 (MQ=21)
                           cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAAACACAACCACAAGGGAGATGGGACCCATCGAGAACAAAt                                                                                                                          >  2:7734226/1‑150 (MQ=25)
                            tGCACCGACGCTACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtt                                                                                                                         <  2:4762457/150‑1 (MQ=255)
                              cACCGCCGCGCCGGACTAAAATTTTTTCTTTTGTTTGGTTTTCGCCCGGGCGGGGGGGGTGGGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                                                       <  2:2862715/150‑1 (MQ=255)
                               cccGACGCGAGGGACCAAAATTTATTGGTAAGGTTGTTTTTCTCCCGCACCCTGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                                                                      <  1:2605190/149‑1 (MQ=255)
                                ccgccgCGACGGACTAAAAAGTATGCGTAAGGTGGGTGTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                     <  1:8199539/150‑1 (MQ=38)
                                ccgccgCGACGGACGAAAATGAATTCTTACGGTTGGTTTTCGCCCTCGCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                     <  2:6683029/150‑1 (MQ=255)
                                          ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCCCACCCAAAACCACACCCCCAACGGGGGGGGGGCCGAACGCGCCCAAAATTTTTGAGATAtatt                                                                                                           >  1:974418/1‑150 (MQ=18)
                                          ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCGCCCCAAACCACAACCACAAAGGAGAGGGGGCCGCGAGAGCACCAAATTTgtgaggtttaat                                                                                                           >  1:3780794/1‑143 (MQ=16)
                                          ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTCCACCCGCACCCAAAAACACAACCAAAAAGGGGGGGGGGCCCATCGAGCCCAAAATTTTATAGATAtatt                                                                                                           >  1:4267577/1‑150 (MQ=16)
                                                 aatttATTCTTATGTTTGTTTTCCACCCTCCCCCTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                                                                    <  2:3737756/146‑1 (MQ=255)
                                                    gTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                                 <  2:6982084/150‑1 (MQ=255)
                                                     tagccGTAAGGTTGGTTTTCGCCCGCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                                                                <  2:4490259/146‑1 (MQ=255)
                                                          gTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                                                           <  2:7210112/150‑1 (MQ=255)
                                                          gTAAGGTTGGTTTTCGCCCCCTCCCCGTGGGTGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                                                           <  2:4934331/150‑1 (MQ=255)
                                                            atggttggttTTCCCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                                                         <  2:615431/148‑1 (MQ=255)
                                                            aaGGTTGTTTTTCTCCCTCTCCCTGTGGTCTTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                                                         <  2:5310313/150‑1 (MQ=255)
                                                                ttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGAACCGCACCCAAAACCACACCCCCCAGGGGGTGGGGCCCGAGCGAGCCCAAAAGTTGTGAGATAAGTTATTGTGTGTTGTTAGAgaagaa                                                                                     >  2:850550/1‑150 (MQ=12)
                                                                ttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAACCACTACCACAAGCGAGGGGGGGCCCAGCGCGCGCAAATTTTGTGAGGTATATTATTTTGTGGTGTTAAagaacta                                                                                     >  1:6907370/1‑147 (MQ=12)
                                                                ttggttTTCTCCCGCTCCCTGGGGGCTTGGGCCGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                                                                     <  1:1869658/150‑1 (MQ=255)
                                                                 ggtttttCTCCCTCGCCCGGTGGGTGCGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                                                                    <  1:8058354/147‑1 (MQ=255)
                                                                    ttttCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCCAACCCACTCCCACAAGGGGGGGGGGACCGCAAGAGCCCAAATTTTGTGGGGTATCTTCATCTGTCTTGTTACAGAAGAACAGa                                                                                 >  2:2118352/1‑150 (MQ=21)
                                                                          ccctctcccgGGGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                                                           <  1:729024/150‑1 (MQ=255)
                                                                                ccctGTGGGCGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                                                                     <  1:3516318/150‑1 (MQ=255)
                                                                                  cggggggAGCGGGCCGGGGTGCGGCCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                                                                   <  1:2315442/146‑1 (MQ=255)
                                                                                     ggggAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                                                                <  2:5093514/149‑1 (MQ=255)
                                                                                      gtgtgagGGCGGGGGTGTGGGAATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                                                               <  2:3607970/146‑1 (MQ=255)
                                                                                         agagGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCCCCCAAAACCGCACCCACAAGGGGGGGGGGTCCCGTAGAGCACAAAATTTGTGGGATATAATATTATGTCTTGTAACACAAAAACAAAAAAAAGGGCGGTTGGGCCAGa                                                            >  2:1372014/1‑150 (MQ=255)
                                                                                              gCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt                                                       <  1:7915150/150‑1 (MQ=255)
                                                                                               ccGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGATCAGAAAATCGGGCTTTTGCGCCTGAATATTg                                                      <  2:6407012/150‑1 (MQ=255)
                                                                                                cGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                                     <  1:7734226/150‑1 (MQ=255)
                                                                                                 ggggTGAGGGCCTCAGAGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                    <  1:629969/150‑1 (MQ=255)
                                                                                                 ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                    <  2:4537206/150‑1 (MQ=255)
                                                                                                 ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                    <  2:4663731/150‑1 (MQ=255)
                                                                                                   ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga                                                  <  2:5246802/150‑1 (MQ=255)
                                                                                                       cggccctCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGTGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGTAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaa                                              <  1:6522288/144‑1 (MQ=255)
                                                                                                         ggCCTCCGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa                                            <  1:2433577/150‑1 (MQ=255)
                                                                                                         gccctccgcgcgcCCGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa                                            <  1:1711066/143‑1 (MQ=255)
                                                                                                                  gcgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCg                                   <  2:6436580/150‑1 (MQ=255)
                                                                                                                   cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg                                  <  1:5937141/150‑1 (MQ=255)
                                                                                                                   cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg                                  <  2:8246550/150‑1 (MQ=255)
                                                                                                                   cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg                                  <  1:6792897/150‑1 (MQ=255)
                                                                                                                     cACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGa                                <  2:5724312/150‑1 (MQ=255)
                                                                                                                     cACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATAGCGGAGAAAAAGCATGGCCGGGa                                <  1:3634479/150‑1 (MQ=255)
                                                                                                                      gctgatgccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGaa                               >  2:1862559/8‑150 (MQ=255)
                                                                                                                       cGTTCACCCTCACCCTAACCCTCTCACTCAAGTGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTTCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAAt                              <  1:557666/150‑1 (MQ=255)
                                                                                                                         ttCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGa                            <  2:7455417/150‑1 (MQ=255)
                                                                                                                          tcaccctcacccTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGAc                           <  1:4702539/150‑1 (MQ=255)
                                                                                                                          tcacccaaacccAAACCCAATCCCTCAATGGAGAGGGGACCGATAGAGCACAAATTTTGAGAGATATATTAATATGTCTTTTAACAGAATAACAGAAAAAAGGGCTTTTTCGCCTGAATATTTCGGAGAAAAATCAAGGCCGGGAAAGaa                           >  2:4164155/1‑149 (MQ=255)
                                                                                                                           caccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGAcc                          <  1:6592533/150‑1 (MQ=255)
                                                                                                                            accctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCCCAAATTTTGAGAGATATCTTCTTCGGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcg                         <  2:2061777/150‑1 (MQ=255)
                                                                                                                            accctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcg                         <  2:4982242/150‑1 (MQ=255)
                                                                                                                               ctcaccctAACCCTCTCCCTCACGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                      <  1:4082529/150‑1 (MQ=255)
                                                                                                                               ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                      >  1:1820740/1‑150 (MQ=255)
                                                                                                                                tcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAg                     <  2:7362430/150‑1 (MQ=255)
                                                                                                                                tcacccgAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAg                     <  1:2899574/150‑1 (MQ=255)
                                                                                                                                    cctAACCCTCTCCCTCCAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                 >  2:3816068/1‑150 (MQ=255)
                                                                                                                                    cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGAGCCaa                 >  1:7130453/1‑150 (MQ=255)
                                                                                                                                    cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGAGCCaa                 >  2:5978989/1‑150 (MQ=255)
                                                                                                                                    cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                 >  1:7031386/1‑150 (MQ=255)
                                                                                                                                    cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                 <  1:245186/150‑1 (MQ=255)
                                                                                                                                    cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAAGGACCGCGAGccca                 >  2:4423167/1‑148 (MQ=255)
                                                                                                                                     ctACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAt                <  1:484327/150‑1 (MQ=255)
                                                                                                                                      tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt               >  1:6694845/1‑150 (MQ=255)
                                                                                                                                      tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt               >  2:2437010/1‑150 (MQ=255)
                                                                                                                                       aaCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg              <  2:4267577/150‑1 (MQ=255)
                                                                                                                                       aaCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg              <  2:6797765/150‑1 (MQ=255)
                                                                                                                                        aCCCTCTCCCTCCAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTgg             <  1:2648461/150‑1 (MQ=255)
                                                                                                                                        aCCCTCTCCCTCATGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTgg             >  2:1266717/1‑150 (MQ=255)
                                                                                                                                        aCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTgg             <  1:5109162/150‑1 (MQ=255)
                                                                                                                                        aCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGAGATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTgg             <  1:3994341/150‑1 (MQ=255)
                                                                                                                                         cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa            <  1:1302563/150‑1 (MQ=255)
                                                                                                                                         cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa            <  1:4444107/150‑1 (MQ=255)
                                                                                                                                          ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc           >  2:7981085/1‑150 (MQ=255)
                                                                                                                                          ccTCTCCCGTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc           >  1:3139035/1‑150 (MQ=255)
                                                                                                                                            tctcCCTTGAGGGAGAGGGTACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc         >  1:6184402/1‑150 (MQ=255)
                                                                                                                                            tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc         >  2:7899252/1‑150 (MQ=255)
                                                                                                                                            tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc         >  2:2637432/1‑150 (MQ=255)
                                                                                                                                                ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGCGCCAATTGGACGCAAAg     >  1:3150729/1‑150 (MQ=255)
                                                                                                                                                ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAg     >  1:1164766/1‑150 (MQ=255)
                                                                                                                                                ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAg     <  2:5672403/150‑1 (MQ=255)
                                                                                                                                                ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAg     >  1:7501311/1‑150 (MQ=255)
                                                                                                                                                ccTCAAGGGAGAGGGGACCGATCGAGCACAAATATTGAGAGATATATTCTTATGTCAAGTAACAGAATAAAAGAAAATCGGGCTTTTCCGAATGAAAATTGCGGAGAAAAAGCATGTCCGGGAATGACAGAGAGCCAATAGGACGAAAAg     >  1:2118658/1‑150 (MQ=255)
                                                                                                                                                 cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg    >  1:459308/1‑150 (MQ=255)
                                                                                                                                                 cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg    <  1:8423366/150‑1 (MQ=255)
                                                                                                                                                 cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTATGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg    >  1:5061975/1‑150 (MQ=255)
                                                                                                                                                  tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                                                        >  1:5784786/1‑65 (MQ=255)
                                                                                                                                                   cAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGa  <  2:6832231/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGA  >  NC_000913/164473‑164769

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: