Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 372,581 G→T 100% A118A (GCG→GCT yaiL → DUF2058 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273372,5810GT97.1% 743.6 / 9.2 239A118A (GCG→GCTyaiLDUF2058 family protein
Reads supporting (aligned to +/- strand):  ref base G (2/5);  new base T (126/106);  total (128/111)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 239 total aligned reads displayed.

AAGCTCAGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGA  >  CP009273/372436‑372724
                                                                                                                                                 |                                                                                                                                               
aaGCTCAGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTAt                                                                                                                                               <  2:2538011/148‑1 (MQ=255)
 aGCTCAGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTAtt                                                                                                                                              >  1:2625672/1‑148 (MQ=255)
      aGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCg                                                                                                                                         <  1:386420/148‑1 (MQ=255)
      aGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCg                                                                                                                                         <  1:1701107/148‑1 (MQ=255)
        gTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTAACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGACCGCt                                                                                                                                       >  1:3081247/1‑148 (MQ=255)
          gAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCttg                                                                                                                                     <  1:2943598/148‑1 (MQ=255)
           aaGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCttgt                                                                                                                                    <  2:1205679/148‑1 (MQ=255)
           aaGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCttgt                                                                                                                                    <  1:1606759/148‑1 (MQ=255)
            agcagcTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCttgtt                                                                                                                                   >  1:2125363/1‑148 (MQ=255)
             gcagcTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCttgttg                                                                                                                                  <  2:917348/148‑1 (MQ=255)
             gcagcTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCttgttg                                                                                                                                  <  2:1874824/148‑1 (MQ=255)
              cagcTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTgg                                                                                                                                 <  1:621093/148‑1 (MQ=255)
              cagcTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTgg                                                                                                                                 >  2:1844987/1‑148 (MQ=255)
               agcTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGt                                                                                                                                <  1:1348836/148‑1 (MQ=255)
                  ttATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGa                                                                                                                             >  2:144628/1‑148 (MQ=255)
                    aTTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATa                                                                                                                           >  2:2205762/1‑148 (MQ=255)
                     ttGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATaa                                                                                                                          >  1:3449184/1‑148 (MQ=255)
                     ttGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGTCGATaa                                                                                                                          >  2:3437467/1‑148 (MQ=255)
                      tGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAAc                                                                                                                         >  1:3753246/1‑148 (MQ=255)
                      tGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAAc                                                                                                                         >  1:1856320/1‑148 (MQ=255)
                       gAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACa                                                                                                                        <  2:2791765/148‑1 (MQ=255)
                        aaaTGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACaa                                                                                                                       <  2:1597387/148‑1 (MQ=255)
                         aaTGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAAt                                                                                                                      >  1:3261447/1‑148 (MQ=255)
                          aTGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATa                                                                                                                     >  2:2820336/1‑148 (MQ=255)
                          aTGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATa                                                                                                                     >  1:589491/1‑148 (MQ=255)
                          aTGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATa                                                                                                                     >  2:584886/1‑148 (MQ=255)
                           tGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAg                                                                                                                    <  1:434410/148‑1 (MQ=255)
                           tGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAg                                                                                                                    >  2:408030/1‑148 (MQ=255)
                           tGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAg                                                                                                                    <  1:3392806/148‑1 (MQ=255)
                           tGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAg                                                                                                                    <  2:2265387/148‑1 (MQ=255)
                           tGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAg                                                                                                                    >  1:2226055/1‑148 (MQ=255)
                            gAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGc                                                                                                                   >  2:2949891/1‑148 (MQ=255)
                             aaCCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCg                                                                                                                  <  2:1636818/148‑1 (MQ=255)
                               ccGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGAGTAATGGTCGTCTGGCTATTTCCCGCTTGTTGGTCGATAACAATAGCGaa                                                                                                                >  1:1051360/1‑148 (MQ=255)
                                cGAATCACCATTGCCAATGGCGATATTGTTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAg                                                                                                               <  2:3308280/148‑1 (MQ=255)
                                cGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAg                                                                                                               >  2:3369390/1‑148 (MQ=255)
                                  aaTCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGt                                                                                                             >  1:26624/1‑148 (MQ=255)
                                   aTCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTg                                                                                                            >  1:2765235/1‑148 (MQ=255)
                                       ccATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATa                                                                                                        <  2:1078991/148‑1 (MQ=255)
                                       ccATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATa                                                                                                        >  1:2427799/1‑148 (MQ=255)
                                        cATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATAc                                                                                                       >  1:2003095/1‑148 (MQ=255)
                                         aTTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACg                                                                                                      >  2:1198195/1‑148 (MQ=255)
                                         aTTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACg                                                                                                      <  2:1474131/148‑1 (MQ=255)
                                            gCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCa                                                                                                   <  2:2957269/148‑1 (MQ=255)
                                            gCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCa                                                                                                   <  2:1884210/148‑1 (MQ=255)
                                               aagggCGATATTGTTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCatta                                                                                                <  1:3548350/145‑1 (MQ=255)
                                               aaTGGCGATATTGGTTTTAACTTCACTGACCGCAATCTGATTAAGAAGAATTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCatta                                                                                                <  1:1724378/148‑1 (MQ=255)
                                                   gCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATCGCGAAGGTGAATACGCCATTATTcc                                                                                            <  1:2992944/148‑1 (MQ=255)
                                                   gCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTcc                                                                                            <  1:1083828/148‑1 (MQ=255)
                                                      atatTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGc                                                                                         <  1:287606/148‑1 (MQ=255)
                                                       tatTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGcc                                                                                        >  1:3252940/1‑148 (MQ=255)
                                                         ttGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAg                                                                                      >  1:3563669/1‑148 (MQ=255)
                                                         ttGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCACGCCAg                                                                                      >  1:3564485/1‑148 (MQ=255)
                                                           ggTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCg                                                                                    <  2:3629535/148‑1 (MQ=255)
                                                            gTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGt                                                                                   >  2:558113/1‑148 (MQ=255)
                                                            gTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGt                                                                                   >  1:2364044/1‑148 (MQ=255)
                                                              tttAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTg                                                                                 >  1:3807977/1‑148 (MQ=255)
                                                                  aCTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGa                                                                             >  2:1327184/1‑148 (MQ=255)
                                                                    ttCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATa                                                                           >  2:1982118/1‑148 (MQ=255)
                                                                       aCTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATaaaa                                                                        >  1:3516042/1‑148 (MQ=255)
                                                                       aCTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATaaaa                                                                        <  1:3695510/148‑1 (MQ=255)
                                                                         tGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAAtt                                                                      >  2:1154878/1‑148 (MQ=255)
                                                                          tagtgataTCTGTTTAATAATTTTTTTGTCGATATGCTCACAAAGGGGAATTTTATTTTTGGTCGTCTGGCTTTNGCCGTATTTTTGGTGGTTAAAAATAGCGAAGGTTAATACGCCATTATTCCCGCCAGCGTTGCCGATATAATTg                                                                     <  1:1935301/141‑1 (MQ=255)
                                                                          gACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTg                                                                     <  1:632227/148‑1 (MQ=255)
                                                                            cGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCt                                                                   <  2:1080084/148‑1 (MQ=255)
                                                                               cAATCTGATTAAGATGTTTTTTGTCGTTAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAg                                                                <  1:995902/148‑1 (MQ=255)
                                                                                 aTCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAgcg                                                              >  1:1504852/1‑148 (MQ=255)
                                                                                  tCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAgcgc                                                             >  1:1765860/1‑148 (MQ=255)
                                                                                  tCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAgcgc                                                             >  1:832570/1‑148 (MQ=255)
                                                                                    tGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGa                                                           <  2:3607768/148‑1 (MQ=255)
                                                                                      aTTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATg                                                         >  2:1767397/1‑148 (MQ=255)
                                                                                       ttAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGc                                                        >  1:1007981/1‑148 (MQ=255)
                                                                                          agaagaTTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCcag                                                     >  2:35210/1‑148 (MQ=255)
                                                                                          agaagaTTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGAGAAAATTGCTCAGCGCGATGCcag                                                     >  2:1969254/1‑148 (MQ=255)
                                                                                                     gTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAgtgt                                          <  1:647540/148‑1 (MQ=255)
                                                                                                       cGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCCGTATAGTGTTg                                        >  1:68378/1‑148 (MQ=255)
                                                                                                         aTAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGca                                      <  2:761613/148‑1 (MQ=255)
                                                                                                           aaGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGcaca                                    >  2:1673577/1‑148 (MQ=255)
                                                                                                           aaGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCAGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGcaca                                    >  2:1911370/1‑148 (MQ=255)
                                                                                                           aaGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCACGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGcaca                                    >  2:687543/1‑148 (MQ=255)
                                                                                                             gCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAgc                                  >  1:851544/1‑148 (MQ=255)
                                                                                                             gCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAgc                                  >  1:2525054/1‑148 (MQ=255)
                                                                                                                  cacaGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGGATAGTGTTGCACAGCGCGCt                             >  2:2775193/1‑148 (MQ=255)
                                                                                                                        cgcAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAgcgc                       <  1:912839/148‑1 (MQ=255)
                                                                                                                        cgcAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAgcgc                       <  1:144628/148‑1 (MQ=255)
                                                                                                                            gTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGaa                   <  2:3261344/148‑1 (MQ=255)
                                                                                                                               gATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAagag                >  1:1887603/1‑148 (MQ=255)
                                                                                                                                 ttAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCa              <  1:1656351/148‑1 (MQ=255)
                                                                                                                                 ttAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCa              >  1:2502181/1‑148 (MQ=255)
                                                                                                                                  tAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAg             >  1:1836305/1‑148 (MQ=255)
                                                                                                                                  tAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAg             >  2:2217042/1‑148 (MQ=255)
                                                                                                                                    aTGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGa           >  2:2350590/1‑148 (MQ=255)
                                                                                                                                     tGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGAt          <  1:1005258/148‑1 (MQ=255)
                                                                                                                                      ggTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATg         >  2:3351035/1‑148 (MQ=255)
                                                                                                                                       gtcgtcTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGa        <  2:243218/148‑1 (MQ=255)
                                                                                                                                        tcgtcTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCa              <  1:13713/141‑1 (MQ=255)
                                                                                                                                        tcgtcTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCa              >  2:13713/1‑141 (MQ=255)
                                                                                                                                         cgtcTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAg      <  2:2125363/148‑1 (MQ=255)
                                                                                                                                             tGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAgacga  <  2:3807977/148‑1 (MQ=255)
                                                                                                                                             tGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAgacga  >  2:254678/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
AAGCTCAGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGA  >  CP009273/372436‑372724

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: