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,2810TA99.2% 834.4 / ‑2.7 251D146E (GAT→GAAydiFputative acetyl‑CoA:acetoacetyl‑CoA transferase: alpha subunit/beta subunit
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (132/117);  total (133/118)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 251 total aligned reads displayed.

AATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTG  >  CP009273/1770136‑1770423
                                                                                                                                                 |                                                                                                                                              
aaTCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAAcc                                                                                                                                              <  1:1089437/148‑1 (MQ=255)
aaTCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAAcc                                                                                                                                              >  1:1951875/1‑148 (MQ=255)
  tCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCAc                                                                                                                                            >  2:2911014/1‑148 (MQ=255)
    gccgcGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGc                                                                                                                                          <  2:41841/148‑1 (MQ=255)
          tATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCaa                                                                                                                                    >  2:1653606/1‑148 (MQ=255)
           aTTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAg                                                                                                                                   >  1:1148303/1‑148 (MQ=255)
            tttCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAgg                                                                                                                                  >  1:193356/1‑148 (MQ=255)
             ttCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAggc                                                                                                                                 >  1:807416/1‑148 (MQ=255)
                   aCTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCaaa                                                                                                                           >  2:2074348/1‑148 (MQ=255)
                   aCTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCaaa                                                                                                                           <  2:3221224/148‑1 (MQ=255)
                   aCTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCaaa                                                                                                                           <  1:3360854/148‑1 (MQ=255)
                      cGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACtg                                                                                                                        >  1:758272/1‑148 (MQ=255)
                        cAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACtgaa                                                                                                                      <  2:1940031/148‑1 (MQ=255)
                           aaCAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACtgaatga                                                                                                                   <  2:469566/148‑1 (MQ=255)
                           aaCAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACtgaatga                                                                                                                   >  1:1480117/1‑148 (MQ=255)
                           aaCAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACtgaatga                                                                                                                   >  1:780533/1‑148 (MQ=255)
                           aaCAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACtgaatga                                                                                                                   <  2:3005858/148‑1 (MQ=255)
                                aataaaatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCa                                                                                                              <  1:2819314/148‑1 (MQ=255)
                                aataaaatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCa                                                                                                              <  1:3095501/148‑1 (MQ=255)
                                   aaaatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTa                                                                                                           >  2:259605/1‑148 (MQ=255)
                                    aaatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTaa                                                                                                          >  2:974872/1‑148 (MQ=255)
                                    aaatTAGTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTaa                                                                                                          <  2:2274186/148‑1 (MQ=255)
                                     aatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTaaa                                                                                                         <  2:793793/148‑1 (MQ=255)
                                     aatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTaaa                                                                                                         <  1:1125106/148‑1 (MQ=255)
                                      attattGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAaag                                                                                                        >  2:1338957/1‑148 (MQ=255)
                                       ttattGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAaaga                                                                                                       >  1:811216/1‑148 (MQ=255)
                                           tGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCTCGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGAc                                                                                                   >  2:1162277/1‑148 (MQ=255)
                                           tGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGAc                                                                                                   >  2:560709/1‑148 (MQ=255)
                                            gCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGAcc                                                                                                  <  1:2906782/148‑1 (MQ=255)
                                               tataACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGa                                                                                               <  1:3423372/148‑1 (MQ=255)
                                               tataACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGa                                                                                               >  1:3408920/1‑148 (MQ=255)
                                                  aaCTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTa                                                                                            <  2:1959389/148‑1 (MQ=255)
                                                   aCTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTaa                                                                                           <  1:2117845/148‑1 (MQ=255)
                                                   aCTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTaa                                                                                           >  2:772354/1‑148 (MQ=255)
                                                    cTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTaaa                                                                                          <  2:3227698/148‑1 (MQ=255)
                                                    cTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTaaa                                                                                          >  1:1085632/1‑148 (MQ=255)
                                                     tACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAAc                                                                                         >  1:2839595/1‑148 (MQ=255)
                                                       cccACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTg                                                                                       >  2:2399486/1‑148 (MQ=255)
                                                         cacaAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGt                                                                                     >  2:301254/1‑148 (MQ=255)
                                                          acaAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTc                                                                                    >  1:3382746/1‑148 (MQ=255)
                                                           caAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCg                                                                                   <  2:2457050/148‑1 (MQ=255)
                                                             aGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACTTGATTAAACTGGTCGAg                                                                                 <  1:18180/148‑1 (MQ=255)
                                                             aGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAg                                                                                 <  1:3187580/148‑1 (MQ=255)
                                                             aGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAg                                                                                 <  1:19485/148‑1 (MQ=255)
                                                              ggTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGt                                                                                >  2:1661788/1‑148 (MQ=255)
                                                                          caAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAg                                                                    >  1:496012/1‑148 (MQ=255)
                                                                           aaaCCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGa                                                                   <  2:1042221/148‑1 (MQ=255)
                                                                           aaaCCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGa                                                                   <  2:326411/148‑1 (MQ=255)
                                                                           aaaCCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGa                                                                   <  2:2796701/148‑1 (MQ=255)
                                                                             aCCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAat                                                                 <  2:1729630/148‑1 (MQ=255)
                                                                             aCCTTACGCGCCGCCGCAGCCCACCAGCATTGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTTAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAat                                                                 <  1:3255047/148‑1 (MQ=255)
                                                                              ccTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGAATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAata                                                                <  2:895450/148‑1 (MQ=255)
                                                                                 tACGCGCCGCCGCATCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATAtct                                                             <  2:491537/148‑1 (MQ=255)
                                                                                 tACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATAtct                                                             <  2:1846649/148‑1 (MQ=255)
                                                                                    gcgcCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTa                                                          <  1:1879479/148‑1 (MQ=255)
                                                                                    gcgcCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTa                                                          <  1:68312/148‑1 (MQ=255)
                                                                                    gcgcCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTa                                                          <  1:1172226/148‑1 (MQ=255)
                                                                                    gcgcCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTa                                                          <  2:1553988/148‑1 (MQ=255)
                                                                                     cgccgccgcAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTAt                                                         >  1:2487286/1‑148 (MQ=255)
                                                                                     cgccgccgcAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTAt                                                         <  2:1783352/148‑1 (MQ=255)
                                                                                     cgccgccgcAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTAt                                                         >  1:10034/1‑148 (MQ=255)
                                                                                      gccgccgcAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTAtt                                                        <  1:459982/148‑1 (MQ=255)
                                                                                         gccgcAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCTAGTTTGATAACAAAGAATATCTCTATTACa                                                     <  1:2154717/148‑1 (MQ=255)
                                                                                         gccgcAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACa                                                     <  2:1269251/148‑1 (MQ=255)
                                                                                               gCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGa                                               <  2:340514/148‑1 (MQ=255)
                                                                                               gCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGa                                               >  1:2989121/1‑148 (MQ=255)
                                                                                                cccACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGAt                                              >  1:3309099/1‑148 (MQ=255)
                                                                                                cccACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGAt                                              >  2:2892077/1‑148 (MQ=255)
                                                                                                 ccaccaGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGAtt                                             <  2:3452541/148‑1 (MQ=255)
                                                                                                  caccaGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTg                                            >  1:2927823/1‑148 (MQ=255)
                                                                                                   accaGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTgc                                           >  1:3061973/1‑148 (MQ=255)
                                                                                                    ccaGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGa                                               <  1:1524058/143‑1 (MQ=255)
                                                                                                     caGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTgcgc                                         <  1:1958021/148‑1 (MQ=255)
                                                                                                       gCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCa                                       >  2:954197/1‑148 (MQ=255)
                                                                                                       gCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCa                                       >  2:2856318/1‑148 (MQ=255)
                                                                                                        ccTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAg                                      <  2:2384171/148‑1 (MQ=255)
                                                                                                          tGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGat                                    <  1:2316969/148‑1 (MQ=255)
                                                                                                          tGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGat                                    >  1:2779332/1‑148 (MQ=255)
                                                                                                          tGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGat                                    <  1:2673426/148‑1 (MQ=255)
                                                                                                           ggTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGata                                   <  1:2085608/148‑1 (MQ=255)
                                                                                                              attattAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTg                                >  1:1717369/1‑148 (MQ=255)
                                                                                                                tattaGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGcc                              <  2:723210/148‑1 (MQ=255)
                                                                                                                  ttaGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCtt                            <  2:1951875/148‑1 (MQ=255)
                                                                                                                   taGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTc                           >  1:2375962/1‑148 (MQ=255)
                                                                                                                    aGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCa                          >  1:2951612/1‑148 (MQ=255)
                                                                                                                    aGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCa                          <  2:2374528/148‑1 (MQ=255)
                                                                                                                    aGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCa                          <  2:1513138/148‑1 (MQ=255)
                                                                                                                        atatTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTcg                      >  1:3469813/1‑148 (MQ=255)
                                                                                                                        atatTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTcg                      <  2:2540453/148‑1 (MQ=255)
                                                                                                                                 tCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTaccac             >  2:3308394/1‑148 (MQ=255)
                                                                                                                                  cGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTaccacc            >  2:2833292/1‑148 (MQ=255)
                                                                                                                                  cGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTaccacc            >  2:1649851/1‑148 (MQ=255)
                                                                                                                                    ggACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTg          >  2:224559/1‑148 (MQ=255)
                                                                                                                                    ggACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTg          >  1:1370345/1‑148 (MQ=255)
                                                                                                                                    ggACATGTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTg          <  2:2899313/148‑1 (MQ=255)
                                                                                                                                        aTTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGAc      >  1:949979/1‑148 (MQ=255)
                                                                                                                                          ttGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAg    >  1:2867483/1‑148 (MQ=255)
                                                                                                                                          ttGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAg    >  1:1419623/1‑148 (MQ=255)
                                                                                                                                            gTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTg  <  2:715858/148‑1 (MQ=255)
                                                                                                                                            gTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTg  <  1:974872/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
AATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTG  >  CP009273/1770136‑1770423

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: