Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,819,917 C→T 8.7% intergenic (‑35/+295) spy ← / ← astE periplasmic ATP‑independent protein refolding chaperone, stress‑induced/succinylglutamate desuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,819,9170CT 8.7% 921.6 / 51.2 274intergenic (‑35/+295)spy/astEperiplasmic ATP‑independent protein refolding chaperone, stress‑induced/succinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base C (118/132);  new base T (16/8);  total (134/140)
Fisher's exact test for biased strand distribution p-value = 8.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 275 total aligned reads displayed.

CATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGA  >  CP009273/1819769‑1820065
                                                                                                                                                    |                                                                                                                                                    
 atcatcGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGttt                                                                                                                                                    >  2:4033037/1‑150 (MQ=255)
 atcatcGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCtt                                                                                                                                                    >  1:2088695/1‑150 (MQ=255)
    atcGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGTTTTCt                                                                                                                                                 >  2:1494592/1‑150 (MQ=255)
    atcGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCt                                                                                                                                                 <  1:3013881/150‑1 (MQ=255)
      cGGCTTCGCGTCAGCCTGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGAATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTa                                                                                                                                               <  2:4091655/150‑1 (MQ=255)
       ggCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGTTTTCTTaa                                                                                                                                              <  1:1672721/150‑1 (MQ=255)
        gCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAg                                                                                                                                             <  1:1043244/150‑1 (MQ=255)
         cTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGt                                                                                                                                            >  1:3726766/1‑150 (MQ=255)
          ttCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTa                                                                                                                                           >  2:1306714/1‑150 (MQ=255)
          ttCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTa                                                                                                                                           >  1:928045/1‑150 (MQ=255)
              cgTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGt                                                                                                                                       >  2:948893/1‑150 (MQ=255)
               gTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTg                                                                                                                                      <  2:2759393/150‑1 (MQ=255)
                 catccGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGcc                                                                                                                                    >  1:2405976/4‑150 (MQ=255)
                 cAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGcc                                                                                                                                    >  2:310664/1‑150 (MQ=255)
                 cAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGcc                                                                                                                                    >  1:3809817/1‑150 (MQ=255)
                 cAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGcc                                                                                                                                    >  1:2888435/1‑150 (MQ=255)
                 cAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGcc                                                                                                                                    >  1:1480393/1‑150 (MQ=255)
                    ccGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGTTTTCTTAAGTAGCGTGCCgtg                                                                                                                                 <  1:3387522/150‑1 (MQ=255)
                      ggTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGtt                                                                                                                               <  2:586099/150‑1 (MQ=255)
                      ggTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGtt                                                                                                                               <  2:4386573/150‑1 (MQ=255)
                       gTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTg                                                                                                                              <  2:1653604/150‑1 (MQ=255)
                       gTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTg                                                                                                                              <  2:3178936/150‑1 (MQ=255)
                        tgctgcGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGa                                                                                                                             >  2:2242227/1‑150 (MQ=255)
                         gctgcGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGAc                                                                                                                            <  1:1220518/150‑1 (MQ=255)
                          ctgcGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACg                                                                                                                           >  1:4010494/1‑150 (MQ=255)
                              ggTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGAt                                                                                                                       <  2:1511228/150‑1 (MQ=255)
                              ggTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGAt                                                                                                                       <  1:590650/150‑1 (MQ=255)
                               gTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGAtc                                                                                                                      <  1:1830335/150‑2 (MQ=255)
                                 aGTGGTGTCTGCGGCATGGGCCAGGTTAGCNGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTaa                                                                                                                    >  1:4283627/1‑150 (MQ=255)
                                  gtggtgTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAAc                                                                                                                   >  1:4057988/1‑150 (MQ=255)
                                    ggtgTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGTTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACtt                                                                                                                 >  1:2158684/1‑150 (MQ=255)
                                         cTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTg                                                                                                            <  2:1816145/150‑1 (MQ=255)
                                          tGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGa                                                                                                           <  2:3519158/150‑1 (MQ=255)
                                            cGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACt                                                                                                         <  2:2932364/150‑1 (MQ=255)
                                             ggCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACtt                                                                                                        >  1:80417/1‑150 (MQ=255)
                                             ggCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACtt                                                                                                        >  2:4596307/1‑150 (MQ=255)
                                               cATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTa                                                                                                      <  2:1683283/150‑1 (MQ=255)
                                               cATGGGCAAGGTTATCCGCGCCTAGAGCCAGGGTAGATGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTa                                                                                                      <  1:279284/150‑1 (MQ=255)
                                                  gggCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCg                                                                                                   >  2:1146377/1‑150 (MQ=255)
                                                       aGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAAt                                                                                              <  2:1008175/150‑1 (MQ=255)
                                                       aGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAAt                                                                                              <  2:842918/150‑1 (MQ=255)
                                                        ggTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAAtt                                                                                             <  1:1374133/150‑1 (MQ=255)
                                                           tAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAAt                                                                                          >  1:1481784/1‑150 (MQ=255)
                                                            aGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATc                                                                                         >  2:448596/1‑150 (MQ=255)
                                                             gCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCa                                                                                        <  1:2238528/150‑1 (MQ=255)
                                                                gcgcCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCagag                                                                                     <  1:2949968/150‑1 (MQ=255)
                                                                  gcCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCa                                                                                   >  1:4578273/1‑150 (MQ=255)
                                                                    cAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAAc                                                                                 >  2:1273373/1‑150 (MQ=255)
                                                                    cAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAAc                                                                                 >  1:1036788/1‑150 (MQ=255)
                                                                      agagCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACgg                                                                               <  2:1710284/150‑1 (MQ=255)
                                                                        agCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTa                                                                             >  1:4522945/1‑150 (MQ=255)
                                                                             gggTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACa                                                                        <  2:1950569/150‑1 (MQ=255)
                                                                                      caacaaACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAgtg                                                               <  1:2636147/150‑1 (MQ=255)
                                                                                       aacaaacaGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAgtgt                                                              <  2:1437847/150‑1 (MQ=255)
                                                                                        acaaacaGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGTTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATAAGAGCAACGGTAAAACAATGAAAGTGTa                                                             >  2:3665044/1‑150 (MQ=255)
                                                                                        acaaacaGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTa                                                             >  1:3262653/1‑150 (MQ=255)
                                                                                           aacaGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTaaaa                                                          <  1:4425875/150‑1 (MQ=255)
                                                                                            acaGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGTTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTaaaaa                                                         >  1:2702009/1‑150 (MQ=255)
                                                                                              aGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAac                                                       >  2:3486850/1‑150 (MQ=255)
                                                                                               gTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAaca                                                      >  1:518931/1‑150 (MQ=255)
                                                                                                  cAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACtt                                                   >  2:2708039/1‑150 (MQ=255)
                                                                                                   aGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACttt                                                  >  1:141826/1‑150 (MQ=255)
                                                                                                    gTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGTTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACtttt                                                 <  1:1997812/150‑1 (MQ=255)
                                                                                                    gTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGTTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACtttt                                                 >  1:2919722/1‑150 (MQ=255)
                                                                                                       aaTTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTgag                                              <  1:1549674/150‑3 (MQ=255)
                                                                                                         tttACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGAGcc                                            <  2:2435737/150‑1 (MQ=255)
                                                                                                          ttACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGAGCCa                                           <  1:967570/150‑1 (MQ=255)
                                                                                                              gCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGAGCCAATTa                                       >  2:4351193/1‑150 (MQ=255)
                                                                                                                atatTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATg                                     <  1:772371/150‑1 (MQ=255)
                                                                                                                atatTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATg                                     <  1:3751251/150‑1 (MQ=255)
                                                                                                                  atTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGAGCCAATTATGga                                   <  2:1455748/150‑1 (MQ=255)
                                                                                                                  atTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGAGCCAATTATGga                                   <  1:37317/150‑1 (MQ=255)
                                                                                                                   ttCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGAGCCAATTATGgag                                  <  1:3412527/150‑1 (MQ=255)
                                                                                                                     cTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGaa                                >  2:2279435/1‑150 (MQ=255)
                                                                                                                        tatCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGTGaaaaa                             <  2:2570013/150‑1 (MQ=255)
                                                                                                                        tatCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGaaaaa                             <  2:6603/150‑1 (MQ=255)
                                                                                                                         atCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGAGCCAATTATGGAGaaaaaa                            >  2:3886718/1‑150 (MQ=255)
                                                                                                                         atCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGAGCCAATTATGGAGaaaaaa                            >  1:132258/1‑150 (MQ=255)
                                                                                                                          tccttcctTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGaaaaaaa                           >  2:2196219/1‑150 (MQ=255)
                                                                                                                            cttccttTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGAGCCAATTATGGAGAAAAAAAGa                         >  1:4499189/1‑150 (MQ=255)
                                                                                                                              tccttTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGAGCCAATTATGGAGAAAAAAAGaaa                       <  2:3761209/150‑1 (MQ=255)
                                                                                                                               ccttTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGaaaa                      <  1:1375234/150‑1 (MQ=255)
                                                                                                                                cttTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAAt                     >  1:2438319/1‑150 (MQ=255)
                                                                                                                                 tttCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAAtt                    <  1:659197/150‑1 (MQ=255)
                                                                                                                                 tttCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAAtt                    >  2:3932082/1‑150 (MQ=255)
                                                                                                                                  ttCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAAttt                   <  2:2609797/150‑1 (MQ=255)
                                                                                                                                   tCACTTATTTATTACGGTTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTg                  <  1:1994291/150‑1 (MQ=255)
                                                                                                                                     aGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGAt                <  2:2326810/150‑1 (MQ=255)
                                                                                                                                      gTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATg               <  1:197193/150‑1 (MQ=255)
                                                                                                                                          tttattACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGATATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAGAGCAATGATAGTGTAAAAAACACTTTTGAGCCAATTATGGAGAAAAAAAGAAAATTTGATGgaga           <  2:467097/150‑1 (MQ=255)
                                                                                                                                          tttattACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGgaga           >  1:4330247/1‑150 (MQ=255)
                                                                                                                                            tattaCGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGt         <  1:4460940/150‑1 (MQ=255)
                                                                                                                                            tattaCGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGt         >  2:3049062/1‑150 (MQ=255)
                                                                                                                                             attaCGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGAGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGtg        >  1:3641219/1‑150 (MQ=255)
                                                                                                                                              ttaCGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGtga       <  1:4193146/150‑1 (MQ=255)
                                                                                                                                              ttaCGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTGTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGAGCCAATAATGGAGAAAAAAAGAAAATTTGATGGAGAGtga       >  1:2895702/1‑150 (MQ=255)
                                                                                                                                               taCGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGAGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGtgat      >  2:548499/1‑150 (MQ=255)
                                                                                                                                                aCGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGtgatg     <  2:3941118/150‑1 (MQ=255)
                                                                                                                                                aCGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGtgatg     >  1:4153988/1‑150 (MQ=255)
                                                                                                                                                 cGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGtgatga    >  2:4580106/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAAAGAAAATTTGATGGAGAGTGATGAGA  >  CP009273/1819769‑1820065

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: