Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,242,509 C→T 100% L26F (CTC→TTC)  ubiC → chorismate‑‑pyruvate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,242,5090CT99.4% 1039.2 / ‑2.9 316L26F (CTC→TTC) ubiCchorismate‑‑pyruvate lyase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (157/157);  total (157/159)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 316 total aligned reads displayed.

ATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTAC  >  CP009273/4242362‑4242655
                                                                                                                                                   |                                                                                                                                                  
     cTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTg                                                                                                                                             <  2:1865519/150‑1 (MQ=255)
     cTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTg                                                                                                                                             >  1:1338115/1‑150 (MQ=255)
     cTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTg                                                                                                                                             >  2:4824454/1‑150 (MQ=255)
      tttACATCTCCCGTAAACACTATTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCACCTGTTCGACTgg                                                                                                                                            >  1:1552340/1‑150 (MQ=255)
       ttACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGc                                                                                                                                           >  2:3535946/1‑150 (MQ=255)
         aCATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTg                                                                                                                                         >  2:1952870/1‑150 (MQ=255)
          cATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGt                                                                                                                                        <  2:69171/150‑1 (MQ=255)
             ctcCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGc                                                                                                                                     >  2:1523073/1‑150 (MQ=255)
                ccGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTgg                                                                                                                                  >  1:721818/1‑150 (MQ=255)
                  gTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTggag                                                                                                                                >  1:2011588/1‑150 (MQ=255)
                   tAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTggagg                                                                                                                               <  2:3777793/150‑1 (MQ=255)
                   tAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGATAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTggagg                                                                                                                               >  2:1416613/1‑150 (MQ=255)
                     aaCACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGAt                                                                                                                             >  1:3230352/1‑150 (MQ=255)
                     aaCACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGAt                                                                                                                             <  2:1195798/150‑1 (MQ=255)
                       cacTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTc                                                                                                                           >  2:2152737/1‑150 (MQ=255)
                          ttttCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCAt                                                                                                                        <  2:4959889/150‑1 (MQ=255)
                          ttttCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCAt                                                                                                                        <  2:4312074/150‑1 (MQ=255)
                           tttCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATg                                                                                                                       <  2:2463524/150‑1 (MQ=255)
                            ttCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGa                                                                                                                      >  2:652272/1‑150 (MQ=255)
                            ttCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGa                                                                                                                      <  1:379843/150‑1 (MQ=255)
                              cTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACa                                                                                                                    <  1:4984586/150‑1 (MQ=255)
                                gCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACaaa                                                                                                                  >  2:4599404/1‑150 (MQ=255)
                                 cGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACaaaa                                                                                                                 >  1:33213/1‑150 (MQ=255)
                                  gATACCATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAAc                                                                                                                <  1:4801022/150‑1 (MQ=255)
                                  gATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAAc                                                                                                                <  1:4538702/150‑1 (MQ=255)
                                  gATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAAc                                                                                                                >  1:1647548/1‑150 (MQ=255)
                                     aCAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGtt                                                                                                             <  2:1100031/150‑1 (MQ=255)
                                      cAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGttt                                                                                                            >  1:1070278/1‑150 (MQ=255)
                                        aTGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTg                                                                                                          <  1:906804/150‑1 (MQ=255)
                                          gCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGaa                                                                                                        >  2:61460/1‑150 (MQ=255)
                                           ccTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAc                                                                                                       >  2:5268641/1‑150 (MQ=255)
                                            cTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAca                                                                                                      <  1:2554538/150‑1 (MQ=255)
                                            cTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAca                                                                                                      <  1:2130609/150‑1 (MQ=255)
                                            cTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAca                                                                                                      >  1:2071110/1‑150 (MQ=255)
                                            cTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAca                                                                                                      <  1:152447/150‑1 (MQ=255)
                                              ttACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAcagc                                                                                                    <  1:1268274/150‑1 (MQ=255)
                                                aCGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAcagcag                                                                                                  <  2:956707/150‑1 (MQ=255)
                                                aCGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAcagcag                                                                                                  <  1:2808677/150‑1 (MQ=255)
                                                aCGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAcagcag                                                                                                  <  1:30852/150‑1 (MQ=255)
                                                  gTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTTCGTGCTATGCGATATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAATTCGTTTTGAACAGAAggg                                                                                                <  2:3107239/150‑1 (MQ=255)
                                                  gTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAggg                                                                                                <  2:2971411/150‑1 (MQ=255)
                                                   ttATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGa                                                                                               <  1:652272/150‑1 (MQ=255)
                                                   ttATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGa                                                                                               >  2:1263265/1‑150 (MQ=255)
                                                    tATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGaa                                                                                              <  1:1735971/150‑1 (MQ=255)
                                                    tATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGaa                                                                                              <  2:4263509/150‑1 (MQ=255)
                                                       gTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGaaaaa                                                                                           >  1:3562735/1‑150 (MQ=255)
                                                         aaCGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACg                                                                                         >  1:1925850/1‑150 (MQ=255)
                                                         aaCGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACg                                                                                         >  1:3597903/1‑150 (MQ=255)
                                                             gagaGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTaa                                                                                     <  2:2312120/150‑1 (MQ=255)
                                                              agagTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAg                                                                                    <  2:2783157/150‑1 (MQ=255)
                                                              agagTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAg                                                                                    <  1:340748/150‑1 (MQ=255)
                                                                agTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCg                                                                                  >  2:5067664/1‑150 (MQ=255)
                                                                 gTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGt                                                                                 >  1:2910327/1‑150 (MQ=255)
                                                                     ggCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACg                                                                             >  2:3280957/1‑150 (MQ=255)
                                                                     ggCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACg                                                                             >  1:2280280/1‑150 (MQ=255)
                                                                      gCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACga                                                                            >  1:941950/1‑150 (MQ=255)
                                                                      gCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACga                                                                            <  2:832574/150‑1 (MQ=255)
                                                                      gCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACga                                                                            >  1:2674090/1‑150 (MQ=255)
                                                                        aTGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACgatg                                                                          >  1:3696600/1‑150 (MQ=255)
                                                                        aTGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACgatg                                                                          <  2:5014843/150‑1 (MQ=255)
                                                                              cacCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCcgc                                                                    >  1:3762283/1‑150 (MQ=255)
                                                                                 cccGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGaa                                                                 <  1:3536002/150‑1 (MQ=255)
                                                                                   cgcgTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAgg                                                               <  1:3919108/150‑1 (MQ=255)
                                                                                    gcgTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAggg                                                              <  1:1842931/150‑1 (MQ=255)
                                                                                    gcgTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAggg                                                              >  2:752467/1‑150 (MQ=255)
                                                                                     cgTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGt                                                             >  2:520293/1‑150 (MQ=255)
                                                                                      gTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGtt                                                            >  1:5162864/1‑150 (MQ=255)
                                                                                      gTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGtt                                                            >  2:3168660/1‑150 (MQ=255)
                                                                                      gTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCATGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACTTTTTGAACAGCAGTGAAAAACGGTAAGAGTGACGATGATCCGCGAAGGGtt                                                            >  1:532861/1‑150 (MQ=255)
                                                                                           cGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCg                                                       >  2:4718857/1‑150 (MQ=255)
                                                                                            gCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGa                                                      >  1:4006630/1‑150 (MQ=255)
                                                                                                cTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAg                                                  <  2:4074060/150‑1 (MQ=255)
                                                                                                cTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAg                                                  >  1:795442/1‑150 (MQ=255)
                                                                                                cTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAg                                                  >  1:510357/1‑150 (MQ=255)
                                                                                                  gcgtgcgCTGCGCTATTGGACAGAGATCCCTGCCCAGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGaa                                                <  2:1977051/150‑1 (MQ=255)
                                                                                                   cgtgcgCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAAt                                               <  2:3552560/150‑1 (MQ=255)
                                                                                                       cgctgcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGaaa                                           >  1:1875604/1‑150 (MQ=255)
                                                                                                       cgctgcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGAGGATCCGCGAAGGGTTTGTCGAGCAGAATGaaa                                           >  1:3116677/1‑150 (MQ=255)
                                                                                                        gctgcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAAt                                          <  1:3517962/150‑1 (MQ=255)
                                                                                                          tgcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATcc                                        <  2:3988659/150‑1 (MQ=255)
                                                                                                              ctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCga                                    >  2:1342476/1‑150 (MQ=255)
                                                                                                                    tAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACt                              <  2:1787567/150‑1 (MQ=255)
                                                                                                                    tAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACt                              <  1:4410353/150‑1 (MQ=255)
                                                                                                                       agagATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGcc                           <  2:4449557/150‑1 (MQ=255)
                                                                                                                        gagaTCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCg                          <  1:3946725/150‑1 (MQ=255)
                                                                                                                          gaTCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGACCTGCCgct                        >  2:5325204/1‑150 (MQ=255)
                                                                                                                          gaTCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCgct                        >  2:2245134/1‑150 (MQ=255)
                                                                                                                          gaTCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCgct                        <  2:2901605/150‑1 (MQ=255)
                                                                                                                              cctgccctgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGctgctg                    >  1:242024/1‑150 (MQ=255)
                                                                                                                                cgccctgGATCCGCAACTGTTCGTCTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGTTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGcc                  <  1:1424552/149‑1 (MQ=255)
                                                                                                                                   cctgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGaa               >  1:4485957/1‑150 (MQ=255)
                                                                                                                                    ctgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGaaa              >  2:3836089/1‑150 (MQ=255)
                                                                                                                                     tgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAag             >  2:2707520/1‑150 (MQ=255)
                                                                                                                                       gATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAagag           <  1:1775340/150‑1 (MQ=255)
                                                                                                                                        aTCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGt          <  2:2592237/150‑1 (MQ=255)
                                                                                                                                           cGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGtctc       <  1:3020023/150‑1 (MQ=255)
                                                                                                                                             cAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGt     >  1:4187019/1‑150 (MQ=255)
                                                                                                                                             cAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGt     >  2:3035969/1‑150 (MQ=255)
                                                                                                                                               aCTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTa   <  2:3597903/150‑1 (MQ=255)
                                                                                                                                               aCTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTa   >  1:465610/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTAC  >  CP009273/4242362‑4242655

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: