Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,739,085 G→C 17.0% A113A (GCC→GCG hutU ← urocanate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,739,0850GC17.0% 509.3 / 90.3 194A113A (GCC→GCGhutUurocanate hydratase
Reads supporting (aligned to +/- strand):  ref base G (78/83);  new base C (18/15);  total (96/98)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 195 total aligned reads displayed.

CGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGC  >  NC_002947/5738937‑5739233
                                                                                                                                                    |                                                                                                                                                    
    gTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACCTCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCa                                                                                                                                                 <  2:3014643/150‑1 (MQ=255)
    gTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCa                                                                                                                                                 <  1:205447/150‑1 (MQ=255)
     taatGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCaa                                                                                                                                                >  2:4241487/4‑149 (MQ=255)
        tGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAgtgt                                                                                                                                             <  2:309770/150‑1 (MQ=255)
        tGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAatgc                                                                                                                                             >  2:3504344/1‑149 (MQ=255)
            ggCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCa                                                                                                                                         <  1:1620008/150‑1 (MQ=255)
            ggCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCa                                                                                                                                         >  1:979187/1‑150 (MQ=255)
              cGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGCAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTGCCCAGTTGGCCCAGTGTGGCAcc                                                                                                                                       <  1:137770/150‑1 (MQ=255)
               gACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCa                                                                                                                                      >  1:1058861/1‑150 (MQ=255)
                   ggCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGtt                                                                                                                                  <  2:4098003/150‑1 (MQ=255)
                    gCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCATCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTg                                                                                                                                 <  1:4364817/150‑1 (MQ=255)
                    gCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTg                                                                                                                                 >  2:1270436/1‑150 (MQ=255)
                      ttCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGa                                                                                                                               <  1:3932345/150‑1 (MQ=255)
                      ttCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGa                                                                                                                               <  1:171685/150‑1 (MQ=255)
                         gacgaAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGtt                                                                                                                            <  1:4210459/150‑1 (MQ=255)
                         gacgaAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGtt                                                                                                                            >  2:1896528/1‑150 (MQ=255)
                         gacgaAGGGTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGtt                                                                                                                            <  2:2121392/150‑1 (MQ=255)
                         cacgaAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAAtt                                                                                                                            <  1:3634626/149‑1 (MQ=255)
                          acgaAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTg                                                                                                                           <  2:5323374/150‑1 (MQ=255)
                           cgaAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTgg                                                                                                                          <  1:3408698/150‑1 (MQ=255)
                           cgaAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTgg                                                                                                                          <  1:1170298/150‑1 (MQ=255)
                           cgaAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGACATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTgg                                                                                                                          >  1:648506/1‑150 (MQ=255)
                            gaAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGc                                                                                                                         >  1:2846065/1‑150 (MQ=255)
                            gaAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGc                                                                                                                         >  1:3492593/1‑150 (MQ=255)
                            gaAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGc                                                                                                                         >  1:2811767/1‑150 (MQ=255)
                            gaAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAAtt                                                                                                                            <  2:2620930/147‑1 (MQ=255)
                                gTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATc                                                                                                                     <  1:1270436/150‑1 (MQ=255)
                                 tttCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCa                                                                                                                    <  2:4573303/150‑1 (MQ=255)
                                       aGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCAcgc                                                                                                              <  1:2657113/150‑1 (MQ=255)
                                          tGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCAcgcgcg                                                                                                           >  2:4860935/1‑150 (MQ=255)
                                          tGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTGGCGATGAGAAcgcgcg                                                                                                           <  2:4173301/150‑1 (MQ=255)
                                           gCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCa                                                                                                                 >  2:4032880/1‑143 (MQ=255)
                                           gCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCa                                                                                                                 >  2:2253135/1‑143 (MQ=255)
                                           gCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCa                                                                                                                 <  1:4032880/143‑1 (MQ=255)
                                              cTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGc                                                                                                       <  2:2502212/150‑1 (MQ=255)
                                              cTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGcc                                                                                                                                                   >  2:702416/1‑106 (MQ=255)
                                                   cGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGc                                                                                                  <  1:639685/150‑1 (MQ=255)
                                                    gATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCt                                                                                                 >  1:809932/1‑150 (MQ=255)
                                                       gCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCtgtg                                                                                              <  1:4078441/150‑1 (MQ=255)
                                                        cccTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTgg                                                                                             <  1:5171402/150‑1 (MQ=255)
                                                           tgtgcaccGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTGGCGATGAGAACGCGCGGTGCATCCTTGTGGGtg                                                                                          >  1:732237/7‑149 (MQ=255)
                                                                gCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGaa                                                                                     <  1:895972/150‑1 (MQ=255)
                                                                accGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTGGCGATGAGAACGCGCGGTGCATCCTTGTGGGtgcggaa                                                                                     >  2:2343858/2‑144 (MQ=255)
                                                                accGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTGGCGATGAGAACGCGCGGTGCATCCTTGTGGGtgcggaa                                                                                     >  2:5037294/2‑144 (MQ=255)
                                                                   gATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAacac                                                                                  >  2:2445859/1‑150 (MQ=255)
                                                                       tAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACaccg                                                                              >  2:46970/1‑150 (MQ=255)
                                                                         gATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACaccgac                                                                            <  1:3865902/150‑1 (MQ=255)
                                                                         gATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACaccgac                                                                            >  2:299327/1‑150 (MQ=255)
                                                                         gATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACaccgac                                                                            >  1:4911098/1‑150 (MQ=255)
                                                                          aTCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACaccgacc                                                                           >  1:3130281/1‑150 (MQ=255)
                                                                           tCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCGTGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACaccgaccg                                                                          <  2:757615/150‑1 (MQ=255)
                                                                              aGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCt                                                                       >  2:3780449/1‑150 (MQ=255)
                                                                              aGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCt                                                                       >  1:3747166/1‑150 (MQ=255)
                                                                              aGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCt                                                                       >  1:1731115/1‑150 (MQ=255)
                                                                               gCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCtt                                                                      <  2:897714/150‑1 (MQ=255)
                                                                                     ggcggTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGGTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAAGCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGa                                                                >  2:3795185/1‑150 (MQ=255)
                                                                                          tCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCa                                                           <  1:4860935/150‑1 (MQ=255)
                                                                                          gggccTGGCGATGTATGGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCa                                                           >  2:1573876/5‑150 (MQ=255)
                                                                                           ggccTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCAc                                                          >  1:3831481/4‑150 (MQ=255)
                                                                                             tCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTggg                                                                                            <  1:2354406/114‑1 (MQ=255)
                                                                                               tGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTGGCGATGAGAACGCGCGGTGCATCCTTGTGGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACcagc                                                      >  2:3478572/1‑150 (MQ=255)
                                                                                                  ccATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTGGCGATGAGAACGCGCGGTGCATCCTTGTGGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACcagcaag                                                   >  2:1724442/1‑148 (MQ=255)
                                                                                                  ccATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAgg                                                   <  2:4188475/150‑1 (MQ=255)
                                                                                                  ccATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAgg                                                   <  2:4303909/150‑1 (MQ=255)
                                                                                                   cATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGAGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAggg                                                  >  2:984874/1‑150 (MQ=255)
                                                                                                    aTACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTGGCGATGAGAACGCGCGGTGCATCCTTGTGGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAGCAAGGt                                                 >  1:2037436/1‑150 (MQ=255)
                                                                                                     tACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTc                                                                                          <  2:3861139/108‑1 (MQ=255)
                                                                                                        aTCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcg                                             <  2:3843481/150‑1 (MQ=255)
                                                                                                        aTCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcg                                             >  2:2368139/1‑150 (MQ=255)
                                                                                                         tCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGTTCACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcgt                                            >  2:2060353/1‑150 (MQ=255)
                                                                                                         tCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAgg                                                   <  2:3087799/143‑1 (MQ=255)
                                                                                                         tCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAgg                                                   >  1:3087799/1‑143 (MQ=255)
                                                                                                         tCGCCAGGCCCTTGGCGACCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcga                                            >  2:3309363/1‑149 (MQ=255)
                                                                                                          cGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcgtc                                           <  1:4275491/150‑1 (MQ=255)
                                                                                                          cGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcgtc                                           <  2:2790727/150‑1 (MQ=255)
                                                                                                          cGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcgtc                                           <  1:5031752/150‑1 (MQ=255)
                                                                                                            ccAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcg                                             >  2:2963993/1‑146 (MQ=255)
                                                                                                            ccAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcg                                             <  1:2963993/146‑1 (MQ=255)
                                                                                                             cAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCgtcgtc                                        <  2:2890039/150‑1 (MQ=255)
                                                                                                              aGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCt                                       >  1:5286535/1‑150 (MQ=255)
                                                                                                               ggCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGGTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCtt                                      >  2:2681609/1‑150 (MQ=255)
                                                                                                                  ccTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCa                                   <  2:5295076/150‑1 (MQ=255)
                                                                                                                       gCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCgg                              <  2:2231196/150‑1 (MQ=255)
                                                                                                                         gTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGt                            <  1:140532/150‑1 (MQ=255)
                                                                                                                           ccAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAacaca                                                                                 <  1:1899852/95‑1 (MQ=255)
                                                                                                                           ccAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAacaca                                                                                 >  2:1899852/1‑95 (MQ=255)
                                                                                                                            cAGTTCGTTGACGTGGTCCCAGTTGGCCCAGTGTGGCACCAGGCTGGAGTTGGCAATCAGCACGCGCGGTGCGCTGCTGTGGGTCTTGAACCCACCGACCGGCTTGCCCGACGGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAg                         <  2:5394533/150‑1 (MQ=255)
                                                                                                                             aGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAgg                        <  1:1573876/150‑1 (MQ=255)
                                                                                                                              gTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAggg                       >  2:1145412/1‑150 (MQ=255)
                                                                                                                               ttCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGt                      >  1:5359378/1‑150 (MQ=255)
                                                                                                                                    tGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTcga                 >  1:4958575/1‑150 (MQ=255)
                                                                                                                                     gAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTcgac                <  1:1010360/150‑1 (MQ=255)
                                                                                                                                      aaGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTcgacg               >  2:4489980/1‑150 (MQ=255)
                                                                                                                                       aGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAgg                                                   <  1:2576207/113‑1 (MQ=255)
                                                                                                                                            tCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTg         <  1:3477750/150‑1 (MQ=255)
                                                                                                                                              ccAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTc       >  1:4469106/1‑150 (MQ=255)
                                                                                                                                               cAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCg      >  1:2291364/1‑150 (MQ=255)
                                                                                                                                                 gTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTa    >  1:2454671/1‑150 (MQ=255)
                                                                                                                                                 gTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTa    >  1:927735/1‑150 (MQ=255)
                                                                                                                                                   tGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGc  <  1:3016421/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGC  >  NC_002947/5738937‑5739233

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: