Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,018,866 C→T 16.3% intergenic (+66/+44) xpt → / ← PP_5266 xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,018,8660CT16.3% 366.5 / 57.3 155intergenic (+66/+44)xpt/PP_5266xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base C (48/81);  new base T (12/13);  total (60/95)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.69e-01

Only 100 of 157 total aligned reads displayed.

GAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCC  >  NC_002947/6018717‑6019014
                                                                                                                                                     |                                                                                                                                                     
gAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGcc                                                                                                                                                       <  1:5997118/150‑1 (MQ=255)
   cTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTa                                                                                                                                                    <  1:109518/150‑1 (MQ=255)
   cTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTa                                                                                                                                                    <  1:5842032/150‑1 (MQ=255)
   cTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTa                                                                                                                                                    <  1:468647/150‑1 (MQ=255)
   cTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTa                                                                                                                                                    <  2:6295085/150‑1 (MQ=255)
   cTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTa                                                                                                                                                    <  2:4554125/150‑1 (MQ=255)
    tGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTAc                                                                                                                                                   <  1:942294/150‑1 (MQ=255)
    tGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTAc                                                                                                                                                   >  1:5025616/1‑150 (MQ=255)
     ggACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCCGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTAcc                                                                                                                                                  <  2:242213/150‑1 (MQ=255)
      gACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCt                                                                                                                                                 <  1:5347909/150‑1 (MQ=255)
       aCAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtg                                                                                                                                                <  2:146912/150‑1 (MQ=255)
        cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgt                                                                                                                                               <  2:495369/150‑1 (MQ=255)
        cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgt                                                                                                                                               <  2:249494/150‑1 (MQ=255)
        cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgt                                                                                                                                               <  1:741871/150‑1 (MQ=25)
        cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgt                                                                                                                                               <  2:2194218/150‑1 (MQ=25)
        cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgt                                                                                                                                               <  1:1822967/150‑1 (MQ=25)
        cAGCCAGGGCTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgt                                                                                                                                               <  1:1822748/150‑1 (MQ=21)
         aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgtg                                                                                                                                              <  1:4257707/150‑1 (MQ=255)
         aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgtg                                                                                                                                              <  2:2412449/150‑1 (MQ=255)
         aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgtg                                                                                                                                              <  1:5429485/150‑1 (MQ=25)
          gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                                             <  2:3928538/150‑1 (MQ=255)
          gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                                             <  2:5488737/150‑1 (MQ=255)
          gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                                             <  1:3234769/150‑1 (MQ=255)
          gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTgg                                                                                                                                             <  2:4559312/150‑1 (MQ=25)
          gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTgg                                                                                                                                             <  1:4696991/150‑1 (MQ=25)
           ccAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACACGTGCAGCACAGGTCTCGAAAGCTGTGCCGGACCGGTggg                                                                                                                                            >  2:5941506/1‑150 (MQ=255)
           ccAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCCCGCCCACAGGTGCAGCACAGGTATCGAACGCGGTGCCCTACCAGTggg                                                                                                                                            >  2:5941695/1‑150 (MQ=255)
            cAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACCCGGGCAGCACAAGGCGCGGACCCGGGGGCGGACCAGTGGGa                                                                                                                                           >  1:6193276/1‑150 (MQ=255)
                gTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCgg                                                                                                                                      <  1:275773/150‑1 (MQ=255)
                gTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCgg                                                                                                                                      <  1:3556380/150‑1 (MQ=255)
                 ttACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCggg                                                                                                                                     <  1:988109/150‑1 (MQ=255)
                           gAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCcgcg                                                                                                                           >  2:5003306/1‑150 (MQ=255)
                             aTCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGCACGCTGTGCCGTACCAGTGGGAGCGGGCGTGCCCGCGaa                                                                                                                         >  2:5533375/1‑150 (MQ=255)
                                  tGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTACAGCGAATAAACGCGACCCCTGCTTCGCGGGCACGCCCGCTCCCCCCGGTGCAGCACTGGTCTCGAAGGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGcc                                                                                                                    <  2:4847596/150‑1 (MQ=37)
                                          gTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCt                                                                                                            >  2:3946405/1‑150 (MQ=255)
                                            gAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAACCGCGACACCGGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCttt                                                                                                          <  1:5122153/150‑1 (MQ=255)
                                                tCGCTGGACGTTGGGGTGGTCAGCTTCATCGAATAAACGCCACACCGGCTTCGCGGGCACGCCCGCTCCCACTGGTGCAGCACCGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATc                                                                                                      <  2:2516336/150‑1 (MQ=37)
                                                tCGCTGGACGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATc                                                                                                      <  2:2053485/150‑1 (MQ=255)
                                                tCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACCCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATc                                                                                                      >  2:6688867/1‑150 (MQ=255)
                                                 cGCTGGTTGGTGGTGTGGTCAGCTTCATCGATTAAACGCTACACCTGCTTCGGGGGCACGCCCGCTCCCACTGGTGCAGCACTGGTCTCGAACGCTGTGCCGTACCTGGGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCt                                                                                                     <  1:5050679/150‑1 (MQ=38)
                                                  gCTGGATGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTg                                                                                                    <  2:5964877/150‑1 (MQ=255)
                                                  gCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGGACCAGTGGGAGCGGGCGGGCCCCCGAATGGCCGCAACCCGTTTATCTg                                                                                                    >  2:1667310/1‑150 (MQ=38)
                                                  gCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGGCTCGAACCCAGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGCCGCAACCCGTTTATATg                                                                                                    >  2:457186/1‑150 (MQ=37)
                                                  gCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACACGTGCAGCAAAGGGGTCGAACGCTGTGCCGGAACGGTGGGAGCGGGCGTGCCCCCGAAAGGCCGCAAACCATTTAACTg                                                                                                    >  2:3466607/1‑150 (MQ=37)
                                                   gtGGATGGTGGTGGGGTCAGCTACATCGAATAATCGCGACACCGGCTTCGCGGGCACGCCCGCTCCCCCTGGTGCAGCACCGGTCCCTACCGGTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTgg                                                                                                   <  2:3069350/149‑1 (MQ=37)
                                                           tggtgtggtCAGCTTCATTGAATAAACGCGACACCTGCTTCGGGGGCACGCCCGCTCCCACCGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTgccgcc                                                                                           <  2:2374967/150‑1 (MQ=255)
                                                                        ttCATCGAATAACCGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACTGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCcaaca                                                                              <  1:6685538/150‑1 (MQ=255)
                                                                        ttCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACCGGTGCAGCACAGGTCACGAACGCTGTGCCGTCCCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCcaaca                                                                              <  1:4546182/150‑1 (MQ=255)
                                                                        ttCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCCCTGCCACCGGTTCACCACACCGCTAGAACGCGGTGCCCGACCTGTGGGAGCGGGCCAGGCACCCAAAGGCCCCCACCCATTTATCTTGTTTCCGCCACCCCCCCccacc                                                                              >  2:2552637/1‑149 (MQ=255)
                                                                                  caaCGCGACACCTGCTTCGCGGGCCCGCCCGCGCCCACCGGTGCAGCACAGGTCTCGAACGCTGTGCCGTCCCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTg                                                                    <  1:3864049/149‑1 (MQ=255)
                                                                                  aatcgcgACACGTGGTTCGCGGGCCGGCCCGCTCCCACCGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCGGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTg                                                                    <  1:3801593/147‑1 (MQ=38)
                                                                                  aaccgcgACACCTGCTTGGCGGGCACGCCCGCTCCCACCGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTg                                                                    <  1:3320212/147‑1 (MQ=255)
                                                                                  aaaCGCGACACCTGCTTCGGGGGCACGCCCGCTCCCACTGGTGCAGCACATGTCTCTAAGGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTg                                                                    <  1:2068920/150‑1 (MQ=38)
                                                                                  aaaCGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACGTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTg                                                                    <  2:1702501/150‑1 (MQ=255)
                                                                                  aaaCGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTg                                                                    <  2:4155588/150‑1 (MQ=255)
                                                                                  aaaCGCGACACCTGCTGCGCGGGCACGCCCGCTCCCACAGGTGCAGCACTGTTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTg                                                                    <  2:746209/150‑1 (MQ=255)
                                                                                       cgACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGATCGCTGTGCCGTAGCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTaaa                                                               <  2:4810089/150‑1 (MQ=255)
                                                                                       cgACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTaaa                                                               <  2:5609211/150‑1 (MQ=255)
                                                                                              tGCTTCGCGGGCACGCCCGCTCCCCCACGTGCAGCACAGGTCTCCCACGCGGTTCCGGACCTGTGGGAGCGGGCGGGCCCCCCAATGGCCGCAACCCGTTTATCTGGTTTCCCCCCCCCCCCCCACCAACCAGCGCGGGGAAAGttaatg                                                        >  2:786506/1‑146 (MQ=255)
                                                                                              tGCTTCGCGGGCACCCCCGCTCCCCCAGGGGCAGCACAGGTCTAGAACGCAGTGCCCTACCAGTGGGACCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCAGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAActtcttg                                                        >  2:4761551/1‑149 (MQ=255)
                                                                                                 ttCGCGGGCACGCCCGCTCCCACCGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCttt                                                     <  2:1528364/150‑1 (MQ=255)
                                                                                                        gCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGGGCTGCCCCGGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCGTGTATCTGGTTGCCGCCAGCCCCGCCCACCACAAGCGCTGGTAAAGTTCGTCTGTAAGcccc                                              >  2:1056387/1‑150 (MQ=17)
                                                                                                              ccGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGc                                        >  1:2642384/1‑150 (MQ=255)
                                                                                                              ccGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGt                                                                  <  1:4953241/124‑1 (MQ=17)
                                                                                                              ccGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGt                                                                  <  1:4953517/124‑1 (MQ=17)
                                                                                                              ccGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGt                                                                  >  2:4953241/1‑124 (MQ=17)
                                                                                                               cGCTCCCACCGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGGGGGAGCGGGCGTGCGCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCt                                       <  2:5025616/150‑1 (MQ=255)
                                                                                                                 cTCCCACGGGTGCAGCACCGGTCTCGAACGCTGTGCCGTACCTGTGGGCGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGc                                     <  1:3108033/150‑1 (MQ=255)
                                                                                                                 cTCCCACCGGTGCAGCACAGGTCTCGAACGCTGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGc                                     <  1:1638564/150‑1 (MQ=255)
                                                                                                                 cTCCCACAGGTGCAGCACAGGTCTCGACCGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCGGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGc                                     <  2:4446880/150‑1 (MQ=255)
                                                                                                                 cTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGc                                     <  1:4467077/150‑1 (MQ=255)
                                                                                                                  tCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGGCGTGCCCGCGAATGGCCGCAACCCTGTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGc                                     >  2:6215788/1‑150 (MQ=255)
                                                                                                                  tCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcc                                    >  2:504869/1‑150 (MQ=21)
                                                                                                                  tCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcc                                    >  2:504970/1‑150 (MQ=21)
                                                                                                                   cccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                                   <  2:5550303/150‑1 (MQ=255)
                                                                                                                    ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc                                  >  1:855199/1‑150 (MQ=255)
                                                                                                                     cacaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCa                                 >  2:2595625/1‑150 (MQ=255)
                                                                                                                      acaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa                                >  2:5918183/1‑150 (MQ=255)
                                                                                                                      acaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa                                >  1:6728261/1‑150 (MQ=255)
                                                                                                                      acaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa                                >  2:3871459/1‑150 (MQ=255)
                                                                                                                      acaGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa                                >  1:2038572/1‑150 (MQ=25)
                                                                                                                      acaGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa                                >  1:1612793/1‑150 (MQ=25)
                                                                                                                        aGGTGCAGCACCGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACt                              <  2:5788186/150‑1 (MQ=255)
                                                                                                                        aGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACt                              >  2:4004003/1‑150 (MQ=255)
                                                                                                                          gTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGc                            >  2:6706682/1‑150 (MQ=255)
                                                                                                                                ctgaGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGc                      <  1:4399097/147‑1 (MQ=255)
                                                                                                                                cacaGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTAACTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGc                      >  2:5920148/1‑150 (MQ=255)
                                                                                                                                   aGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGcccc                   >  2:3193240/1‑150 (MQ=255)
                                                                                                                                     gTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAg                 <  2:6416982/150‑1 (MQ=255)
                                                                                                                                     gTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAg                 >  2:6547703/1‑150 (MQ=255)
                                                                                                                                      tctcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAgg                >  1:1489505/1‑150 (MQ=255)
                                                                                                                                        tcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGttc                                                           >  2:4339763/1‑105 (MQ=255)
                                                                                                                                        tcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGttc                                                           <  1:4339763/105‑1 (MQ=255)
                                                                                                                                        tcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTg              >  1:3314246/1‑150 (MQ=255)
                                                                                                                                         cgagaccTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTgg             >  2:979537/7‑150 (MQ=255)
                                                                                                                                         cGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTgg             <  2:408348/150‑1 (MQ=255)
                                                                                                                                         cGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTgg             >  1:6030699/1‑150 (MQ=255)
                                                                                                                                          gAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAgc                                                                        >  1:2785912/1‑90 (MQ=255)
                                                                                                                                          gAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGc            >  2:2178477/1‑150 (MQ=255)
                                                                                                                                                    ccGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATAccc  >  1:5510963/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCC  >  NC_002947/6018717‑6019014

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