Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 647,234 C→T 100% E27K (GAA→AAA)  citC ← [citrate [pro‑3S]‑lyase] ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273647,2340CT96.0% 951.1 / 20.0 273E27K (GAA→AAA) citC[citrate [pro‑3S]‑lyase] ligase
Reads supporting (aligned to +/- strand):  ref base C (1/10);  new base T (98/164);  total (99/174)
Fisher's exact test for biased strand distribution p-value = 6.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
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 280 total aligned reads displayed.

TCAGCGCCAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACT  >  CP009273/647087‑647380
                                                                                                                                                   |                                                                                                                                                   
 cAGCGCCAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATg                                                                                                                                                  >  1:3000492/1‑150 (MQ=255)
  aGCGCCAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGc                                                                                                                                                 >  2:4161877/1‑150 (MQ=255)
    cgcCAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAg                                                                                                                                               >  2:3200721/1‑150 (MQ=255)
    cgcCAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAg                                                                                                                                               <  2:2785073/150‑1 (MQ=255)
     gcCAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAgg                                                                                                                                              >  1:1582662/1‑150 (MQ=255)
     gcCAGTCCTTCACCGCGGACGGATTCACTGATAGCAAAGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAgg                                                                                                                                              <  1:4033829/150‑1 (MQ=255)
         gTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGAGAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAAtt                                                                                                                                          >  1:3875272/1‑150 (MQ=255)
           ccTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTgg                                                                                                                                        <  2:3784618/150‑1 (MQ=255)
           ccTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTgg                                                                                                                                        <  2:2201145/150‑1 (MQ=255)
              tCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCg                                                                                                                                     >  1:2161560/1‑150 (MQ=255)
                aCCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGAt                                                                                                                                   <  1:1353626/150‑1 (MQ=255)
                     ggacggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCg                                                                                                                              <  2:4120095/150‑1 (MQ=255)
                      gacggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTTTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGc                                                                                                                             <  1:2703993/150‑1 (MQ=255)
                         ggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAt                                                                                                                         >  2:3758116/1‑150 (MQ=255)
                          gagtCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAtt                                                                                                                        <  2:3727463/147‑1 (MQ=255)
                          gaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAtt                                                                                                                        <  2:4169234/150‑1 (MQ=255)
                          gaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAtt                                                                                                                        <  1:3690637/150‑1 (MQ=255)
                            ttCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAtttt                                                                                                                      <  2:1973635/150‑1 (MQ=255)
                             tCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAttttt                                                                                                                     <  1:3287048/150‑1 (MQ=255)
                             tCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCAttttt                                                                                                                     <  2:4462430/150‑1 (MQ=255)
                              cACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAtttttt                                                                                                                    >  2:3470490/1‑150 (MQ=255)
                              cACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAtttttt                                                                                                                    >  1:4081688/1‑150 (MQ=255)
                              cACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCAATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCAtttttt                                                                                                                    >  2:3519984/1‑150 (MQ=255)
                                 tGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTAt                                                                                                                 <  1:741273/150‑1 (MQ=255)
                                   aTAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTAttt                                                                                                               <  2:1008491/150‑1 (MQ=255)
                                       cAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTg                                                                                                           >  2:1820748/1‑150 (MQ=255)
                                       cAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTg                                                                                                           <  2:1242728/150‑1 (MQ=255)
                                        aaCGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGa                                                                                                          <  2:1061893/150‑1 (MQ=255)
                                          cGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAAc                                                                                                        <  2:1944481/150‑1 (MQ=255)
                                          cGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAAc                                                                                                        <  1:2873142/150‑1 (MQ=255)
                                           gCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACg                                                                                                       <  1:3025573/150‑1 (MQ=255)
                                           gCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACg                                                                                                       >  2:132501/1‑150 (MQ=255)
                                            cATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGt                                                                                                      >  1:4553952/1‑150 (MQ=255)
                                             aTTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGtt                                                                                                     <  2:1361508/150‑1 (MQ=255)
                                             aTTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGtt                                                                                                     <  2:2643913/150‑1 (MQ=255)
                                             aTTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGtt                                                                                                     <  1:3612715/150‑1 (MQ=255)
                                              tttAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGttt                                                                                                    <  1:292275/150‑1 (MQ=255)
                                              tttAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGttt                                                                                                    <  1:1739151/150‑1 (MQ=255)
                                               ttAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGtttt                                                                                                   >  1:2502901/1‑150 (MQ=255)
                                               ttAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGtttt                                                                                                   <  2:1023914/150‑1 (MQ=255)
                                                taataaTATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTa                                                                                                  >  1:4149958/1‑150 (MQ=255)
                                                        tttCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTg                                                                                          <  1:3198772/150‑1 (MQ=255)
                                                        tttCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTg                                                                                          <  2:1670382/150‑1 (MQ=255)
                                                        tttCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTg                                                                                          >  2:1387698/1‑150 (MQ=255)
                                                         ttCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGa                                                                                         <  2:2398719/150‑1 (MQ=255)
                                                          tCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGaa                                                                                        >  2:4322500/1‑150 (MQ=255)
                                                           ccGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGaaa                                                                                       <  1:2237213/150‑1 (MQ=255)
                                                             ggCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAat                                                                                     >  1:679387/1‑150 (MQ=255)
                                                              gCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAata                                                                                    <  2:3314701/150‑1 (MQ=255)
                                                              gCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAata                                                                                    <  1:3681887/150‑1 (MQ=255)
                                                               cAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAatat                                                                                   >  1:2512048/1‑150 (MQ=255)
                                                                aaTTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATc                                                                                  >  1:4338843/1‑150 (MQ=255)
                                                                    ccACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTg                                                                              <  2:3049378/150‑1 (MQ=255)
                                                                     cACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGc                                                                             >  2:4460307/1‑150 (MQ=255)
                                                                      aCCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGcc                                                                            <  2:2934410/150‑1 (MQ=255)
                                                                        cGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGa                                                                          >  2:2559728/1‑150 (MQ=255)
                                                                         gCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGaa                                                                         >  1:1117881/1‑150 (MQ=255)
                                                                          cACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAAc                                                                        <  2:2735245/150‑1 (MQ=255)
                                                                           aCGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACa                                                                       >  1:3355893/1‑150 (MQ=255)
                                                                            cgcgATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACat                                                                      <  2:3620036/150‑1 (MQ=255)
                                                                               gATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACataat                                                                   <  2:3739551/150‑1 (MQ=255)
                                                                                 tAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTTTTTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACataata                                                                  >  2:4527540/1‑150 (MQ=255)
                                                                                  aaGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaa                                                                >  2:67381/1‑150 (MQ=255)
                                                                                  aaGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaa                                                                <  1:125527/150‑1 (MQ=255)
                                                                                   aGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaat                                                               >  2:4196060/1‑150 (MQ=255)
                                                                                    gCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaata                                                              <  1:632054/150‑1 (MQ=255)
                                                                                      tttcCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGt                                                            <  2:2279924/150‑1 (MQ=255)
                                                                                        ttCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTAt                                                          <  2:1478161/150‑1 (MQ=255)
                                                                                          cATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATcc                                                        <  2:932612/150‑1 (MQ=255)
                                                                                           aTCTCGGGTTACTGTAATAAATACTTCGACTGTGTTGTAAACGCTCAAATCATTTTTATGCTGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCt                                                       <  2:451159/150‑1 (MQ=255)
                                                                                           aTCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCt                                                       <  1:1186199/150‑1 (MQ=255)
                                                                                              gcgGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGAGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGaa                                                    <  2:1201989/150‑1 (MQ=255)
                                                                                               cgGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAg                                                   <  2:51106/150‑1 (MQ=255)
                                                                                                  gTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTg                                                >  2:3862047/1‑150 (MQ=255)
                                                                                                    tACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCa                                              <  1:2015444/150‑1 (MQ=255)
                                                                                                          aataaataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtgttgt                                        <  1:1514867/150‑1 (MQ=255)
                                                                                                           ataaataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATgttgtt                                       <  2:3819523/150‑1 (MQ=255)
                                                                                                                 aCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGAt                                 <  2:2659347/150‑1 (MQ=255)
                                                                                                                 aCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGAt                                 <  1:2309339/150‑1 (MQ=255)
                                                                                                                    tCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTg                              <  2:3025862/150‑1 (MQ=255)
                                                                                                                      gACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCa                            >  1:1869162/1‑150 (MQ=255)
                                                                                                                         tgtgGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACg                         <  2:4074274/150‑1 (MQ=255)
                                                                                                                            ggtgTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAAt                      <  1:1021955/150‑1 (MQ=255)
                                                                                                                             gtgtCAACGCTCAAATCATTTTTATGCAGGATTTGGGCGATTTCCGCCATTTTTTTATTTTCTGTACGTTTTACGCGGGTGAAAATATCATTGCCGAAAATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAAtg                     <  1:3798310/150‑1 (MQ=255)
                                                                                                                             gtgtCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAAtg                     >  1:2288039/1‑150 (MQ=255)
                                                                                                                                tCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGTTGCATGTTGTTATCGATTTGCAACGAAtgttg                  <  1:842907/150‑1 (MQ=255)
                                                                                                                                  aaCGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATgttgtt                <  1:4165120/150‑1 (MQ=255)
                                                                                                                                    cGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCa              <  1:342511/150‑1 (MQ=255)
                                                                                                                                    cGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCa              <  1:3935936/150‑1 (MQ=255)
                                                                                                                                    cGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCa              <  1:4496739/150‑1 (MQ=255)
                                                                                                                                       tCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATg           >  2:1834730/1‑150 (MQ=255)
                                                                                                                                        cAACTCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGt          <  1:1647238/150‑1 (MQ=255)
                                                                                                                                        cAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGt          >  2:2075104/1‑150 (MQ=255)
                                                                                                                                          aaTCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTg        <  2:2679204/150‑1 (MQ=255)
                                                                                                                                           aTCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGc       >  2:623583/1‑150 (MQ=255)
                                                                                                                                            tCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCa      <  1:4416134/150‑1 (MQ=255)
                                                                                                                                              aTTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCaaa    <  2:2177932/150‑1 (MQ=255)
                                                                                                                                               ttttTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAAc   <  1:978322/150‑1 (MQ=255)
                                                                                                                                                ttttATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACt  <  2:926974/150‑1 (MQ=255)
                                                                                                                                                ttttATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACt  <  2:927024/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCAGCGCCAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACT  >  CP009273/647087‑647380

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: