Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,143,605 G→T 8.6% intergenic (‑579/+19) ppc ← / ← argE phosphoenolpyruvate carboxylase/acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,143,6050GT 8.6% 522.9 / 20.3 150intergenic (‑579/+19)ppc/argEphosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base G (75/62);  new base T (11/2);  total (86/64)
Fisher's exact test for biased strand distribution p-value = 4.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 156 total aligned reads displayed.

AACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCG  >  CP009273/4143458‑4143751
                                                                                                                                                   |                                                                                                                                                  
 aCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAgcgc                                                                                                                                                   >  1:1230777/1‑148 (MQ=39)
  cAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGcc                                                                                                                                                  >  1:1368382/1‑148 (MQ=39)
  cAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGcc                                                                                                                                                  <  1:1532426/148‑1 (MQ=39)
   aGTGCCGGAAGCGGCGCGAGCGCCTTATGCGGCCTACTGTTGGGCATCTATTGAGTCACTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCt                                                                                                                                                 <  2:51607/148‑1 (MQ=11)
     tGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTa                                                                                                                                               <  2:3206215/148‑1 (MQ=39)
     tGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTa                                                                                                                                               >  1:2756534/1‑148 (MQ=39)
          gATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCgg                                                                                                                                          >  2:520841/1‑148 (MQ=255)
                gcgcGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCTCTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACg                                                                                                                                    <  2:1107591/148‑1 (MQ=255)
                gcgcGAGCGCCTTATCCGGCCTACGGTTGGACATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACg                                                                                                                                    >  2:3271470/1‑148 (MQ=255)
                 cgcgAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGt                                                                                                                                   <  1:1287413/148‑1 (MQ=40)
                         ccTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGc                                                                                                                           <  2:71905/148‑1 (MQ=255)
                          cTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGcc                                                                                                                          >  1:1764516/1‑148 (MQ=255)
                          cTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGcc                                                                                                                          >  2:1220553/1‑148 (MQ=255)
                          cTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGcc                                                                                                                          >  1:91218/1‑148 (MQ=255)
                                   ccTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAAt                                                                                                                 <  1:975141/148‑1 (MQ=255)
                                   ccTACGGCTGGGCATCGTTCGAGTCACTGTAGGTCGGATAAGATGCGATAGCATAGCATCCGACATTAGTGCGGAACTGCTTTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAAt                                                                                                                 <  1:3173987/148‑1 (MQ=17)
                                     tACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATgg                                                                                                               <  2:2527101/148‑1 (MQ=255)
                                          ttGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAAt                                                                                                          >  2:2249765/1‑148 (MQ=255)
                                          ttGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAAt                                                                                                          >  2:1427026/1‑148 (MQ=255)
                                           tGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAAtt                                                                                                         >  1:3009007/1‑148 (MQ=255)
                                                 tCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTgg                                                                                                   >  1:1776885/1‑148 (MQ=255)
                                                  cGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTggg                                                                                                  >  2:240211/1‑148 (MQ=255)
                                                     ttGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGACTTACCTGGGTTa                                                                                               >  1:580810/1‑148 (MQ=255)
                                                         gTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAg                                                                                           <  1:692991/148‑1 (MQ=255)
                                                         gTCACTGTAGGTCGGATAAGAGGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAg                                                                                           <  1:692988/148‑1 (MQ=255)
                                                          tCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTa                                                                                                                                               <  2:1568218/95‑1 (MQ=40)
                                                           cACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGtt                                                                                         >  2:1863214/1‑148 (MQ=255)
                                                           cACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGtt                                                                                         <  2:2672239/148‑1 (MQ=255)
                                                             cTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTcg                                                                                       >  2:1131702/1‑148 (MQ=255)
                                                                tAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCgg                                                                                    >  1:1962181/1‑148 (MQ=255)
                                                                  ggTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGt                                                                                  <  2:1858516/148‑1 (MQ=255)
                                                                      ggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCTAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGc                                                                              <  1:2528070/148‑1 (MQ=255)
                                                                      ggATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGc                                                                              <  2:670840/148‑1 (MQ=255)
                                                                      ggATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGc                                                                              <  1:1263702/148‑1 (MQ=255)
                                                                        aTAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCtt                                                                            >  2:2284845/1‑148 (MQ=255)
                                                                         tAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTg                                                                           >  2:1447729/1‑148 (MQ=255)
                                                                          aaGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGa                                                                          <  1:1235947/148‑1 (MQ=255)
                                                                          aaGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGa                                                                          >  1:1525355/1‑148 (MQ=255)
                                                                          aaGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGa                                                                          >  1:433583/1‑148 (MQ=255)
                                                                          aaGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGa                                                                          >  2:1416533/1‑148 (MQ=255)
                                                                           aGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGAt                                                                         <  2:530952/148‑1 (MQ=255)
                                                                           aGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGAt                                                                         <  1:1937202/148‑1 (MQ=255)
                                                                            gATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATa                                                                        >  1:2168749/1‑148 (MQ=255)
                                                                            gATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATa                                                                        >  2:3184386/1‑148 (MQ=255)
                                                                             aTGCGCAAGCATCGCCTCCGACATTATTGCGGCACTGGAGTTTGTCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATaa                                                                       <  2:2060578/148‑1 (MQ=255)
                                                                             aTGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATaa                                                                       >  1:1443858/1‑148 (MQ=255)
                                                                             aTGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATaa                                                                       >  1:2841900/1‑148 (MQ=255)
                                                                             aTGCGCAAGCATAGCATCCGACATTATCGCTGCACTGGAGTTTGGAATCAGCGCCGGATGCGTCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATaa                                                                       <  2:2677454/148‑1 (MQ=25)
                                                                               gcgcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAAc                                                                     >  1:586007/1‑148 (MQ=255)
                                                                               gcgcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTATTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAAc                                                                     >  2:52385/1‑148 (MQ=255)
                                                                                 gcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCg                                                                   <  1:2332233/148‑1 (MQ=255)
                                                                                 gcAACCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCg                                                                   <  1:301790/148‑1 (MQ=255)
                                                                                  cAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgt                                                                  <  1:2115788/148‑1 (MQ=255)
                                                                                    aGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgtgt                                                                >  2:1243095/1‑148 (MQ=255)
                                                                                    aGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgtgt                                                                <  2:728460/148‑1 (MQ=255)
                                                                                     gcatcgcatcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGt                                                                                  >  2:2371863/1‑129 (MQ=255)
                                                                                     gcatcgcatcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGtt                                                               >  1:906593/1‑148 (MQ=255)
                                                                                      catcgcatcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGttt                                                              <  2:709874/148‑1 (MQ=255)
                                                                                        tcgcatcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTcc                                                            >  2:580116/1‑148 (MQ=255)
                                                                                         cgcatcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAgcg                                                                                                                                                    <  1:1181766/59‑1 (MQ=35)
                                                                                         cgcatcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAgcg                                                                                                                                                    <  1:1182138/59‑1 (MQ=35)
                                                                                         cgcatcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAgcg                                                                                                                                                    <  1:2198125/59‑1 (MQ=35)
                                                                                            atcCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTa                                                                                                                               >  2:2272457/1‑77 (MQ=17)
                                                                                            atcCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTa                                                                                                                               <  1:2272457/77‑1 (MQ=17)
                                                                                             tcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGata                                                       <  1:329439/148‑1 (MQ=255)
                                                                                             tcCGACATCATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGata                                                       <  1:2098603/148‑1 (MQ=255)
                                                                                              ccGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGatat                                                      >  2:3070943/1‑148 (MQ=255)
                                                                                              ccGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGatat                                                      >  2:2981104/1‑148 (MQ=255)
                                                                                              ccGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGatat                                                      >  2:1124979/1‑148 (MQ=255)
                                                                                               cGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATAtt                                                     <  1:2414769/148‑1 (MQ=255)
                                                                                                 aCATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATtca                                                   >  2:2443701/1‑148 (MQ=255)
                                                                                                 aCATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATtca                                                   >  1:1095310/1‑148 (MQ=255)
                                                                                                   attattGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAATCCGTGTTTCCAGATATtcatc                                                 >  2:88725/1‑148 (MQ=255)
                                                                                                   attattGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATtcatc                                                 >  2:1680718/1‑148 (MQ=255)
                                                                                                    ttattGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATtcatca                                                >  1:219714/1‑148 (MQ=255)
                                                                                                       ttGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGt                                             >  1:2393181/1‑148 (MQ=255)
                                                                                                       gtGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGt                                             <  1:2159998/147‑1 (MQ=255)
                                                                                                        tGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGtt                                            >  2:1034974/1‑148 (MQ=255)
                                                                                                         gCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTtg                                           <  1:520841/148‑1 (MQ=255)
                                                                                                                  gAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTg                                  >  1:3154262/1‑148 (MQ=255)
                                                                                                                  gAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTg                                  >  1:3154359/1‑148 (MQ=255)
                                                                                                                  gAGTTGGGCCACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTg                                  <  1:2769389/148‑1 (MQ=255)
                                                                                                                     tttGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGagt                               >  1:1511838/1‑146 (MQ=255)
                                                                                                                     tttGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGAtt                               >  2:2827999/1‑148 (MQ=255)
                                                                                                                      ttGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTa                              <  1:2443701/148‑1 (MQ=255)
                                                                                                                          cAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAAt                                                                                                                 >  2:3173974/1‑61 (MQ=255)
                                                                                                                          cAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTg                          <  1:1011905/148‑1 (MQ=255)
                                                                                                                          cAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGAGATTCATCAGGTTGATGAGCCTGATTAATTg                          <  2:2646741/148‑1 (MQ=255)
                                                                                                                             cAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGc                       >  2:1723225/1‑148 (MQ=255)
                                                                                                                             cAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGc                       >  2:748088/1‑148 (MQ=255)
                                                                                                                             cAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGc                       <  1:2374513/148‑1 (MQ=255)
                                                                                                                             cAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGc                       <  1:2316793/148‑1 (MQ=255)
                                                                                                                              aGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGcc                      <  2:2970241/148‑1 (MQ=255)
                                                                                                                                cgcCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAg                    >  1:2034970/1‑148 (MQ=255)
                                                                                                                                cgcCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTCGATGAGCCTGATTAATTGAGCCAg                    >  2:580372/1‑148 (MQ=255)
                                                                                                                                       tGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCaa             <  1:2414831/148‑1 (MQ=255)
                                                                                                                                       tGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCaa             <  1:345260/148‑1 (MQ=255)
                                                                                                                                           gcgcGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACAcc         <  2:226671/148‑1 (MQ=255)
                                                                                                                                            cgcgAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCa        <  2:2134408/148‑1 (MQ=255)
                                                                                                                                                aGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGt    >  2:1601297/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCG  >  CP009273/4143458‑4143751

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: