Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,244,663 G→T 100% H578N (CAC→AAC)  plsB ← glycerol‑3‑phosphate O‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,244,6630GT95.3% 553.2 / 17.7 193H578N (CAC→AAC) plsBglycerol‑3‑phosphate O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (7/2);  new base T (102/82);  total (109/84)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 193 total aligned reads displayed.

ATTGCCGCCATCAGCGAAGGCAGCACCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGT  >  CP009273/4244518‑4244807
                                                                                                                                                 |                                                                                                                                                
  tGCCGCCATCAGCGAAGGCAGCACCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCg                                                                                                                                              <  1:359492/148‑1 (MQ=255)
  tGCCGCCATCAGCGAAGGCAGCACCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCg                                                                                                                                              <  1:447404/148‑1 (MQ=255)
   gccgccATCAGCGAAGGCAGCACCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGa                                                                                                                                             <  1:2293027/148‑1 (MQ=255)
   gccgccATCAGCGAAGGCAGCACCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGa                                                                                                                                             <  1:2053858/148‑1 (MQ=255)
       ccATCAGCGAAGGCAGCACCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAg                                                                                                                                         <  1:1690890/148‑1 (MQ=255)
        cATCAGCGAAGGCAGCACCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGc                                                                                                                                        >  2:2813894/1‑148 (MQ=255)
         aTCAGCGAAGGCAGCACCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGAGTGTGTCTTTCTCGACTTCAATCTTGTTCATTTGCAGCGCGTGATCGATAAGCt                                                                                                                                       >  2:209968/1‑148 (MQ=255)
          tCAGCGAAGGCAGCACCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTc                                                                                                                                      >  2:3485599/1‑148 (MQ=255)
          tCAGCGAAGGCAGCACCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTc                                                                                                                                      >  1:973803/1‑148 (MQ=255)
          tCAGCGAAGGCAGCACCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTc                                                                                                                                      >  1:2730504/1‑148 (MQ=255)
           cAGCGAAGGCAGCACCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGGTATCGATAAGCTcg                                                                                                                                     >  2:2574478/1‑148 (MQ=255)
            aGCGAAGGCAGCACCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATACGCTcgc                                                                                                                                    >  2:2002983/1‑148 (MQ=255)
            aGCGAAGGCAGCACCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTcgc                                                                                                                                    <  2:505983/148‑1 (MQ=255)
                  ggCAGCACCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTcgctggcgc                                                                                                                              >  1:3067068/1‑148 (MQ=255)
                   gcagcaCCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTcgctggcgct                                                                                                                             >  2:1930808/1‑148 (MQ=255)
                      gcaCCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCgctggcgctggc                                                                                                                          <  2:683898/148‑1 (MQ=255)
                      gcaCCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCgctggcgctggc                                                                                                                          >  1:594027/1‑148 (MQ=255)
                       caCCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCgctggct                                                                                                                         >  1:2973840/1‑148 (MQ=255)
                       caCCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCgctggct                                                                                                                         >  1:2923261/1‑148 (MQ=255)
                       caCCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCgctggct                                                                                                                         >  1:353717/1‑148 (MQ=255)
                       caCCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCgctggct                                                                                                                         >  2:633461/1‑148 (MQ=255)
                         ccAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGa                                                                                                                       >  2:822104/1‑148 (MQ=255)
                          cAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGaa                                                                                                                      >  2:1885784/1‑148 (MQ=255)
                            aCATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAgg                                                                                                                    >  1:1877878/1‑148 (MQ=255)
                             cATCTGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGa                                                                                                                   <  2:3172462/148‑1 (MQ=255)
                              atatGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGaa                                                                                                                  <  1:2831310/148‑1 (MQ=255)
                               tatGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAAc                                                                                                                 >  2:1774621/1‑148 (MQ=255)
                                atGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACg                                                                                                                >  1:2808195/1‑148 (MQ=255)
                                   cgcAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTa                                                                                                             >  2:22219/1‑148 (MQ=255)
                                    gcAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTag                                                                                                            <  2:2047810/148‑1 (MQ=255)
                                    gcAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTag                                                                                                            <  2:2047838/148‑1 (MQ=255)
                                      aaTGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTagag                                                                                                          <  2:2994425/148‑1 (MQ=255)
                                        tgttgtTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTc                                                                                                        <  2:1702614/148‑1 (MQ=255)
                                        tgttgtTGCGATAGTAGGTCATCAGCACCGCGTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTc                                                                                                        <  2:69172/148‑1 (MQ=255)
                                           tgttgCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGt                                                                                                     <  1:955709/148‑1 (MQ=255)
                                              tgCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGa                                                                                                  <  2:2716319/148‑1 (MQ=255)
                                                 gATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTa                                                                                               <  1:2895733/148‑1 (MQ=255)
                                                   tagtagGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAgg                                                                                             >  2:286886/1‑148 (MQ=255)
                                                   tagtagGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTCTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAgg                                                                                             >  2:3401922/1‑148 (MQ=255)
                                                    agtagGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAggg                                                                                            >  2:858757/1‑148 (MQ=255)
                                                    agtagGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAggg                                                                                            <  1:1100593/148‑1 (MQ=255)
                                                      tagGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCa                                                                                          >  2:3205671/1‑148 (MQ=255)
                                                      tagGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCa                                                                                          >  2:3205372/1‑148 (MQ=255)
                                                      tagGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCa                                                                                          >  1:2408474/1‑148 (MQ=255)
                                                         gTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGt                                                                                       >  2:1761498/1‑148 (MQ=255)
                                                         gTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGt                                                                                       >  2:553520/1‑148 (MQ=255)
                                                           catcaGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTg                                                                                     <  2:1683515/148‑1 (MQ=255)
                                                             tcaGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTgcg                                                                                   >  2:1441632/1‑148 (MQ=255)
                                                                gCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCAt                                                                                >  1:2178208/1‑148 (MQ=255)
                                                                gCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCAt                                                                                <  1:1250358/148‑1 (MQ=255)
                                                                gCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCAt                                                                                <  2:1589965/148‑1 (MQ=255)
                                                                 cACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATc                                                                               >  2:13281/1‑148 (MQ=255)
                                                                   ccGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAg                                                                             >  1:3379411/1‑148 (MQ=255)
                                                                    cGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTGTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGa                                                                            <  2:1876749/148‑1 (MQ=255)
                                                                      cTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATc                                                                          >  1:2026120/1‑148 (MQ=255)
                                                                      cTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATc                                                                          >  1:182458/1‑148 (MQ=255)
                                                                       ttGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATcc                                                                         >  1:3166794/1‑148 (MQ=255)
                                                                           tCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGt                                                                     <  1:3444288/148‑1 (MQ=255)
                                                                                cgGCAGAATGATGATGTCGCCGCTTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCag                                                                <  2:1815521/148‑1 (MQ=255)
                                                                                cgGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCag                                                                <  2:1815201/148‑1 (MQ=255)
                                                                                 ggCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCagt                                                               <  2:3248588/148‑1 (MQ=255)
                                                                                   cAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCagttg                                                             >  2:332441/1‑148 (MQ=255)
                                                                                   cAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCagttg                                                             >  1:2510528/1‑148 (MQ=255)
                                                                                     gAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGGAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCagttgag                                                           >  1:1815138/1‑148 (MQ=255)
                                                                                     gAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCagttgag                                                           <  1:121805/148‑1 (MQ=255)
                                                                                     gAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCagttgag                                                           <  1:633461/148‑1 (MQ=255)
                                                                                      aaTGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCagttgagt                                                          >  2:1095966/1‑148 (MQ=255)
                                                                                        tgatgatgTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCagttgagttg                                                        <  2:3185658/148‑1 (MQ=255)
                                                                                        tgatgatgTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCagttgagttg                                                        <  1:1649861/148‑1 (MQ=255)
                                                                                          atgatgTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCt                                                      <  1:2577935/148‑1 (MQ=255)
                                                                                            gatgTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCg                                                    >  1:1665706/1‑148 (MQ=255)
                                                                                             atgTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCgg                                                   <  1:2695975/148‑1 (MQ=255)
                                                                                               gTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATACAGGTAGCAGTTGAGTTGCTCGGtt                                                 >  1:1111457/1‑148 (MQ=255)
                                                                                                     gATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGc                                           >  2:2941929/1‑148 (MQ=255)
                                                                                                      aTTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCt                                          <  2:353717/148‑1 (MQ=255)
                                                                                                           gtCTTTCTCGACTTCAAACTTGTTCATGTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCgg                                     <  1:148566/148‑1 (MQ=255)
                                                                                                            tCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCggg                                    >  1:211753/1‑148 (MQ=255)
                                                                                                             cTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGATCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCTTGGAGTAGGGCACGTTGCGAATCAGATACAGGTAGCAGTTGAGTTGCTAGGTTAAATGCTCGCGGgt                                   >  1:1020771/1‑148 (MQ=255)
                                                                                                              tttCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGgtg                                  <  1:718369/148‑1 (MQ=255)
                                                                                                               ttCTCGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGgtga                                 <  2:2255000/148‑1 (MQ=255)
                                                                                                                tctcGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGgtgag                                <  1:1580622/148‑1 (MQ=255)
                                                                                                                tctcGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGgtgag                                >  2:1171791/1‑148 (MQ=255)
                                                                                                                  tcGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGgtgagtg                              <  2:2814913/148‑1 (MQ=255)
                                                                                                                   cGACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGgtgagtga                             >  1:2811347/1‑148 (MQ=255)
                                                                                                                   cGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGgtgagtga                             >  1:2549607/1‑148 (MQ=255)
                                                                                                                    gACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGtgagtgag                            <  1:1905268/148‑1 (MQ=255)
                                                                                                                    gACTTCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGtgagtgag                            <  1:129553/148‑1 (MQ=255)
                                                                                                                       ttCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAgcgc                         >  1:2410608/1‑148 (MQ=255)
                                                                                                                        tCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCt                        >  2:1974083/1‑148 (MQ=255)
                                                                                                                        tCAAACTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCt                        >  1:3045677/1‑148 (MQ=255)
                                                                                                                            aCTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACg                    >  1:3023830/1‑148 (MQ=255)
                                                                                                                            aCTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTAGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACg                    >  2:2657419/1‑148 (MQ=255)
                                                                                                                             cTTGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGt                   >  2:2378209/1‑148 (MQ=255)
                                                                                                                               tGTTCATTTGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGa                 >  1:3112775/1‑148 (MQ=255)
                                                                                                                                    aTTTGCAGCGCGTTATCGAGAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCa            <  2:2512041/148‑1 (MQ=255)
                                                                                                                                      ttGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGt          <  1:2787557/148‑1 (MQ=255)
                                                                                                                                      ttGCAGCGCGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGt          <  1:1039654/148‑1 (MQ=255)
                                                                                                                                           gcgcGTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAgcgc     >  2:114735/1‑148 (MQ=255)
                                                                                                                                            cgcgTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTTCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCa    >  1:531628/1‑148 (MQ=255)
                                                                                                                                            cgcgTTATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCa    <  2:437973/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ATTGCCGCCATCAGCGAAGGCAGCACCAACATATGCGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGT  >  CP009273/4244518‑4244807

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