Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,621 A→G 6.7% 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 6.7% 826.1 / 35.8 274intergenic (+87/‑109)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (77/170);  new base G (8/10);  total (93/181)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-01

Only 100 of 276 total aligned reads displayed.

ATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAA  >  NC_000913/164473‑164770
                                                                                                                                                    |                                                                                                                                                     
 tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACACACACCCAAACCCACTCCCACAAg                                                                                                                                                     >  1:6300272/1‑150 (MQ=38)
   tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGGGCACCCGCACCCAAACCCACACCCGCCAggg                                                                                                                                                   >  2:8171921/1‑150 (MQ=37)
      ggCCGCCCGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATGCGTACGGTTGGTTTTCTCCCCCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGgaga                                                                                                                                                <  1:6103928/150‑1 (MQ=38)
       gCCGGACGACCCGGAAACTCCTGCACCGACGCTACGCACGAAAATTTATTCGTATGGTTGTTTTACCCCCTCCCCCTGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                                                               <  1:4873490/150‑1 (MQ=255)
       gCCGGACGACCCGCCAAATCCGGCACCGACGCGACGGACGAAAAAGTATGCGTATGGTGGGTTTTCTCCCTCACCCTGTGGGTGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                                               <  1:8389784/150‑1 (MQ=255)
       cccGGACGACCCGGCCCCTACTGCACCCACGCGCCGGACGACAACGTATGCGTATGGTTGGTTTTCCCCCCCCCCCGGGGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                                               <  2:795209/149‑1 (MQ=255)
        ccGGACGCCCCGGCAAATCCTGCCCCGCCGCGACGGACGAAAACTTTTTCGTAAGGTTGGTTTCCCCCCTCGCCCGGGGGGAGAGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                                                              <  2:6524299/150‑1 (MQ=255)
          ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTGAACCACACCCAAACCCACACCCCCCACGGGGAgggg                                                                                                                                            >  1:1365179/1‑150 (MQ=37)
           gccgccccGGCAAACCCTGCCCCGCCGCGACGGCCTAAACAGTATGCGTCTGGTTGGTTTTCGCCCGCCCCCGGGGGGCGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAGCCCTCTCCCTCAAGGGAGAGGGGa                                                                                                                                           <  2:7139958/145‑1 (MQ=255)
           gacgacCCGGAAAATACGGCACCGCGGCGCCGGACTAAAATTTTTTCTTATGGTGGTTTTTCCCCCCCGCCCGGTGGGGGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                                                           <  1:1957795/150‑1 (MQ=255)
             cgccccGGCAAATACTGCACCGACGCGCCGGCCTAAAAAGTAGGCGTATGGTTGTTTTCCGCCCCCGCCGGGGGGGTGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttag                                                                                                                                         <  2:7210627/147‑5 (MQ=255)
                cccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCCCACCCAAACCCCCACCCCCCAGGGAGGGGGGGCCgcg                                                                                                                                      >  1:6356193/1‑148 (MQ=25)
                  cGGCAAATACTGCCCCGACGCGACGGACGAAAACGAATTCGTATGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCg                                                                                                                                    <  2:4508326/150‑1 (MQ=38)
                  cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCACACCCCAACCCACTCCCACAAGGGCGGGGGGGCCCATCg                                                                                                                                    >  2:7342754/1‑150 (MQ=25)
                     cAAATCCTGCCCCGACGCGCCGGCCTAAAAAGTATTCGTAAGGTTGGTTTTCGCCCTCCCCCTGGGGGTGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGc                                                                                                                                 <  2:6810620/150‑1 (MQ=255)
                           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)
                                  gacgcgacgGACGAAAAAGTATTGTTAAGGTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                                                                    <  2:6695220/150‑1 (MQ=39)
                                          ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCGCCCCAAACCACAACCACAAAGGAGAGGGGGCCGCGAGAGCACCAAATTTgtgaggtttaat                                                                                                            >  1:3780794/1‑143 (MQ=16)
                                          ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTCCACCCGCACCCAAAAACACAACCAAAAAGGGGGGGGGGCCCATCGAGCCCAAAATTTTATAGATAtatt                                                                                                            >  1:4267577/1‑150 (MQ=16)
                                                 aatttATTCTTATGTTTGTTTTCCACCCTCCCCCTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                                                                     <  2:3737756/146‑1 (MQ=255)
                                                    gTATTCGTATGGTTGGTTTTCGCCCTCGCCCTGTGGGAGCGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                                  <  2:5025023/150‑1 (MQ=255)
                                                     tATTGGTATGGTGGGTTTACGCCCTCTCCCTGTGGTTGCGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                                                                 <  2:1761942/150‑1 (MQ=255)
                                                       ttCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGCGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa                                                                                               <  2:4568610/150‑1 (MQ=255)
                                                        tCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                                                              <  2:6305760/150‑1 (MQ=255)
                                                          gTAAGGTTGGTTTTCTCCCTCTCCCCGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                                                            <  2:4934525/150‑1 (MQ=255)
                                                           tttggggggttTTGCCCCGCGCCCGGGGGGAGCGGGCCGGGGTGCGGGCCTCAGCGCGCCCGGTCCCCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACag                                                                                           <  2:5311753/143‑1 (MQ=37)
                                                             tggttgtttTTCTCCCTCGCCGTGTGGGCTTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                                                         <  2:3199411/149‑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)
                                                                  ggttTTCGCCCGCGCCCGGTGGGCGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                                                                    <  2:1830336/150‑1 (MQ=255)
                                                                         tccctctccctGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                                                             <  2:5490067/150‑1 (MQ=255)
                                                                          ccctctcccgGGGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                                                            <  1:729024/150‑1 (MQ=255)
                                                                           cctctccctGTGGTAGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                                                           <  2:115883/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)
                                                                                               ccGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGATCAGAAAATCGGGCTTTTGCGCCTGAATATTg                                                       <  2:6407012/150‑1 (MQ=255)
                                                                                               ccGGGGTGAGGGCCTCAGCGCGCACGTGCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGAGATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg                                                       <  1:7972200/150‑1 (MQ=255)
                                                                                                cGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                                      <  1:7734226/150‑1 (MQ=255)
                                                                                                 ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                     <  2:4663731/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)
                                                                                                           cATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc                                           <  2:419167/150‑1 (MQ=255)
                                                                                                           cATCAGCGCGCACGTTCACACTCACCCTATCCCTCTGCCTAAAGGGAGAGGGGACCGATCGAGCACTAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGTAATAGc                                           <  2:6026056/150‑1 (MQ=255)
                                                                                                                  gcgcCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCg                                    <  1:7135257/150‑1 (MQ=255)
                                                                                                                   cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg                                   <  1:6792897/150‑1 (MQ=255)
                                                                                                                   cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg                                   <  1:5937141/150‑1 (MQ=255)
                                                                                                                   cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg                                   <  2:8246550/150‑1 (MQ=255)
                                                                                                                     cACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGa                                 <  2:5724611/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)
                                                                                                                           caccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGAcc                           <  1:6592533/150‑1 (MQ=255)
                                                                                                                            accctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCCCAAATTTTGAGAGATATCTTCTTCGGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcg                          <  2:2061777/150‑1 (MQ=255)
                                                                                                                              cctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcgcg                        >  1:5796894/1‑150 (MQ=255)
                                                                                                                               ctcaccctAACCCTCTCCCTCACGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                       <  1:4082529/150‑1 (MQ=255)
                                                                                                                               ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                       >  1:1820740/1‑150 (MQ=255)
                                                                                                                                tcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAg                      <  2:2119520/150‑1 (MQ=255)
                                                                                                                                tcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAg                      <  2:7362107/150‑1 (MQ=255)
                                                                                                                                tcacccgAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAg                      <  1:2899574/150‑1 (MQ=255)
                                                                                                                                   ccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCa                   >  2:2453560/1‑150 (MQ=255)
                                                                                                                                    cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                  <  1:245186/150‑1 (MQ=255)
                                                                                                                                    cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                  >  1:7031386/1‑150 (MQ=255)
                                                                                                                                    cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAAGGACCGCGAGccca                  >  2:4423167/1‑148 (MQ=255)
                                                                                                                                       tccccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg               >  2:390097/3‑150 (MQ=255)
                                                                                                                                       aaCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg               <  2:4267577/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:2979721/150‑1 (MQ=255)
                                                                                                                                          ccTCTCCCGTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc            >  1:3139035/1‑150 (MQ=255)
                                                                                                                                           ctctCCCTCAAGGGAGAGGGGACCGATCGCGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAACATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg           <  2:5290799/150‑1 (MQ=255)
                                                                                                                                           ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAGTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg           <  1:6719833/150‑1 (MQ=255)
                                                                                                                                            tctcCCTTGAGGGAGAGGGTACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc          >  1:6184402/1‑150 (MQ=255)
                                                                                                                                            tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc          >  2:2637432/1‑150 (MQ=255)
                                                                                                                                            tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc          >  2:7899252/1‑150 (MQ=255)
                                                                                                                                              tcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCaa        <  2:4913708/150‑1 (MQ=255)
                                                                                                                                               cccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCaaa       <  1:952595/150‑1 (MQ=255)
                                                                                                                                                ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGCGCCAATTGGACGCAAAg      >  1:3150729/1‑150 (MQ=255)
                                                                                                                                                ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAg      <  2:5672403/150‑1 (MQ=255)
                                                                                                                                                ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAg      >  1:1164766/1‑150 (MQ=255)
                                                                                                                                                ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAg      >  1:7501311/1‑150 (MQ=255)
                                                                                                                                                ccTCAAGGGAGAGGGGACCGATCGAGCACAAATATTGAGAGATATATTCTTATGTCAAGTAACAGAATAAAAGAAAATCGGGCTTTTCCGAATGAAAATTGCGGAGAAAAAGCATGTCCGGGAATGACAGAGAGCCAATAGGACGAAAAg      >  1:2118658/1‑150 (MQ=255)
                                                                                                                                                 cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg     >  1:459308/1‑150 (MQ=255)
                                                                                                                                                 cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTATGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg     >  1:5061975/1‑150 (MQ=255)
                                                                                                                                                   cAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGa   <  2:6832231/150‑1 (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: