Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,888,584 G→T 90.7% I185I (ATC→ATA ygcQ ← putative flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,888,5840GT90.7% 942.7 / 84.7 334I185I (ATC→ATAygcQputative flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (19/12);  new base T (153/150);  total (172/162)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 334 total aligned reads displayed.

TGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTCAGTTGTTGCAT  >  CP009273/2888437‑2888730
                                                                                                                                                   |                                                                                                                                                  
tgcGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTAt                                                                                                                                                  <  2:102925/150‑1 (MQ=255)
  cgcCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTAttt                                                                                                                                                <  1:1610344/150‑1 (MQ=255)
     cAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCt                                                                                                                                             >  1:826482/1‑150 (MQ=255)
     cAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCt                                                                                                                                             >  1:2130876/1‑150 (MQ=255)
           cgcCAACAACAATACACCCCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTGGCCCACTTCGGCTCCCAGCTGTGCTGCCAGCCTAGCGATGTCCTgggtct                                                                                                                                       >  2:4231623/1‑146 (MQ=255)
           cgcCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTAt                                                                                                                                       >  2:4489678/1‑150 (MQ=255)
           cgcCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTAt                                                                                                                                       >  1:26073/1‑150 (MQ=255)
           cgcCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTAt                                                                                                                                       >  1:2330633/1‑150 (MQ=255)
             ccAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCt                                                                                                                                     <  1:4531980/150‑1 (MQ=255)
              caacaacaaTACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTg                                                                                                                                    <  2:581680/150‑1 (MQ=255)
               aacaacaaTACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGc                                                                                                                                   <  1:1282746/150‑1 (MQ=255)
                acaacaaTACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGcc                                                                                                                                  <  2:2000679/150‑1 (MQ=255)
                  aacaaTACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCctc                                                                                                                                >  1:3918389/1‑150 (MQ=255)
                  aacaaTACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCctc                                                                                                                                <  2:1637690/150‑1 (MQ=255)
                  aacaaTACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCctc                                                                                                                                >  1:4349803/1‑150 (MQ=255)
                   acaaTACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCctct                                                                                                                               <  2:61172/150‑1 (MQ=255)
                   acaaTACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCctct                                                                                                                               >  2:2992657/1‑150 (MQ=255)
                    caaTACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCtctc                                                                                                                              >  2:3299594/1‑150 (MQ=255)
                     aaTACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCACTCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTcc                                                                                                                             <  2:1858341/150‑1 (MQ=255)
                     aaTACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTcc                                                                                                                             <  1:398511/150‑1 (MQ=255)
                      aTACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCg                                                                                                                            <  2:77540/150‑1 (MQ=255)
                       tACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGc                                                                                                                           >  2:3337699/1‑150 (MQ=255)
                          acacCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCt                                                                                                                        <  1:888313/150‑1 (MQ=255)
                          acacCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCt                                                                                                                        <  1:237517/150‑1 (MQ=255)
                          acacCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCt                                                                                                                        <  2:1515838/150‑1 (MQ=255)
                           cacCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTg                                                                                                                       <  1:776192/150‑1 (MQ=255)
                            acCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGc                                                                                                                      <  2:1777721/150‑1 (MQ=255)
                            acCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGc                                                                                                                      >  2:3177033/1‑150 (MQ=255)
                              cTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGccc                                                                                                                    >  2:1564413/1‑150 (MQ=255)
                                 aGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCCTAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCccac                                                                                                                 >  2:1777460/1‑150 (MQ=255)
                                         gcagATGTCCGGAAATACCAATCACTTTTTCGGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCAAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCa                                                                                                         <  2:1901691/150‑1 (MQ=255)
                                         gcagATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCa                                                                                                         <  2:4153624/150‑1 (MQ=255)
                                          cagATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCAc                                                                                                        >  1:3400468/1‑150 (MQ=255)
                                            gATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCAccc                                                                                                      >  2:1948166/1‑150 (MQ=255)
                                             aTGTCCGGAAATACCAATCACTTTTTCAGCAACAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCg                                                                                                     <  2:3258951/150‑1 (MQ=255)
                                              tGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTCTTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGt                                                                                                    >  2:479144/1‑150 (MQ=255)
                                                 ccGGACATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCCGCACCCGTCTg                                                                                                 <  1:3725840/150‑1 (MQ=255)
                                                 ccGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTg                                                                                                 >  1:2336558/1‑150 (MQ=255)
                                                  cGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTgg                                                                                                >  2:146845/1‑150 (MQ=255)
                                                   ggAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGc                                                                                               >  2:3059381/1‑150 (MQ=255)
                                                   ggAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGc                                                                                               >  2:3498639/1‑150 (MQ=255)
                                                            aaTCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTaa                                                                                      <  1:2212030/150‑1 (MQ=255)
                                                              tCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATg                                                                                    >  1:4360761/1‑150 (MQ=255)
                                                              tCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATg                                                                                    <  1:3247996/150‑1 (MQ=255)
                                                                aCTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATggg                                                                                  <  1:4081001/150‑1 (MQ=255)
                                                                aCTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATggg                                                                                  <  2:840546/150‑1 (MQ=255)
                                                                aCTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATggg                                                                                  >  1:323433/1‑150 (MQ=255)
                                                                       cAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCggg                                                                           >  2:1619022/1‑150 (MQ=255)
                                                                         gCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTa                                                                         >  1:484742/1‑150 (MQ=255)
                                                                           aTCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTCTTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAAc                                                                       >  2:822517/1‑150 (MQ=255)
                                                                              aaCGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACAtt                                                                    >  1:1607350/1‑150 (MQ=255)
                                                                              aaCGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACAtt                                                                    >  1:3468264/1‑150 (MQ=255)
                                                                               aCGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACAttt                                                                   <  1:1560580/150‑1 (MQ=255)
                                                                                cGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACAtttt                                                                  >  2:2241096/1‑150 (MQ=255)
                                                                                cGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACAtttt                                                                  >  2:1122126/1‑150 (MQ=255)
                                                                                 gCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACAttttt                                                                 >  1:1002409/1‑150 (MQ=255)
                                                                                   cACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCa                                                               >  2:2828880/1‑150 (MQ=255)
                                                                                    aCCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAg                                                              >  2:2729015/1‑150 (MQ=255)
                                                                                      cGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGAt                                                            <  1:1733511/150‑1 (MQ=255)
                                                                                      cGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGAt                                                            <  1:2122319/150‑1 (MQ=255)
                                                                                        ttCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATcc                                                          >  2:4014205/1‑150 (MQ=255)
                                                                                           atcaCTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCctct                                                       <  2:2330633/150‑1 (MQ=255)
                                                                                            tcaCTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCtg                                                      <  2:1124583/150‑1 (MQ=255)
                                                                                             caCTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCtgt                                                     >  2:3745587/1‑150 (MQ=255)
                                                                                              aCTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCtgtg                                                    <  1:4489678/150‑1 (MQ=255)
                                                                                              aCTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCtgtg                                                    <  1:2953098/150‑1 (MQ=255)
                                                                                               cTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGc                                                   <  2:2412364/150‑1 (MQ=255)
                                                                                               cTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGc                                                   >  2:3433307/1‑150 (MQ=255)
                                                                                               cTCGCGCCCGACAGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGc                                                   <  2:1197888/150‑1 (MQ=255)
                                                                                                tCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCt                                                  <  2:3557316/150‑1 (MQ=255)
                                                                                                 cgcgCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTa                                                 <  2:646665/150‑1 (MQ=255)
                                                                                                  gcgcCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTCATGGGTCGCGGGTAACATTTTTCAGATCCACTGTGCaaa                                                >  2:2552999/1‑150 (MQ=255)
                                                                                                      ccGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAAccag                                            >  2:949149/1‑150 (MQ=255)
                                                                                                      ccGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAAccag                                            >  2:4273997/1‑150 (MQ=255)
                                                                                                        gACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAAccagcc                                          <  1:4000843/150‑1 (MQ=255)
                                                                                                          cTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAAccagccag                                        >  2:187949/1‑150 (MQ=255)
                                                                                                                cccACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGa                                  >  1:1659847/1‑150 (MQ=255)
                                                                                                                 ccACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGaa                                 <  1:2100355/150‑1 (MQ=255)
                                                                                                                 ccACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCATGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTATTCAGATCCTCTGTGCTAACCAGACAGTCCGGaa                                 >  2:2046931/1‑150 (MQ=255)
                                                                                                                  cACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAg                                >  2:4333208/1‑150 (MQ=255)
                                                                                                                  cACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAg                                <  1:3295465/150‑1 (MQ=255)
                                                                                                                   aCTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAgc                               <  2:2908113/150‑1 (MQ=255)
                                                                                                                     ttCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAgcgc                             <  2:1091587/150‑1 (MQ=255)
                                                                                                                       cGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCAc                           >  2:1467729/1‑150 (MQ=255)
                                                                                                                       cGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCAc                           >  2:1314958/1‑150 (MQ=255)
                                                                                                                       cGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCAc                           >  1:1806554/1‑150 (MQ=255)
                                                                                                                        ggCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCAcc                          <  1:2272866/150‑1 (MQ=255)
                                                                                                                        ggCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCAcc                          <  1:146845/150‑1 (MQ=255)
                                                                                                                          cTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCg                        >  1:159960/1‑150 (MQ=255)
                                                                                                                           tCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCgg                       <  1:654661/150‑1 (MQ=255)
                                                                                                                           tCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCgg                       >  2:544418/1‑150 (MQ=255)
                                                                                                                               aGCTTTTCTGCCAGCAGAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACa                   <  2:3918389/150‑1 (MQ=255)
                                                                                                                                  ttttCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATg                >  1:1150398/1‑150 (MQ=255)
                                                                                                                                    ttCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGtt              <  2:3078926/150‑1 (MQ=255)
                                                                                                                                     tCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTc             <  2:1601635/150‑1 (MQ=255)
                                                                                                                                      cTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTCa            >  1:2474461/1‑150 (MQ=255)
                                                                                                                                        gCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTCAgt          >  1:2847144/1‑150 (MQ=255)
                                                                                                                                        gCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTCAgt          >  1:2687201/1‑150 (MQ=255)
                                                                                                                                               tAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTCAGTTGTTGCa   >  2:3619477/1‑150 (MQ=255)
                                                                                                                                               tAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTCAGTTGTTGCa   >  1:717095/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTCAGTTGTTGCAT  >  CP009273/2888437‑2888730

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: