Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 322,795 A→C 12.0% L55R (CTA→CGA)  betA ← choline dehydrogenase, a flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273322,7950AC12.0% 402.1 / 55.2 167L55R (CTA→CGA) betAcholine dehydrogenase, a flavoprotein
Reads supporting (aligned to +/- strand):  ref base A (99/48);  new base C (8/12);  total (107/60)
Fisher's exact test for biased strand distribution p-value = 2.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 169 total aligned reads displayed.

AGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGA  >  CP009273/322648‑322942
                                                                                                                                                   |                                                                                                                                                   
 gaTCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAg                                                                                                                                                    >  2:2539342/1‑148 (MQ=255)
    ccAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCgg                                                                                                                                                 >  2:2732040/1‑148 (MQ=255)
    ccAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCgg                                                                                                                                                 <  1:3298431/148‑1 (MQ=255)
      aGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTCGCGGGa                                                                                                                                               >  1:851949/1‑148 (MQ=255)
       gcgcATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTCGCGGGaa                                                                                                                                              <  1:1782295/148‑1 (MQ=255)
       gcgcATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGaa                                                                                                                                              <  2:2092856/148‑1 (MQ=255)
       gcgcATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGaa                                                                                                                                              >  2:237451/1‑148 (MQ=255)
       gcgcATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGaa                                                                                                                                              <  1:1535197/148‑1 (MQ=255)
       gcgcATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGaa                                                                                                                                              >  1:185622/1‑148 (MQ=255)
        cgcATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTCGCGGGAAt                                                                                                                                             <  1:2745661/148‑1 (MQ=255)
           aTTGCCACGGATGTATCACATGCCGTTGATCAGCGAAGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCATGTTCCGTTTCATAGGCCCAGTTGTATCTTTTACCCTGTAGCGGGAATGcc                                                                                                                                          <  1:1047292/148‑1 (MQ=255)
            ttGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCa                                                                                                                                         <  2:3255813/148‑1 (MQ=255)
            ttGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCa                                                                                                                                         <  1:245340/148‑1 (MQ=255)
             tGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAg                                                                                                                                        <  2:1951256/148‑1 (MQ=255)
             tGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAg                                                                                                                                        >  2:74506/1‑148 (MQ=255)
             tGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAg                                                                                                                                        <  2:902283/148‑1 (MQ=255)
             tGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAg                                                                                                                                        <  2:900685/148‑1 (MQ=255)
               ccACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAggg                                                                                                                                      >  2:2002151/1‑148 (MQ=255)
                cACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTCGCGGGAATGCCAGGGc                                                                                                                                     <  2:955407/148‑1 (MQ=255)
                cACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGc                                                                                                                                     <  2:2474692/148‑1 (MQ=255)
                cACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGc                                                                                                                                     <  2:2474696/148‑1 (MQ=255)
                    gATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTCGCGGGAATGCCAGGGCAGCg                                                                                                                                 >  2:2988571/1‑148 (MQ=255)
                    gATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTCGCGGGAATGCCAGGGCAGCg                                                                                                                                 >  2:2989360/1‑148 (MQ=255)
                    gATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCg                                                                                                                                 >  1:451914/1‑148 (MQ=255)
                     aTGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCgg                                                                                                                                >  1:867920/1‑148 (MQ=255)
                         aGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCAt                                                                                                                            <  1:1409251/148‑1 (MQ=255)
                          gCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATc                                                                                                                           >  2:1546771/1‑148 (MQ=255)
                           cacaTGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCt                                                                                                                          >  2:444202/1‑148 (MQ=255)
                            acaTGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTg                                                                                                                         <  2:259004/148‑1 (MQ=255)
                             caTGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTgg                                                                                                                        >  1:1313626/1‑148 (MQ=255)
                             caTGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTgg                                                                                                                        >  1:1314096/1‑148 (MQ=255)
                               tGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGt                                                                                                                      >  2:2630178/1‑148 (MQ=255)
                                   gTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCgg                                                                                                                  >  2:424996/1‑148 (MQ=255)
                                    ttGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGa                                                                                                                 <  2:574452/148‑1 (MQ=255)
                                       aTCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGACTGCCAGGGCAGCGGGCAGCTGGGTGCGGAAGt                                                                                                              >  1:3136419/1‑148 (MQ=255)
                                        tCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTCGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTc                                                                                                             <  1:3303177/148‑1 (MQ=255)
                                        tCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTc                                                                                                             >  1:1209663/1‑148 (MQ=255)
                                          aGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTCGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCaa                                                                                                           <  2:2471730/148‑1 (MQ=255)
                                          aGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCaa                                                                                                           <  1:1756324/148‑1 (MQ=255)
                                          aGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCaa                                                                                                           <  1:2633321/148‑1 (MQ=255)
                                          aGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCaa                                                                                                           <  1:2843023/148‑1 (MQ=255)
                                            cgacgaTCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAg                                                                                                         >  2:236919/1‑148 (MQ=255)
                                                 aTCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATa                                                                                                    >  2:3013320/1‑148 (MQ=255)
                                                 aTCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATa                                                                                                    <  1:606922/148‑1 (MQ=255)
                                                 aTCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATa                                                                                                    <  1:3138866/148‑1 (MQ=255)
                                                  tCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAg                                                                                                   >  1:2065282/1‑148 (MQ=255)
                                                   ccACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGt                                                                                                  >  1:3118026/1‑148 (MQ=255)
                                                   ccACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGt                                                                                                  >  2:1619359/1‑148 (MQ=255)
                                                    cACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTc                                                                                                 <  1:2443519/148‑1 (MQ=255)
                                                      cccAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCg                                                                                               >  2:759032/1‑148 (MQ=255)
                                                      cccAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCg                                                                                               >  2:256879/1‑148 (MQ=255)
                                                         aGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGgc                                                                                            >  1:739506/1‑148 (MQ=255)
                                                          gACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGgcc                                                                                           >  1:610633/1‑148 (MQ=255)
                                                            ccTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGgccgc                                                                                         >  2:673556/1‑148 (MQ=255)
                                                              tttACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGccgccc                                                                                       >  1:2021501/1‑148 (MQ=255)
                                                               ttACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCAGCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGccgcccg                                                                                      >  1:1970913/1‑148 (MQ=255)
                                                               ttACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAAAGGTTGAGGTTACGTTTCATAGGCCCAGTTGTAGCGTTTAACCTGTAGCGGGAATGACAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAATCGATATTCCGGGccgcccg                                                                                      >  1:2921678/1‑148 (MQ=255)
                                                                 aCCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCt                                                                                    >  1:526741/1‑148 (MQ=255)
                                                                 aCCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCt                                                                                    >  1:2189327/1‑148 (MQ=255)
                                                                 aCCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCt                                                                                    >  1:526855/1‑148 (MQ=255)
                                                                                 aTGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTCGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACgg                                                                    <  1:2988571/148‑1 (MQ=255)
                                                                                 aTGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTCGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACgg                                                                    <  1:2989315/148‑1 (MQ=255)
                                                                                     ggcggTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGt                                                                <  1:1908849/148‑1 (MQ=255)
                                                                                     ggcggTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGt                                                                <  1:692740/148‑1 (MQ=255)
                                                                                      gcggTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTCGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTa                                                               >  1:2529258/1‑148 (MQ=255)
                                                                                      gcggTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTCGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTa                                                               >  1:2531148/1‑148 (MQ=255)
                                                                                       cggTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTAt                                                              >  1:1150267/1‑148 (MQ=255)
                                                                                         gTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTc                                                            >  2:971742/1‑148 (MQ=255)
                                                                                         gTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTc                                                            >  2:2769336/1‑148 (MQ=255)
                                                                                           tATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCgg                                                          >  1:2997490/1‑148 (MQ=255)
                                                                                             ttCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGAt                                                        >  2:1889406/1‑148 (MQ=255)
                                                                                             ttCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGAt                                                        >  2:2158935/1‑148 (MQ=255)
                                                                                                aTAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCtt                                                     >  1:1648071/1‑148 (MQ=255)
                                                                                                    aCGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTtcagt                                                 >  1:2975536/1‑148 (MQ=255)
                                                                                                            ggTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACggg                                         >  1:3148119/1‑148 (MQ=255)
                                                                                                               tCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAg                                      <  1:3054614/148‑1 (MQ=255)
                                                                                                                   tttCATAGGCCCAGTTGTAGCGTTTACCCTGTCGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCgag                                  >  1:2794362/1‑148 (MQ=255)
                                                                                                                      cATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAAc                               >  1:934932/1‑148 (MQ=255)
                                                                                                                        tAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGt                             <  2:1540068/148‑1 (MQ=255)
                                                                                                                        tAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGt                             <  1:983680/148‑1 (MQ=255)
                                                                                                                         aGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGtt                            >  1:2820521/1‑148 (MQ=255)
                                                                                                                          ggCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTg                           <  2:1356272/148‑1 (MQ=255)
                                                                                                                          ggCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTg                           <  2:1356258/148‑1 (MQ=255)
                                                                                                                          ggCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTg                           >  2:1933311/1‑148 (MQ=255)
                                                                                                                              cAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCgg                       <  1:1889406/148‑1 (MQ=255)
                                                                                                                              cAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGAAGCGAGAACGTTGCCgg                       >  2:147649/1‑148 (MQ=255)
                                                                                                                                 ttGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTg                    <  1:567136/148‑1 (MQ=255)
                                                                                                                                 ttGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTg                    >  1:737913/1‑148 (MQ=255)
                                                                                                                                  tGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGa                   >  1:2511116/1‑148 (MQ=255)
                                                                                                                                     aGCGTTTACCCTGTCGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGcc                <  2:572449/148‑1 (MQ=255)
                                                                                                                                      gCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCg               >  1:2069784/1‑148 (MQ=255)
                                                                                                                                         tttACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCa            >  2:1876508/1‑148 (MQ=255)
                                                                                                                                         tttACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCa            >  2:420022/1‑148 (MQ=255)
                                                                                                                                          ttACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCAc           >  2:3150339/1‑148 (MQ=255)
                                                                                                                                             cccTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCaa        <  2:847283/148‑1 (MQ=255)
                                                                                                                                                tGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCaataa     <  1:3215844/148‑1 (MQ=255)
                                                                                                                                                tGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCaataa     >  2:1153352/1‑148 (MQ=255)
                                                                                                                                                 gTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCaataat    >  1:2946187/1‑148 (MQ=255)
                                                                                                                                                 gTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCaataat    >  1:434112/1‑148 (MQ=255)
                                                                                                                                                  tAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAatg   <  1:1546771/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGA  >  CP009273/322648‑322942

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