Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,770,281 T→A 100% D146E (GAT→GAA ydiF → putative acetyl‑CoA:acetoacetyl‑CoA transferase: alpha subunit/beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,770,2810TA95.6% 895.9 / 22.5 250D146E (GAT→GAAydiFputative acetyl‑CoA:acetoacetyl‑CoA transferase: alpha subunit/beta subunit
Reads supporting (aligned to +/- strand):  ref base T (3/8);  new base A (141/97);  total (145/105)
Fisher's exact test for biased strand distribution p-value = 5.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 250 total aligned reads displayed.

CAATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGG  >  CP009273/1770135‑1770427
                                                                                                                                                  |                                                                                                                                                  
       cgcgTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGc                                                                                                                                          >  1:4089924/1‑150 (MQ=255)
        gcgTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCa                                                                                                                                         >  1:181785/1‑150 (MQ=255)
         cgTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCaa                                                                                                                                        >  1:4468838/1‑150 (MQ=255)
            aTTTCTGAACTCGCAGAACAAAATAAACTTAGTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAggc                                                                                                                                     <  1:2263468/150‑1 (MQ=255)
            aTTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAggc                                                                                                                                     >  1:1777493/1‑150 (MQ=255)
              ttCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAggcgg                                                                                                                                   >  1:4135079/1‑150 (MQ=255)
               tCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATGAGTGATATTGGCATCGGGACATTTGTCGACCCACGCCAGCAAggcggc                                                                                                                                  >  1:3546594/1‑150 (MQ=255)
                 tGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCaa                                                                                                                                >  1:732206/1‑150 (MQ=255)
                  gAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCTCGCCAGCAAGGCGGCaca                                                                                                                               >  2:2810747/1‑148 (MQ=255)
                  gAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCaaa                                                                                                                               >  2:2810621/1‑150 (MQ=255)
                         cAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACtgaatg                                                                                                                        >  1:2555764/1‑150 (MQ=255)
                          aGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACtgaatga                                                                                                                       >  2:1383083/1‑150 (MQ=255)
                            aaCAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAg                                                                                                                     <  1:974710/150‑1 (MQ=255)
                            aaCAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAg                                                                                                                     <  1:4698759/150‑1 (MQ=255)
                             aCAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGt                                                                                                                    >  2:2333294/1‑150 (MQ=255)
                              cAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTc                                                                                                                   >  1:4177797/1‑150 (MQ=255)
                               aaaataaaatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCa                                                                                                                  <  1:3077340/150‑1 (MQ=255)
                                aaataaaatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCAc                                                                                                                 >  1:1526730/1‑150 (MQ=255)
                                 aataaaatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACt                                                                                                                <  1:2349343/150‑1 (MQ=255)
                                  ataaaatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTa                                                                                                               >  1:789081/1‑150 (MQ=255)
                                    aaaatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTaaa                                                                                                             <  1:1191023/150‑1 (MQ=255)
                                       attattGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAaagaa                                                                                                          >  1:846491/1‑150 (MQ=255)
                                         tattGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAagaaga                                                                                                        >  1:2602815/1‑150 (MQ=255)
                                         tattGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAagaaga                                                                                                        >  1:578080/1‑150 (MQ=255)
                                           ttGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGAcc                                                                                                      <  2:3343121/150‑1 (MQ=255)
                                           ttGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGAcc                                                                                                      >  2:151010/1‑150 (MQ=255)
                                            tGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCt                                                                                                     >  2:1945555/1‑150 (MQ=255)
                                            tGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCt                                                                                                     <  2:2086491/150‑1 (MQ=255)
                                             gCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTg                                                                                                    >  1:2169518/1‑150 (MQ=255)
                                              cTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGa                                                                                                   >  1:3163030/1‑150 (MQ=255)
                                                tataACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGAtt                                                                                                 >  1:1699078/1‑150 (MQ=255)
                                                 ataACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTa                                                                                                <  2:3637264/150‑1 (MQ=255)
                                                 ataACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTa                                                                                                >  2:1695916/1‑150 (MQ=255)
                                                   aaCTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTaaa                                                                                              >  1:3224028/1‑150 (MQ=255)
                                                   aaCTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTaaa                                                                                              >  1:2118394/1‑150 (MQ=255)
                                                    acgaCCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAAc                                                                                             <  2:2547899/147‑1 (MQ=255)
                                                    aCTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAAc                                                                                             <  1:2891601/150‑1 (MQ=255)
                                                    aCTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAAc                                                                                             >  1:2514507/1‑150 (MQ=255)
                                                    aCTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAAc                                                                                             <  2:177494/150‑1 (MQ=255)
                                                        cccACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGt                                                                                         <  1:4011836/150‑1 (MQ=255)
                                                         ccACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTc                                                                                        <  1:4310772/150‑1 (MQ=255)
                                                         ccACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTc                                                                                        >  1:3864822/1‑150 (MQ=255)
                                                          cacaAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCg                                                                                       <  1:550550/150‑1 (MQ=255)
                                                            caAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAg                                                                                     >  1:1560718/1‑150 (MQ=255)
                                                            caAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAg                                                                                     >  1:3295243/1‑150 (MQ=255)
                                                              aGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGtt                                                                                   <  1:1012686/150‑1 (MQ=255)
                                                                 tgtACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGa                                                                                >  2:438108/1‑150 (MQ=255)
                                                                 tgtACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGa                                                                                >  1:3421271/1‑150 (MQ=255)
                                                                  gtACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGAt                                                                               >  1:283542/1‑150 (MQ=255)
                                                                    aCTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATaa                                                                             >  2:3097163/1‑150 (MQ=255)
                                                                     cTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAAc                                                                            <  2:4448293/150‑1 (MQ=255)
                                                                      ttACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACa                                                                           <  2:4189243/150‑1 (MQ=255)
                                                                        acacAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACtaa                                                                         <  1:2772507/150‑1 (MQ=255)
                                                                        acacAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACaaa                                                                         <  1:3030701/150‑1 (MQ=255)
                                                                         cacaAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCTTCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAg                                                                        <  2:1399648/150‑1 (MQ=255)
                                                                              aCCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAatat                                                                   <  2:3470575/150‑1 (MQ=255)
                                                                                 ttACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATAtctc                                                                >  1:1644806/1‑150 (MQ=255)
                                                                                 ttACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATAtctc                                                                >  1:1259298/1‑150 (MQ=255)
                                                                                   aCGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTa                                                              <  2:26022/150‑1 (MQ=255)
                                                                                    cgcgCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTAt                                                             <  1:3612411/150‑1 (MQ=255)
                                                                                       gccgccgcAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTAc                                                          <  1:3256005/150‑1 (MQ=255)
                                                                                            cgcAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGc                                                     >  2:3258206/1‑150 (MQ=255)
                                                                                            cgcAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGc                                                     >  2:4019379/1‑150 (MQ=255)
                                                                                            cgcAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGc                                                     >  2:1632678/1‑150 (MQ=255)
                                                                                             gcAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCg                                                    >  1:2014410/1‑150 (MQ=255)
                                                                                              cAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGa                                                   >  2:4398536/1‑150 (MQ=255)
                                                                                               aGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGAt                                                  <  2:1936074/150‑1 (MQ=255)
                                                                                                    accaGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCTCGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTgcgc                                             <  2:4686700/150‑1 (MQ=255)
                                                                                                    accaGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTgcgc                                             >  1:1939630/1‑150 (MQ=255)
                                                                                                    accaGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTgcgc                                             >  2:2786543/1‑150 (MQ=255)
                                                                                                      caGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCa                                           >  2:598955/1‑150 (MQ=255)
                                                                                                      caGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCa                                           <  1:264242/150‑1 (MQ=255)
                                                                                                      caGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCa                                           <  1:206422/150‑1 (MQ=255)
                                                                                                       aGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAg                                          <  2:1431934/150‑1 (MQ=255)
                                                                                                       aGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAg                                          >  2:4658052/1‑150 (MQ=255)
                                                                                                       aGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAg                                          <  2:3822072/150‑1 (MQ=255)
                                                                                                         ccTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTCAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGat                                        >  2:4554379/1‑150 (MQ=255)
                                                                                                              tattatTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGc                                   <  1:3069453/150‑1 (MQ=255)
                                                                                                              tattatTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGc                                   >  1:3779417/1‑150 (MQ=255)
                                                                                                              tattatTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGc                                   >  1:1388548/1‑150 (MQ=255)
                                                                                                                ttattaGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCt                                 >  1:3214767/1‑150 (MQ=255)
                                                                                                                 tattaGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCtt                                >  2:1612300/1‑150 (MQ=255)
                                                                                                                 tattaGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCtt                                >  2:448674/1‑150 (MQ=255)
                                                                                                                   ttaGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCa                              >  1:2324752/1‑150 (MQ=255)
                                                                                                                   ttaGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCa                              >  1:1216578/1‑150 (MQ=255)
                                                                                                                    taGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCAt                             <  2:3081009/150‑1 (MQ=255)
                                                                                                                    taGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCAt                             >  2:3406339/1‑150 (MQ=255)
                                                                                                                    taGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCAt                             <  2:215104/150‑1 (MQ=255)
                                                                                                                     aGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAATGAAGACCTGTTTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCAtt                            <  1:3957237/150‑1 (MQ=255)
                                                                                                                         atatTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTcgcg                        >  1:193155/1‑150 (MQ=255)
                                                                                                                             tGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTac                    >  1:103555/1‑150 (MQ=255)
                                                                                                                             tGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTac                    >  2:4665988/1‑150 (MQ=255)
                                                                                                                             tGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTac                    >  1:3934756/1‑150 (MQ=255)
                                                                                                                             tGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTac                    >  1:177138/1‑150 (MQ=255)
                                                                                                                                    gggACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGc             >  2:2720942/1‑150 (MQ=255)
                                                                                                                                    gggACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGc             >  2:1240128/1‑150 (MQ=255)
                                                                                                                                        cATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTGTTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACa         >  2:2632927/1‑150 (MQ=255)
                                                                                                                                           ttGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTg      <  2:652645/150‑1 (MQ=255)
                                                                                                                                            tGTCGAACCACGCCAGCAAGGCGACAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGa     >  2:1118735/1‑150 (MQ=255)
                                                                                                                                               cGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAgg  >  1:438522/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CAATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGG  >  CP009273/1770135‑1770427

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: