Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,099,026 A→C 6.1% intergenic (+70/+28) tpm → / ← htpX thiopurine S‑methyltransferase/protease HtpX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,099,0260AC 6.1% 1087.3 / 41.1 319intergenic (+70/+28)tpm/htpXthiopurine S‑methyltransferase/protease HtpX
Reads supporting (aligned to +/- strand):  ref base A (162/137);  new base C (10/10);  total (172/147)
Fisher's exact test for biased strand distribution p-value = 8.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 327 total aligned reads displayed.

GTGATTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGCATGT  >  NC_002947/2098879‑2099175
                                                                                                                                                   |                                                                                                                                                     
      gAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGccc                                                                                                                                               >  1:5305071/1‑150 (MQ=255)
      gAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGccc                                                                                                                                               >  2:4074200/1‑150 (MQ=255)
       aGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCt                                                                                                                                              >  1:6974258/1‑150 (MQ=255)
         gACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCttt                                                                                                                                            <  1:7871758/150‑1 (MQ=255)
           cAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGt                                                                                                                                          <  1:1407531/150‑1 (MQ=255)
            aGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGtt                                                                                                                                         <  1:2673775/150‑1 (MQ=255)
             gCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAACAAAGGGCGACTTCGGTCGCCCTtttttt                                                                                                                                        <  2:7984792/150‑1 (MQ=255)
             gCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGttt                                                                                                                                        <  2:7610081/150‑1 (MQ=255)
              cccGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACTTCGGTCGCCCTtttttta                                                                                                                                       <  1:6444901/150‑2 (MQ=255)
                cGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTa                                                                                                                                     >  2:4635714/1‑150 (MQ=255)
                 gAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTAt                                                                                                                                    >  1:98090/1‑150 (MQ=255)
                 gAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTAt                                                                                                                                    >  2:3366502/1‑150 (MQ=255)
                     tttGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAgg                                                                                                                                <  2:6808537/150‑1 (MQ=255)
                      ttGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGt                                                                                                                               <  2:758227/150‑1 (MQ=255)
                       tgtgCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTg                                                                                                                              <  2:2345203/150‑1 (MQ=255)
                        gtgCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTgg                                                                                                                             >  1:8198726/1‑150 (MQ=255)
                         tgCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTACCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGa                                                                                                                            >  2:175402/1‑150 (MQ=255)
                           cAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATc                                                                                                                          >  1:2274543/1‑150 (MQ=255)
                            aGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCa                                                                                                                         >  2:3084825/1‑150 (MQ=255)
                             ggCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAg                                                                                                                        >  1:3810415/1‑150 (MQ=255)
                              gCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGc                                                                                                                       >  1:1823314/1‑150 (MQ=255)
                               ccGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGcc                                                                                                                      <  2:2417045/150‑1 (MQ=255)
                                cGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCcg                                                                                                                     >  2:1198197/1‑150 (MQ=255)
                                 ggggCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCcgc                                                                                                                    <  1:2925191/150‑1 (MQ=255)
                                    gCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCt                                                                                                                 <  1:5860193/150‑1 (MQ=255)
                                    gCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCt                                                                                                                 <  2:1351265/150‑1 (MQ=255)
                                    gCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAACCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCt                                                                                                                 <  2:5338328/150‑1 (MQ=255)
                                      gTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGa                                                                                                               <  2:7039938/150‑1 (MQ=255)
                                       tCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGAc                                                                                                              >  1:4303129/1‑150 (MQ=255)
                                       tCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGAc                                                                                                              >  1:4485453/1‑150 (MQ=255)
                                         gAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGc                                                                                                            >  2:1974876/1‑150 (MQ=255)
                                         gAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGc                                                                                                            >  2:5490744/1‑150 (MQ=255)
                                           gTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAg                                                                                                          <  1:3121253/150‑1 (MQ=255)
                                                   gAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTAGa                                                                                                  <  1:4224291/150‑1 (MQ=255)
                                                     gCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATg                                                                                                >  1:2505598/1‑150 (MQ=255)
                                                       gggTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCg                                                                                              <  2:8798473/150‑1 (MQ=255)
                                                        ggTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCgg                                                                                             <  2:6597742/150‑1 (MQ=255)
                                                             aTCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTc                                                                                        <  2:7337423/150‑1 (MQ=255)
                                                             aTCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTc                                                                                        <  1:4847229/150‑1 (MQ=255)
                                                             aTCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTc                                                                                        <  1:3989319/150‑1 (MQ=255)
                                                              tCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTcc                                                                                       >  1:5097022/1‑150 (MQ=255)
                                                              tCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTcc                                                                                       >  1:8641576/1‑150 (MQ=255)
                                                               cAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCa                                                                                      <  1:8278825/150‑1 (MQ=255)
                                                                 ggTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGc                                                                                    <  2:2129702/150‑1 (MQ=255)
                                                                 ggTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGc                                                                                    >  2:6707331/1‑150 (MQ=255)
                                                                 ggTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGc                                                                                    >  2:2961582/1‑150 (MQ=255)
                                                                 ggTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGc                                                                                    >  2:3670361/1‑150 (MQ=255)
                                                                    cAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggt                                                                                 <  1:1160316/150‑1 (MQ=255)
                                                                      gCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggtgg                                                                               <  1:7408973/150‑1 (MQ=255)
                                                                        cGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggtgggt                                                                             <  1:6101795/150‑1 (MQ=255)
                                                                         gATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggtgggtg                                                                            <  2:8547041/150‑1 (MQ=255)
                                                                          aTAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGgtgggtgg                                                                           >  1:3192674/1‑150 (MQ=255)
                                                                          aTAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGgtgggtgg                                                                           <  1:668932/150‑1 (MQ=255)
                                                                                aGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtcatc                                                                     <  2:4360185/150‑1 (MQ=255)
                                                                                aGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtcatc                                                                     >  2:2920231/1‑150 (MQ=255)
                                                                                aGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtcatc                                                                     <  1:1304001/150‑1 (MQ=255)
                                                                                 ggTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtcatca                                                                    <  2:2271084/150‑1 (MQ=255)
                                                                                      ggCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAAc                                                               <  2:3896725/150‑1 (MQ=255)
                                                                                             tGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCa                                                        <  1:3366502/150‑1 (MQ=255)
                                                                                               aCCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAgg                                                      <  1:4937544/150‑1 (MQ=255)
                                                                                               aCCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAgg                                                      >  2:6999337/1‑150 (MQ=255)
                                                                                                cccGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGc                                                     <  2:5658624/150‑1 (MQ=255)
                                                                                                 ccGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGcc                                                    >  1:3502310/1‑150 (MQ=255)
                                                                                                 ccGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGcc                                                    >  1:956739/1‑150 (MQ=255)
                                                                                                    cTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTg                                                 <  2:1440745/150‑1 (MQ=255)
                                                                                                     tCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGc                                                <  2:5176895/150‑1 (MQ=255)
                                                                                                       ccACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCtt                                              >  1:8063845/1‑150 (MQ=255)
                                                                                                        cacaGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTg                                             >  1:173712/1‑150 (MQ=255)
                                                                                                         acaGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGa                                            >  2:3469885/1‑150 (MQ=255)
                                                                                                           aGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAgg                                          <  1:8695615/150‑1 (MQ=255)
                                                                                                            gAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGgc                                         <  2:6435289/150‑1 (MQ=255)
                                                                                                                gCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGgccgcc                                     >  1:5183556/1‑150 (MQ=255)
                                                                                                                 cTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCa                                    >  2:1793946/1‑150 (MQ=255)
                                                                                                                 cTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTGGAGGCCGCCa                                    >  2:7058781/1‑150 (MQ=255)
                                                                                                                 cTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGAGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCa                                    >  1:8291826/1‑150 (MQ=255)
                                                                                                                  tGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCAt                                   >  1:4355026/1‑150 (MQ=255)
                                                                                                                   gCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCAtt                                  <  2:5121410/150‑1 (MQ=255)
                                                                                                                      tAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGAt                               >  1:6847838/1‑150 (MQ=255)
                                                                                                                      tAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGAt                               <  2:3866644/150‑1 (MQ=255)
                                                                                                                        aGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGc                             <  1:7112531/150‑1 (MQ=255)
                                                                                                                           ccGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCTGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGa                          >  2:611259/1‑150 (MQ=255)
                                                                                                                           ccGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGa                          >  2:6659529/1‑150 (MQ=255)
                                                                                                                            cGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGaa                         >  1:8413658/1‑150 (MQ=255)
                                                                                                                             gggAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCg                                                                                              <  1:614581/80‑1 (MQ=255)
                                                                                                                                aaTAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCGTGAGGCCGCCATTGATGCCGAAGGcc                     >  2:4499993/1‑150 (MQ=255)
                                                                                                                                  aaaacaaaGGGCGACTTCGGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCtt                   >  1:2794469/2‑150 (MQ=255)
                                                                                                                                   aaacaaaGGGCGACTTCGGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTc                  >  1:512756/1‑150 (MQ=255)
                                                                                                                                   aaaaaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTc                  <  1:1658891/150‑1 (MQ=255)
                                                                                                                                    aacaaaGGGCGACTTCGGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCa                 >  2:3539455/4‑150 (MQ=255)
                                                                                                                                     aaaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCGTGAGGCCGCCATTGATGCCGAAGGCCTTCAt                >  1:5501219/1‑150 (MQ=255)
                                                                                                                                      caaaGGGCGACTTCGGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATg               >  2:4940042/2‑150 (MQ=255)
                                                                                                                                      aaaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATg               >  2:8354451/1‑150 (MQ=255)
                                                                                                                                       aaaGGGCGACTTCGGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATgg              >  1:1483418/1‑150 (MQ=255)
                                                                                                                                       aaaGGGCGACTTCGGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATgg              >  2:8207864/1‑150 (MQ=255)
                                                                                                                                        aaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGgt             <  1:1459731/150‑1 (MQ=255)
                                                                                                                                        aaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGgt             >  2:4147336/1‑150 (MQ=255)
                                                                                                                                           ggCGACCGAAGTCGCCCTATGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGcc                                                    <  2:8647415/108‑1 (MQ=255)
                                                                                                                                            gCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTcc         >  2:7712046/1‑150 (MQ=255)
                                                                                                                                              gACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCgg       >  2:1600151/1‑150 (MQ=255)
                                                                                                                                                ccGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGCa     <  2:8451532/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTGATTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAAAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGCATGT  >  NC_002947/2098879‑2099175

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: