Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,330,955 G→T 43.6% D290Y (GAT→TAT)  acnA → aconitate hydratase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,330,9550GT43.6% 269.1 / inf 290D290Y (GAT→TAT) acnAaconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base G (85/78);  new base T (57/70);  total (142/148)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 295 total aligned reads displayed.

CCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTG  >  CP009273/1330808‑1331104
                                                                                                                                                   |                                                                                                                                                     
ccGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAAGTTTATGGTGAt                                                                                                                                                     >  2:4493150/1‑150 (MQ=255)
 cGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATg                                                                                                                                                    <  2:1391432/150‑1 (MQ=255)
   gTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGt                                                                                                                                                  >  2:4035907/1‑150 (MQ=255)
     ttCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCt                                                                                                                                                >  1:112866/1‑150 (MQ=255)
     ttCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCt                                                                                                                                                <  2:1767942/150‑1 (MQ=255)
      tCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTg                                                                                                                                               <  1:557308/150‑1 (MQ=255)
       ccATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTgg                                                                                                                                              <  2:3962722/150‑1 (MQ=255)
        cATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGa                                                                                                                                             <  2:424546/150‑1 (MQ=255)
           gCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTc                                                                                                                                          >  1:106390/1‑150 (MQ=255)
           gCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTc                                                                                                                                          <  1:2499867/150‑1 (MQ=255)
            cTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCa                                                                                                                                         >  2:3262311/1‑150 (MQ=255)
            cTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCa                                                                                                                                         <  2:4360735/150‑1 (MQ=255)
            cTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCa                                                                                                                                         >  2:595618/1‑150 (MQ=255)
             ttATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCAc                                                                                                                                        <  1:4570165/150‑1 (MQ=255)
             ttATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCAc                                                                                                                                        <  2:1621037/150‑1 (MQ=255)
               aTCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTa                                                                                                                                      >  2:1014858/1‑150 (MQ=255)
                tCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTAc                                                                                                                                     <  2:4416503/150‑1 (MQ=255)
                 cccGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTAcc                                                                                                                                    <  1:4434395/150‑1 (MQ=255)
                  ccGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCg                                                                                                                                   <  1:1229288/150‑1 (MQ=255)
                   cGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGt                                                                                                                                  >  1:1423264/1‑150 (MQ=255)
                   cGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGt                                                                                                                                  >  1:1310249/1‑150 (MQ=255)
                    ggATGTAGTGGGCTTCAAACTGACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGtt                                                                                                                                 <  1:3685504/150‑1 (MQ=255)
                     gATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTg                                                                                                                                <  1:4237132/150‑1 (MQ=255)
                     gATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTg                                                                                                                                <  1:1832081/150‑1 (MQ=255)
                      aTGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTgg                                                                                                                               >  2:3310731/1‑150 (MQ=255)
                          aGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGa                                                                                                                           <  1:380867/150‑1 (MQ=255)
                            tGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATc                                                                                                                         <  1:1570421/150‑1 (MQ=255)
                            tGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATc                                                                                                                         >  1:3162370/1‑150 (MQ=255)
                             gggCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATcg                                                                                                                        <  2:2721179/150‑1 (MQ=255)
                              ggCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATcgc                                                                                                                       <  2:777741/150‑1 (MQ=255)
                                   cAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGccac                                                                                                                  >  1:1791150/1‑150 (MQ=255)
                                   cAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGccac                                                                                                                  >  1:379427/1‑150 (MQ=255)
                                    aaaCTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGccacc                                                                                                                 >  2:681569/1‑150 (MQ=255)
                                     aaCTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGccacca                                                                                                                >  2:1903084/1‑150 (MQ=255)
                                      aCTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCAt                                                                                                               >  2:3041793/1‑150 (MQ=255)
                                       cTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCAtt                                                                                                              <  1:837627/150‑1 (MQ=255)
                                        ttACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTg                                                                                                             <  2:4689451/150‑1 (MQ=255)
                                            cGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCaa                                                                                                         <  2:1078417/150‑1 (MQ=255)
                                              gAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAata                                                                                                       >  1:2104530/1‑150 (MQ=255)
                                               aaaaTTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAatat                                                                                                      >  1:1191307/1‑150 (MQ=255)
                                                aaaTTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATg                                                                                                     <  1:1701873/150‑1 (MQ=255)
                                                aaaTTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATg                                                                                                     >  2:1819344/1‑150 (MQ=255)
                                                 aaTTACGTGAAGGTATTACCGCAACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGTCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGt                                                                                                    <  2:2153692/150‑1 (MQ=255)
                                                  aTTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTc                                                                                                   <  2:445249/150‑1 (MQ=255)
                                                  aTTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTc                                                                                                   <  1:1604566/150‑1 (MQ=255)
                                                      cGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCa                                                                                               <  2:4141327/150‑1 (MQ=255)
                                                       gTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAg                                                                                              <  1:297691/150‑1 (MQ=255)
                                                         gAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGaa                                                                                            <  1:2686537/150‑1 (MQ=255)
                                                          aaGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAat                                                                                           >  1:1456941/1‑150 (MQ=255)
                                                              tATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATgg                                                                                       <  2:2199804/150‑1 (MQ=255)
                                                              tATTACCGCCACAGACCTGGATCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATgg                                                                                       <  1:1956458/150‑1 (MQ=255)
                                                               aTTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGt                                                                                      <  1:934711/150‑1 (MQ=255)
                                                               aTTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGt                                                                                      >  2:425663/1‑150 (MQ=255)
                                                                ttACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGACTTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTg                                                                                     >  2:1506680/1‑150 (MQ=255)
                                                                ttACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTg                                                                                     >  2:905947/1‑150 (MQ=255)
                                                                  aCCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGcc                                                                                   >  2:701289/1‑150 (MQ=255)
                                                                  aCCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGcc                                                                                   >  2:2462566/1‑150 (MQ=255)
                                                                       cacaGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCtg                                                                              >  1:3996502/1‑150 (MQ=255)
                                                                         caGACCTGGTTCTCACTGTTACCCAAAGGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCtgtg                                                                            <  2:2903227/150‑1 (MQ=255)
                                                                             ccTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGctt                                                                        <  1:914510/150‑1 (MQ=255)
                                                                             ccTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGctt                                                                        >  1:4036984/1‑150 (MQ=255)
                                                                               tGGTTCTCACTGTGACTCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGcttct                                                                      <  2:3671678/150‑1 (MQ=255)
                                                                                 gTTCTCACTGATACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATCTTATGGTTATGGTCTGGATTCACTACCGATGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCttcttc                                                                    >  2:3644818/1‑150 (MQ=255)
                                                                                   tctcACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTccc                                                                  >  2:222277/1‑150 (MQ=255)
                                                                                    ctcACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCa                                                                 <  1:1871835/150‑1 (MQ=255)
                                                                                    ctcACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCa                                                                 <  1:2096497/150‑1 (MQ=255)
                                                                                    ctcACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCa                                                                 <  1:3973023/150‑1 (MQ=255)
                                                                                             aCCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCt                                                        <  2:2652786/150‑1 (MQ=255)
                                                                                               ccAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGt                                                      <  1:4531743/150‑1 (MQ=255)
                                                                                                cAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTa                                                     >  2:1999863/1‑150 (MQ=255)
                                                                                                cAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTa                                                     >  2:4254055/1‑150 (MQ=255)
                                                                                                    tgctgcGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAAccc                                                 <  2:1995807/150‑1 (MQ=255)
                                                                                                      ctgcGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTc                                               <  2:1676076/150‑1 (MQ=255)
                                                                                                       tgcGCAAACCTGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCg                                              <  2:834468/150‑1 (MQ=255)
                                                                                                        gcgcAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGa                                             <  2:3443473/150‑1 (MQ=255)
                                                                                                        gcgcAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTACCAAACGATGCTGTAACCCACGa                                             >  1:2580489/1‑150 (MQ=255)
                                                                                                        gcgcAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGa                                             >  1:1770202/1‑150 (MQ=255)
                                                                                                          gcAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGAtt                                           >  1:4059428/1‑150 (MQ=255)
                                                                                                           cAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTa                                          <  1:3364708/150‑1 (MQ=255)
                                                                                                                    cGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGttt                                 <  1:2637198/150‑1 (MQ=255)
                                                                                                                       ggtggGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAg                              <  2:3819957/150‑1 (MQ=255)
                                                                                                                           gggAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCggg                          >  2:1225521/1‑150 (MQ=255)
                                                                                                                           gggAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCggg                          <  1:2828370/150‑1 (MQ=255)
                                                                                                                           gggAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCggg                          >  2:2289915/1‑150 (MQ=255)
                                                                                                                           gggAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACGACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCggg                          <  2:4588527/150‑1 (MQ=255)
                                                                                                                             gAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGgcg                        >  1:655167/1‑150 (MQ=255)
                                                                                                                               aaTTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCa                      >  1:218265/1‑150 (MQ=255)
                                                                                                                               aaTTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCa                      >  1:3516679/1‑150 (MQ=255)
                                                                                                                                    gtcgAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGaa                 >  2:3536929/1‑150 (MQ=255)
                                                                                                                                    gtcgAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGaa                 <  1:3636486/150‑1 (MQ=255)
                                                                                                                                      cgAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGa               >  1:1793124/1‑150 (MQ=255)
                                                                                                                                        aaTTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATc             >  2:1347160/1‑150 (MQ=255)
                                                                                                                                         aTTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCa            >  2:699339/1‑150 (MQ=255)
                                                                                                                                          ttttATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAg           >  1:2256358/1‑150 (MQ=255)
                                                                                                                                          ttttATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAg           >  1:1646669/1‑150 (MQ=255)
                                                                                                                                          tgttatggttatggTCTGGAGTCACTACCGGTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCACTCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAg           <  1:2457508/148‑1 (MQ=255)
                                                                                                                                            ttATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGt         <  1:2748983/150‑1 (MQ=255)
                                                                                                                                            ttATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGt         >  1:2880819/1‑150 (MQ=255)
                                                                                                                                              aTGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCg       <  1:425679/150‑1 (MQ=255)
                                                                                                                                                 gTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGt    >  1:2298668/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTG  >  CP009273/1330808‑1331104

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: