Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP012306524,2970CT23.0% 503.3 / 59.0 204M1M (ATG→ATA) †citDcitrate lyase, acyl carrier (gamma) subunit
Reads supporting (aligned to +/- strand):  ref base C (39/118);  new base T (47/0);  total (86/118)
Fisher's exact test for biased strand distribution p-value = 1.05e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 231 total aligned reads displayed.

TCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGT  >  AP012306/524148‑524446
                                                                                                                                                     |                                                                                                                                                      
tCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTc                                                                                                                                                        <  1:1957904/150‑1 (MQ=255)
    gAATGGTGGTGCGAATTGCATCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCAttt                                                                                                                                                    <  1:2219181/150‑1 (MQ=255)
     aaTGGTGGTGCGAATTGCATCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCAtttt                                                                                                                                                   <  1:2950997/150‑1 (MQ=255)
     aaTGGTGGTGCGAATTGCATCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTTTTTTTAtttt                                                                                                                                                   >  2:3132379/1‑150 (MQ=255)
      aTGGTGGTGCGAATTGCATCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCAttttt                                                                                                                                                  >  1:1439989/1‑150 (MQ=255)
       tggtggTGCGAATTGCATCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTc                                                                                                                                                 <  1:625619/150‑1 (MQ=255)
         gtggtgCGAATTGCATCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTTATTTTTTAc                                                                                                                                               >  1:3333295/1‑150 (MQ=255)
         gtggtgCGAATTGCATCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTTTTTTTAtttttttc                                                                                                                                               >  2:2418932/1‑147 (MQ=255)
           ggtgCGAATTGCATCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCt                                                                                                                                             <  1:198960/150‑1 (MQ=255)
           ggtgCGAATTGCATCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTTTTTTTATTTTTTACCt                                                                                                                                             >  1:3019161/1‑150 (MQ=255)
               cGAATTGCATCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGttt                                                                                                                                         <  2:1014241/150‑1 (MQ=255)
               cGAATTGCATCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTTTTTTTATTTTTTACCTGgtt                                                                                                                                         >  2:1737275/1‑150 (MQ=255)
                gAATTGCATCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTTCTTTTTTACCCGTTtc                                                                                                                                        >  1:1288957/1‑150 (MQ=255)
                  aTTGCATCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTtctc                                                                                                                                      <  1:780175/150‑1 (MQ=255)
                   ttGCATCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCCGGGTTTTTTTTATTTTTTACCTCTTTCTCa                                                                                                                                     >  2:1983429/1‑150 (MQ=255)
                     gCATCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGGGGGGTTTTTTTCATTTTTTAGCTCTTTCTTATg                                                                                                                                   >  2:3742691/1‑150 (MQ=255)
                     gCATCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTTTTTTTATTTTTTACCTGTTTTTCTTg                                                                                                                                   >  2:3612988/1‑150 (MQ=255)
                      cATCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTTTTTTTATTTTTTCCCTGTTtctccttc                                                                                                                                  >  2:3282987/1‑146 (MQ=255)
                       aTCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCg                                                                                                                                 <  1:1910742/150‑1 (MQ=255)
                       aTCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCg                                                                                                                                 <  1:3851461/150‑1 (MQ=255)
                       aTCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTTTTTTCATTTTTTACCCCTTTCTTCTGCg                                                                                                                                 >  2:2658503/1‑150 (MQ=255)
                         cGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTTTTTTTATTTTTTTCCTTTTTTTCCTTCggg                                                                                                                               >  1:1329509/1‑150 (MQ=255)
                         cGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCCGGGTTATTTTTTTTTTTTACCCTTTTTTTTTGCggg                                                                                                                               >  2:866870/1‑150 (MQ=255)
                          gCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCgggg                                                                                                                              <  2:1568084/150‑1 (MQ=255)
                          gCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCGGGTTTTTTTTCATTTTTTTCCTCTTTTTCCTGCgggg                                                                                                                              >  2:3551636/1‑150 (MQ=255)
                            cAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGt                                                                                                                            <  1:2618857/150‑1 (MQ=255)
                            cAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCCGGTTTATTTTTATTTTTTACCCGTTTTTTAAGCgggggg                                                                                                                            >  2:638462/1‑149 (MQ=255)
                             aaaCTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTc                                                                                                                           <  1:1001063/150‑1 (MQ=255)
                             aaaCTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTc                                                                                                                           <  1:4056657/150‑1 (MQ=255)
                               aCTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTTATTTTTTTCCTGTTTTTTCTTCggggggggt                                                                                                                         >  2:1042748/1‑146 (MQ=255)
                                cTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCttt                                                                                                                        <  2:1442897/150‑1 (MQ=255)
                                 tGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCtttt                                                                                                                       >  1:1873094/1‑150 (MQ=255)
                                  gTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTg                                                                                                                      <  1:2062614/150‑1 (MQ=255)
                                  gTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTg                                                                                                                      <  1:2370133/150‑1 (MQ=255)
                                   tttCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTTATTTTTTTCCTGTTTCTCATGCGGGGGGCtttttt                                                                                                                     >  1:2162373/1‑148 (MQ=255)
                                    ttCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGAc                                                                                                                    <  2:2168803/150‑1 (MQ=255)
                                    ttCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGAc                                                                                                                    <  1:4089558/150‑1 (MQ=255)
                                    ttCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTTATTTTTTACCTGTTTTTCATGCGGGGGGCTTTTTAc                                                                                                                    >  1:400271/1‑150 (MQ=255)
                                     tctcAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACg                                                                                                                   <  1:2433121/150‑1 (MQ=255)
                                       tcAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTTTTTTTATTTTTTACCCGTTTTTCCTTCGGGGGGCTTTTTACGAg                                                                                                                 >  1:3629043/1‑150 (MQ=255)
                                        cAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTTATTTTTTACCTGTTTCTCCTTCGGGGGGCTTTTGACGAGc                                                                                                                >  1:526179/1‑150 (MQ=255)
                                         aaCGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTTATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCt                                                                                                               >  1:3495101/1‑150 (MQ=255)
                                         aaCGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCt                                                                                                               <  1:312937/150‑1 (MQ=255)
                                             ctgctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGGTGGGTTTTTTTCATTTTTTCCCCGTTTCTCATGGGGGGGGCTTTTTGCGAGATGCcg                                                                                                           >  2:3809956/1‑150 (MQ=255)
                                               gctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTTATTTTTTACCCGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCcgcg                                                                                                         >  1:1013200/1‑150 (MQ=255)
                                               gctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTTATTTTTTACCTGTTTCTCATGCGGGGGGGTTTTGGGGAGGTGCcgcg                                                                                                         >  2:2360049/1‑150 (MQ=255)
                                               gctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTTACCTGTTTTTCATTCGGGGGGCTTTTTACGAGCTGCcgcg                                                                                                         >  2:2645391/1‑150 (MQ=255)
                                                ctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCCGGTTTTTTTTCCTTTTTCACCTGCTTTTCTTCCGGGGGGCTTTTGATGAGCTGCCGCGt                                                                                                        >  2:3747282/1‑150 (MQ=255)
                                                 tATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTTATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTc                                                                                                       >  1:3473897/1‑150 (MQ=255)
                                                 tATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTc                                                                                                       <  1:544527/150‑1 (MQ=255)
                                                     gATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTgg                                                                                                   <  2:835465/150‑1 (MQ=255)
                                                      aTTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTTTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGc                                                                                                  >  1:2563373/1‑150 (MQ=255)
                                                      aTTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGc                                                                                                  <  2:1636544/150‑1 (MQ=255)
                                                         tGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCggg                                                                                               <  1:2256081/150‑1 (MQ=255)
                                                           cAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTTTTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCggggg                                                                                             >  1:2352401/1‑148 (MQ=255)
                                                             ggTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTTTTTTTATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTTGCGAGCTGCCGCGTCCTGGCggggggg                                                                                           >  1:3580382/1‑146 (MQ=255)
                                                                cGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTc                                                                                        <  2:1161430/150‑1 (MQ=255)
                                                                 gATATACTGCGTATCGAGTGGGGCGATGCGTATCATCACCTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGGTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCa                                                                                       <  2:2572865/150‑1 (MQ=255)
                                                                        tGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGt                                                                               <  2:1964895/150‑1 (MQ=255)
                                                                         gCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTTATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGGGGTCCAG‑CAgggt                                                                              >  1:681300/1‑148 (MQ=255)
                                                                            tATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTg                                                                           <  1:373173/150‑1 (MQ=255)
                                                                              tCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCa                                                                         <  2:1213987/150‑1 (MQ=255)
                                                                               cGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCaa                                                                        <  2:3497465/150‑1 (MQ=255)
                                                                                   tGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAaataa                                                                    <  1:1064459/150‑1 (MQ=255)
                                                                                      ggCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGc                                                                 <  2:2297848/150‑1 (MQ=255)
                                                                                            gCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTTTTTTTTTTTTTTACCTGTTTCTCCTGCCGGGGGCTTTTTTCGAGCTGCCGCCTCCTGGCGGGGGGGCCCCAGCCAGGTTTTCCAAA‑AAAGCAGCggt                                                           >  2:778565/1‑148 (MQ=255)
                                                                                             cGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGa                                                          <  1:3670477/150‑1 (MQ=255)
                                                                                              gTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGAc                                                         <  2:2689168/150‑1 (MQ=255)
                                                                                                    tcaCATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAgg                                                   <  2:858469/150‑1 (MQ=255)
                                                                                                         tCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCa                                              <  1:1693977/150‑1 (MQ=255)
                                                                                                          ccccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAg                                             <  2:2238876/150‑1 (MQ=255)
                                                                                                           cccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgc                                            <  2:1708784/150‑1 (MQ=255)
                                                                                                             cAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgcgg                                          <  2:1157947/150‑1 (MQ=255)
                                                                                                                 cTCAAGGGTGCCTGCAACGGCGGGCTGGTTTTTTTTTATTTTTTACCTGTTTTTCATGCGGGGGGCTTTTTTCGGGCTGCCGCGTCCTTGCGGGGGGGGTCCAG‑GAGGGTTTGGTAAAAAAACCGCGGGGCTGCAGGGGCCAGCGGggcg                                      >  2:3080445/1‑150 (MQ=255)
                                                                                                                  tCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGa                                     <  1:2694395/150‑1 (MQ=255)
                                                                                                                   cAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGAt                                    <  2:728263/150‑1 (MQ=255)
                                                                                                                   cAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGAt                                    <  1:1249891/150‑1 (MQ=255)
                                                                                                                    aaGGGTGCCTGCAACGGCGGGCTGGTTTTTTTTTATTTTTTTCCTCTTTCTCCTGCGGGGGGCTTTTTTCGAGCTGCCGCGTCCTGGCGGGGGGGGTCCCA‑CAAGGTCTGCAAAAAAAATAGCGCGGCTGCCGGGGGCAGCGGcgcgcca                                   >  2:4256202/1‑147 (MQ=255)
                                                                                                                    aaGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATa                                   <  1:1101240/150‑1 (MQ=255)
                                                                                                                     aGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAg                                  <  2:4429/150‑1 (MQ=255)
                                                                                                                      gggTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGc                                 <  2:3359444/150‑1 (MQ=255)
                                                                                                                              gCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGAt                         <  2:4304622/150‑1 (MQ=255)
                                                                                                                                aaCGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCg                       <  1:894261/150‑1 (MQ=255)
                                                                                                                                 aCGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGt                      <  2:1941597/150‑1 (MQ=255)
                                                                                                                                    gcggGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGtttt                   >  1:2537839/1‑150 (MQ=255)
                                                                                                                                    gcggGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGtttt                   <  2:2942744/150‑1 (MQ=255)
                                                                                                                                     cggGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGttttt                  >  1:1989960/1‑150 (MQ=255)
                                                                                                                                      gggCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTc                 <  2:2208070/150‑1 (MQ=255)
                                                                                                                                      gggCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGGTTTTc                 >  1:1368084/1‑150 (MQ=255)
                                                                                                                                        gCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGc               >  1:1579193/1‑150 (MQ=255)
                                                                                                                                          tGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGGCCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCca             >  2:4149884/1‑150 (MQ=255)
                                                                                                                                           ggTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcag            <  1:3882222/150‑1 (MQ=255)
                                                                                                                                             tttATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagca          <  2:1846106/150‑1 (MQ=255)
                                                                                                                                              ttATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagcag         >  2:110092/1‑150 (MQ=255)
                                                                                                                                              ttATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagcag         >  2:3828257/1‑150 (MQ=255)
                                                                                                                                               tATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGt        >  1:462474/1‑150 (MQ=255)
                                                                                                                                                aTTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGtt       >  2:530495/1‑150 (MQ=255)
                                                                                                                                                aTTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGtt       >  2:2955709/1‑150 (MQ=255)
                                                                                                                                                  tttCATTTTTCACCTGTTTCTCATACGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGa     <  2:584112/150‑1 (MQ=255)
                                                                                                                                                     cATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATACAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGt  >  2:1549941/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTTCTCAACGCTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAG‑CAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGT  >  AP012306/524148‑524446

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑