Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,270,479 C→A 21.8% G365G (GGC→GGA PP_4643 → xanthine/uracil permease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,270,4790CA21.8% 387.6 / 101.7 174G365G (GGC→GGAPP_4643xanthine/uracil permease family protein
Reads supporting (aligned to +/- strand):  ref base C (53/83);  new base A (21/17);  total (74/100)
Fisher's exact test for biased strand distribution p-value = 9.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 175 total aligned reads displayed.

CTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGA  >  NC_002947/5270330‑5270627
                                                                                                                                                     |                                                                                                                                                    
     ggCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGt                                                                                                                                                 >  1:3799716/1‑150 (MQ=255)
       ccGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTgg                                                                                                                                               <  1:4984703/150‑1 (MQ=255)
        cGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTggc                                                                                                                                              <  1:1915071/150‑1 (MQ=255)
             ccGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGAACCGTGGCGGCCa                                                                                                                                         >  1:3713892/1‑150 (MQ=255)
              cGGCGTGAAAAGCCGTTATGTGGTGGCAACGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAg                                                                                                                                        <  2:2821603/150‑1 (MQ=255)
              cGGCGTGAAAAGCCGTTATGTGGTGGCAACGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAg                                                                                                                                        >  1:4664708/1‑150 (MQ=255)
                gggTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGTCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCg                                                                                                                                      <  1:4224928/150‑1 (MQ=255)
                gggTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCg                                                                                                                                      <  1:2104048/150‑1 (MQ=255)
                gggTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCg                                                                                                                                      <  1:1252754/150‑1 (MQ=255)
                     aaaagCCGTTATGTGGTGGCAACGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATc                                                                                                                                 >  2:4054926/4‑150 (MQ=255)
                     aaGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATc                                                                                                                                 >  1:2029429/1‑150 (MQ=255)
                      agagCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATcc                                                                                                                                <  2:1561163/150‑1 (MQ=255)
                      agagCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATcc                                                                                                                                <  1:1727396/150‑1 (MQ=255)
                      agagCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATcc                                                                                                                                <  1:1518954/150‑1 (MQ=255)
                       gagCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCg                                                                                                                               <  1:3428298/150‑1 (MQ=255)
                           cGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCAcc                                                                                                                           <  1:2112837/150‑1 (MQ=255)
                               atgtggtgGCAACGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACttt                                                                                                                       <  2:2136742/148‑3 (MQ=255)
                                 gtggtgGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTcc                                                                                                                     <  1:1498972/150‑1 (MQ=255)
                                   ggtggCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCaa                                                                                                                   <  2:2075984/150‑1 (MQ=255)
                                   ggtggCAACGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCaaa                                                                                                                   <  1:768400/150‑1 (MQ=255)
                                   ggtggCAACGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCaaa                                                                                                                   <  1:2939756/150‑1 (MQ=255)
                                   ggtggCAACGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCaaa                                                                                                                   <  2:4840779/150‑1 (MQ=255)
                                       gcaaCGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTc                                                                                                               <  2:5067644/147‑1 (MQ=255)
                                        cTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCa                                                                                                              <  2:1746752/150‑1 (MQ=255)
                                         tACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGTGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGACAg                                                                                                             >  2:1291679/1‑150 (MQ=255)
                                         aaCGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTcga                                                                                                             >  2:3738023/2‑148 (MQ=255)
                                              gcggcgTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACa                                                                                                        >  2:828671/1‑150 (MQ=255)
                                              gcggcgTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACa                                                                                                        <  1:4314563/150‑1 (MQ=255)
                                              gcggcgTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACa                                                                                                        <  1:3842165/150‑1 (MQ=255)
                                               cggcgTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACaa                                                                                                       <  2:1668913/150‑1 (MQ=255)
                                                  cgTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGat                                                                                                    <  2:4819094/150‑2 (MQ=255)
                                                   gTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGaac                                                                                                   <  1:4546045/150‑1 (MQ=255)
                                                   gTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGaac                                                                                                   >  2:5072919/1‑150 (MQ=255)
                                                      aTTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGaacaac                                                                                                <  2:2195285/150‑1 (MQ=255)
                                                          tggtggTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCa                                                                                            >  1:3957203/1‑150 (MQ=255)
                                                          tggtggTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCa                                                                                            >  2:4272855/1‑150 (MQ=255)
                                                            gtggtACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAAc                                                                                          <  2:4380798/150‑1 (MQ=255)
                                                             tcacaCTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCaatc                                                                                         >  1:4037168/5‑148 (MQ=255)
                                                              cacaCTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCt                                                                                        <  1:2603116/147‑1 (MQ=255)
                                                                tACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGa                                                                                      <  1:2894281/150‑1 (MQ=255)
                                                                 aCTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGat                                                                                     >  1:810029/1‑150 (MQ=255)
                                                                 aCTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGat                                                                                     >  2:3187726/1‑150 (MQ=255)
                                                                 aCTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCTTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGat                                                                                     >  2:1159488/1‑150 (MQ=255)
                                                                 aCTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGat                                                                                     >  1:1668128/1‑150 (MQ=255)
                                                                 aCTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGat                                                                                     >  2:4526020/1‑150 (MQ=255)
                                                                  cTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGatc                                                                                    >  1:1667635/1‑150 (MQ=255)
                                                                  cTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGatc                                                                                    >  1:4558170/1‑150 (MQ=255)
                                                                   tCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGatca                                                                                   >  2:4947319/1‑150 (MQ=255)
                                                                     ggCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGatcatc                                                                                 >  1:4856913/1‑150 (MQ=255)
                                                                     ggCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGatcatc                                                                                 <  1:1358443/150‑1 (MQ=255)
                                                                      gCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCAtcg                                                                                >  2:4661375/1‑150 (MQ=255)
                                                                        ttGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCAtcgtc                                                                              <  2:2532517/150‑1 (MQ=255)
                                                                         ttctgcCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCAttgtgg                                                                             >  1:1344580/3‑148 (MQ=255)
                                                                          tctgcCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTcgc                                                                            <  2:2875540/149‑1 (MQ=255)
                                                                            tgcCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCgg                                                                          >  1:2460541/1‑150 (MQ=255)
                                                                              ccAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGcc                                                                        <  2:729641/150‑1 (MQ=255)
                                                                              ccAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGcc                                                                        >  1:1997780/1‑150 (MQ=255)
                                                                                aGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTc                                                                      <  1:2203702/150‑1 (MQ=255)
                                                                                  tCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGc                                                                    >  2:43450/1‑150 (MQ=255)
                                                                                  tCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGc                                                                    <  1:3172046/150‑1 (MQ=255)
                                                                                   cATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCt                                                                   <  2:1757485/150‑1 (MQ=255)
                                                                                     tGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTgg                                                                 <  2:1724488/150‑1 (MQ=255)
                                                                                     tGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTgg                                                                 <  2:199322/150‑1 (MQ=255)
                                                                                     tGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTgg                                                                 >  2:3818725/1‑150 (MQ=255)
                                                                                      gggACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCg                                                                                                                               <  1:947836/87‑1 (MQ=255)
                                                                                            ggTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCgg                                                          <  1:5093440/150‑1 (MQ=255)
                                                                                               gATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCAt                                                       <  2:929494/150‑1 (MQ=255)
                                                                                                aTTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATg                                                      <  1:3703536/150‑1 (MQ=255)
                                                                                                    ccgccgTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATcc                                                  >  1:1339580/1‑150 (MQ=255)
                                                                                                        cgTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGAt                                              <  1:1169313/150‑1 (MQ=255)
                                                                                                            ccTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGcc                                          >  1:3110636/1‑150 (MQ=255)
                                                                                                             cTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCg                                         <  1:3284562/150‑1 (MQ=255)
                                                                                                                ccccGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCAc                                      <  2:1496456/150‑1 (MQ=255)
                                                                                                                 agcagtgcTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCAcc                                     >  2:2276473/8‑150 (MQ=255)
                                                                                                                 agcagtgcTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCAcc                                     >  2:2276481/8‑150 (MQ=255)
                                                                                                                  ccGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCg                                    >  1:1181783/1‑150 (MQ=255)
                                                                                                                     gtgcTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCAccagaa                                 >  2:716111/4‑146 (MQ=255)
                                                                                                                     gTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGAcc                                 <  2:698733/150‑1 (MQ=255)
                                                                                                                       gcTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCAccagaatt                               <  2:1070928/149‑7 (MQ=255)
                                                                                                                         tGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCt                             <  1:3194464/150‑1 (MQ=255)
                                                                                                                         tGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCt                             <  1:3194471/150‑1 (MQ=255)
                                                                                                                          gggCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTa                            <  1:614301/150‑1 (MQ=255)
                                                                                                                           ggcggcGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTac                           <  2:292890/150‑1 (MQ=255)
                                                                                                                            gcggcgCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTacc                          <  2:3575522/150‑1 (MQ=255)
                                                                                                                              ggcgCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTaccac                        <  1:2660206/150‑1 (MQ=255)
                                                                                                                              ggcgCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTaccac                        <  1:4833517/150‑1 (MQ=255)
                                                                                                                                tgcCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTAccatca                      <  2:848947/149‑1 (MQ=255)
                                                                                                                                   cGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGtt                   >  1:4793105/1‑150 (MQ=255)
                                                                                                                                    ggTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTc                  <  2:1806056/150‑1 (MQ=255)
                                                                                                                                     gTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTcc                 <  2:2663985/150‑1 (MQ=255)
                                                                                                                                         cGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGaa             <  1:2624323/150‑1 (MQ=255)
                                                                                                                                           tGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCggcct           >  1:1697375/1‑146 (MQ=255)
                                                                                                                                            gCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTg          <  1:2332768/150‑1 (MQ=255)
                                                                                                                                            gCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCggcctg          <  2:3131713/150‑6 (MQ=255)
                                                                                                                                              tGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGt        >  2:2976275/1‑150 (MQ=255)
                                                                                                                                              tGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGt        <  2:2064861/150‑1 (MQ=255)
                                                                                                                                               gTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCACCAGAATTCTACCATCAGTTCCCGGCCTGGtt       >  1:664868/1‑150 (MQ=255)
                                                                                                                                                ttCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTc      <  2:3167839/150‑1 (MQ=255)
                                                                                                                                                    gCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCgaga  <  2:4872465/150‑1 (MQ=255)
                                                                                                                                                    gCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAAAGTCAAACTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCgaga  <  1:2744090/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGA  >  NC_002947/5270330‑5270627

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: