Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,994,247 A→T 100% V57E (GTG→GAG)  yigE ← DUF2233 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,994,2470AT99.0% 657.1 / ‑2.4 199V57E (GTG→GAG) yigEDUF2233 family protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (91/106);  total (92/107)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 199 total aligned reads displayed.

GATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCT  >  CP009273/3994102‑3994390
                                                                                                                                                 |                                                                                                                                               
  tGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCg                                                                                                                                             <  1:2687806/148‑1 (MQ=255)
     aCAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGctc                                                                                                                                          <  2:1733442/148‑1 (MQ=255)
       aaaCCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTg                                                                                                                                        <  1:1633621/148‑1 (MQ=255)
        aaCCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGt                                                                                                                                       >  2:3553234/1‑148 (MQ=255)
         aCCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGtt                                                                                                                                      <  2:3757776/148‑1 (MQ=255)
         aCCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGtt                                                                                                                                      >  1:3280430/1‑148 (MQ=255)
         aCCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGtt                                                                                                                                      >  2:1205647/1‑148 (MQ=255)
         aCCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGtt                                                                                                                                      >  1:2201457/1‑148 (MQ=255)
         aCCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGtt                                                                                                                                      <  2:3396280/148‑1 (MQ=255)
          ccGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGttt                                                                                                                                     >  1:597709/1‑148 (MQ=255)
           cGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATTTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTg                                                                                                                                    >  2:3351851/1‑148 (MQ=255)
           cGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTg                                                                                                                                    <  1:608657/148‑1 (MQ=255)
              gcggcgCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAgg                                                                                                                                 <  1:600504/148‑1 (MQ=255)
                  cgcATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTa                                                                                                                             <  2:1557671/148‑1 (MQ=255)
                  cgcATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTa                                                                                                                             <  1:1956998/148‑1 (MQ=255)
                    cATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAAc                                                                                                                           >  2:652150/1‑148 (MQ=255)
                    cATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAAc                                                                                                                           >  2:652639/1‑148 (MQ=255)
                    cATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAAc                                                                                                                           <  1:725628/148‑1 (MQ=255)
                    cATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACAATTAATATAAGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAATACATTTTCTCACGCTCAGTTTGAGGATTAAc                                                                                                                           >  1:3602004/1‑148 (MQ=255)
                     aTAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACg                                                                                                                          <  2:2529446/148‑1 (MQ=255)
                      tAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACgg                                                                                                                         <  1:1894209/148‑1 (MQ=255)
                      tAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACgg                                                                                                                         >  1:2131341/1‑148 (MQ=255)
                         cTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGtat                                                                                                                      >  1:2448658/1‑148 (MQ=255)
                           ttCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATAc                                                                                                                    >  2:3581290/1‑148 (MQ=255)
                           ttCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATAc                                                                                                                    >  2:1638517/1‑148 (MQ=255)
                             catcatAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGc                                                                                                                  <  2:24149/148‑1 (MQ=255)
                             catcatAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGc                                                                                                                  >  1:290063/1‑148 (MQ=255)
                              atcatAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGcc                                                                                                                 <  1:1054894/148‑1 (MQ=255)
                              atcatAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGcc                                                                                                                 <  1:1054888/148‑1 (MQ=255)
                                catAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGTCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTg                                                                                                               >  1:1924249/1‑148 (MQ=255)
                                  tAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTa                                                                                                             >  2:2753903/1‑148 (MQ=255)
                                  tAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTa                                                                                                             <  2:1650121/148‑1 (MQ=255)
                                     aTGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACgg                                                                                                          >  1:3033649/1‑148 (MQ=255)
                                             gTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTc                                                                                                  <  2:2533988/148‑1 (MQ=255)
                                              ttCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCg                                                                                                 >  1:851850/1‑148 (MQ=255)
                                               tCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCAAGTAACGTTCCCCACGCTTAGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCgg                                                                                                <  2:867055/148‑1 (MQ=255)
                                                cATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGa                                                                                               <  1:774791/148‑1 (MQ=255)
                                                cATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGa                                                                                               <  2:2466563/148‑1 (MQ=255)
                                                cATCGCCATCTGCACCTGACCCTGACTATTAAGATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGa                                                                                               <  2:2466616/148‑1 (MQ=255)
                                                 aTCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTACCGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGAt                                                                                              <  1:1101584/148‑1 (MQ=255)
                                                 aTCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGAt                                                                                              <  2:564158/148‑1 (MQ=255)
                                                      cATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTgag                                                                                         >  1:3633928/1‑148 (MQ=255)
                                                       aTCTTCACCTGACACTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCTCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTgaga                                                                                        <  1:22597/148‑1 (MQ=255)
                                                       aTCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTgaga                                                                                        >  2:1058215/1‑148 (MQ=255)
                                                        tCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGagag                                                                                       >  2:569352/1‑148 (MQ=255)
                                                            cACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTgcg                                                                                   >  1:2377010/1‑148 (MQ=255)
                                                              cctgaccctgacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCa                                                                                 <  1:202508/148‑1 (MQ=255)
                                                                   ccctgacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAatcat                                                                            <  1:3668993/148‑1 (MQ=255)
                                                                   ccctgacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAatcat                                                                            <  1:1952565/148‑1 (MQ=255)
                                                                   ccctgacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAatcat                                                                            >  1:3247566/1‑148 (MQ=255)
                                                                     ctgacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAAtcatca                                                                          >  2:1881161/1‑148 (MQ=255)
                                                                     ctgacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAAtcatca                                                                          <  2:851850/148‑1 (MQ=255)
                                                                       gacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGc                                                                        <  2:4274/148‑1 (MQ=255)
                                                                        acTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCg                                                                       >  2:3016974/1‑148 (MQ=255)
                                                                            ttAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATAGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCaa                                                                   >  1:1180500/1‑148 (MQ=255)
                                                                             tAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAAc                                                                  >  1:2488358/1‑148 (MQ=255)
                                                                               atatCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAg                                                                <  2:2916511/148‑1 (MQ=255)
                                                                                          cagAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGa                                                     >  1:637597/1‑148 (MQ=255)
                                                                                             agCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGt                                                  >  1:570997/1‑148 (MQ=255)
                                                                                             agCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGt                                                  >  2:1910553/1‑148 (MQ=255)
                                                                                               cATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTaa                                                <  1:1225306/148‑1 (MQ=255)
                                                                                                 tGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAgg                                              <  1:1405497/148‑1 (MQ=255)
                                                                                                  gTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGt                                             >  1:2379056/1‑148 (MQ=255)
                                                                                                   tAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTg                                            >  2:503435/1‑148 (MQ=255)
                                                                                                    aaCGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGa                                           <  1:1447350/148‑1 (MQ=255)
                                                                                                    aaCGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGa                                           <  1:267111/148‑1 (MQ=255)
                                                                                                     aCGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAg                                          >  1:3579352/1‑148 (MQ=255)
                                                                                                     aCGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAg                                          <  1:3414048/148‑1 (MQ=255)
                                                                                                     aCGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAg                                          <  1:1638489/148‑1 (MQ=255)
                                                                                                     aCGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAg                                          <  2:2784722/148‑1 (MQ=255)
                                                                                                      cGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCACTCTGTTTGAGGATTAACGGTATACTCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAgg                                         >  1:1741970/1‑148 (MQ=255)
                                                                                                         tCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCg                                      <  2:3323461/148‑1 (MQ=255)
                                                                                                         tCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCg                                      <  1:650281/148‑1 (MQ=255)
                                                                                                         tCCCCACGCTTCGCCATTGGCTTTTTGCCAGGACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCg                                      <  1:2303314/148‑1 (MQ=255)
                                                                                                             cACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGaa                                  >  2:482968/1‑148 (MQ=255)
                                                                                                              aCGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGaaa                                 >  2:888135/1‑148 (MQ=255)
                                                                                                              aCGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGaaa                                 <  2:32212/148‑1 (MQ=255)
                                                                                                               cGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGTGAGaaaa                                <  1:3768379/148‑1 (MQ=255)
                                                                                                                 cTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAAt                              <  1:1912124/148‑1 (MQ=255)
                                                                                                                 cTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAAt                              <  2:3170438/148‑1 (MQ=255)
                                                                                                                 cTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAAt                              <  2:3690336/148‑1 (MQ=255)
                                                                                                                   tCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATcc                            >  1:1739445/1‑148 (MQ=255)
                                                                                                                   tCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATcc                            >  1:20681/1‑148 (MQ=255)
                                                                                                                     gCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGt                          >  2:923200/1‑148 (MQ=255)
                                                                                                                     gCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGt                          >  2:2949182/1‑148 (MQ=255)
                                                                                                                     gCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGt                          <  2:2260192/148‑1 (MQ=255)
                                                                                                                       cATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGttt                        >  1:843689/1‑148 (MQ=255)
                                                                                                                         ttGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTga                      >  1:1261095/1‑148 (MQ=255)
                                                                                                                          tGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTgag                     <  1:2927141/148‑1 (MQ=255)
                                                                                                                          tGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTgag                     <  2:199658/148‑1 (MQ=255)
                                                                                                                                tttGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCaa               >  1:1578223/1‑148 (MQ=255)
                                                                                                                                     cAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGa          >  2:2978946/1‑148 (MQ=255)
                                                                                                                                      aGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGAt         <  1:2429813/148‑1 (MQ=255)
                                                                                                                                      aGTACATTTTCTCCCGCTATGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGAt         <  2:3279793/148‑1 (MQ=255)
                                                                                                                                         aCATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGACAAATCCGTTTGAGATTCAAGGTGATCAt      >  1:2840661/1‑148 (MQ=255)
                                                                                                                                          cATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCAtt     <  1:503435/148‑1 (MQ=255)
                                                                                                                                            ttttCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTcc   <  2:3816652/148‑1 (MQ=255)
                                                                                                                                            ttttCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTcc   <  2:2640830/148‑1 (MQ=255)
                                                                                                                                             tttCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCt  >  1:3798826/1‑148 (MQ=255)
                                                                                                                                             tttCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCt  >  1:972475/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCT  >  CP009273/3994102‑3994390

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: ‑