Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,513,300 A→T 7.5% Q6L (CAG→CTG)  curA → curcumin/dihydrocurcumin reductase, NADPH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,513,3000AT 7.5% 651.4 / 35.4 214Q6L (CAG→CTG) curAcurcumin/dihydrocurcumin reductase, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base A (102/96);  new base T (8/8);  total (110/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 219 total aligned reads displayed.

ACAGCCCGAAGCAACGGATTCCTCGATGAAATTATTGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATAT  >  CP009273/1513153‑1513444
                                                                                                                                                   |                                                                                                                                                
aCAGCCCGAAGCAACGGATTCCTCGATGAAATTATTGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAagca                                                                                                                                                  >  2:1832540/1‑148 (MQ=255)
aCAGCCCGAAGCAACGGATTCCTCGATGAAATTATTGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAagca                                                                                                                                                  <  2:3038457/148‑1 (MQ=255)
    cccGAAGCAACGGATTCCTCGATGAAATTATTGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGt                                                                                                                                              >  1:471487/1‑148 (MQ=255)
    cccGAAGCAACGGATTCCTCGATGAAATTATTGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGt                                                                                                                                              >  2:609358/1‑148 (MQ=255)
    cccGAAGCAACGGATTCCTCGATGAAATTATTGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGt                                                                                                                                              >  2:185132/1‑148 (MQ=255)
      cGAAGCAACGGATTCCTCGATGAAATTATTGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTaa                                                                                                                                            >  1:310602/1‑148 (MQ=255)
      cGAAGCAACGGAGTCCTCGATGAAATTATTGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTaa                                                                                                                                            <  2:3119434/148‑1 (MQ=255)
       gAAGCAACGGATTCCTCGATGAAATTATTGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAAt                                                                                                                                           <  1:3214379/148‑1 (MQ=255)
         aGCAACGGATTCCTCGATGAAATTATTGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCg                                                                                                                                         <  2:916621/148‑1 (MQ=255)
           cAACGGATTCCTCGATGAAATTATTGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCTGCGTAATCGAc                                                                                                                                       >  1:1373711/1‑148 (MQ=255)
              cGGATTCCTCGATGAAATTATTGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGtt                                                                                                                                    >  1:503207/1‑148 (MQ=255)
                 aTTCCTCGATGAAATTATTGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTggg                                                                                                                                 >  2:1173731/1‑148 (MQ=255)
                 aTTCCTCGATGAAATTATTGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTggg                                                                                                                                 <  2:2662000/148‑1 (MQ=255)
                  ttCCTCGATGAAATTATTGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGt                                                                                                                                >  1:2652842/1‑148 (MQ=255)
                  ttCCTCGATGAAATTATTGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGt                                                                                                                                <  1:1445535/148‑1 (MQ=255)
                  ttCCTCGATGAAATTATTGCGTGTTGCAGAGGGGATGTGACGGCGGCAAACAAAACCAAGCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGt                                                                                                                                <  1:1956626/148‑1 (MQ=255)
                     cTCGATGAAATTATTGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCt                                                                                                                             <  1:3188748/148‑1 (MQ=255)
                            aaaTTATTGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTcg                                                                                                                      >  2:607855/1‑148 (MQ=255)
                            aaaTTATTGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTcg                                                                                                                      >  1:309875/1‑148 (MQ=255)
                               ttattGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCTGCGTAATCGACGTTGGGTTCTGGCCTCGCGt                                                                                                                   <  2:2897043/148‑1 (MQ=255)
                                  ttGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCca                                                                                                                >  2:377114/1‑148 (MQ=255)
                                   tGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCcac                                                                                                               >  2:3097842/1‑148 (MQ=255)
                                    gCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCcaca                                                                                                              >  2:541311/1‑148 (MQ=255)
                                    gCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCcaca                                                                                                              <  1:2970275/148‑1 (MQ=255)
                                    gCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCcaca                                                                                                              <  2:2708379/148‑1 (MQ=255)
                                     cGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACAt                                                                                                             <  1:2438830/148‑1 (MQ=255)
                                       tgtTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATgg                                                                                                           >  2:833520/1‑148 (MQ=255)
                                              gagGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCt                                                                                                    >  2:1913137/1‑148 (MQ=255)
                                                   gATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTcc                                                                                               >  1:876243/1‑148 (MQ=255)
                                                   gATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTcc                                                                                               >  1:3225503/1‑148 (MQ=255)
                                                     tgtgACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCgg                                                                                             <  1:199976/148‑1 (MQ=255)
                                                     tgtgACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCgg                                                                                             <  1:1661961/148‑1 (MQ=255)
                                                       tgACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGgag                                                                                           <  1:1575662/148‑1 (MQ=255)
                                                        gACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGgaga                                                                                          <  2:2404915/148‑1 (MQ=255)
                                                          cGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAAt                                                                                        <  2:1053013/148‑1 (MQ=255)
                                                            gCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAAttt                                                                                      >  2:2228083/1‑148 (MQ=255)
                                                                caaacaaaACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGt                                                                                  <  2:1708636/148‑1 (MQ=255)
                                                                 aaacaaaACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTc                                                                                 <  1:1098801/148‑1 (MQ=255)
                                                                    caaaACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTg                                                                              <  2:1660968/148‑1 (MQ=255)
                                                                    caaaACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTg                                                                              >  1:1265689/1‑148 (MQ=255)
                                                                       aaCCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTgaag                                                                           <  2:1368680/148‑1 (MQ=255)
                                                                        aCCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTgaaga                                                                          <  1:389032/148‑1 (MQ=255)
                                                                          cAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGaagaag                                                                        >  2:2370312/1‑148 (MQ=255)
                                                                           aaTCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAagaaga                                                                       >  1:91385/1‑148 (MQ=255)
                                                                             tCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAgatg                                                                     >  1:1470752/1‑148 (MQ=255)
                                                                             tCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAgatg                                                                     <  2:2055881/148‑1 (MQ=255)
                                                                             tCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAgatg                                                                     >  1:562262/1‑148 (MQ=255)
                                                                                cGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGatgatg                                                                  <  1:63886/148‑1 (MQ=255)
                                                                                           tactaTTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCTGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCg                                                       <  2:1373711/148‑1 (MQ=255)
                                                                                           tactaTTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCg                                                       <  1:2595763/148‑1 (MQ=255)
                                                                                              taTTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGt                                                    >  2:1464315/1‑148 (MQ=255)
                                                                                              taTTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGt                                                    >  2:3095182/1‑148 (MQ=255)
                                                                                               aTTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTg                                                   >  2:329115/1‑148 (MQ=255)
                                                                                                ttGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGa                                                  >  1:1613639/1‑148 (MQ=255)
                                                                                                ttGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGa                                                  >  1:1944995/1‑148 (MQ=255)
                                                                                                 tGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTACCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGaa                                                 >  2:513621/1‑148 (MQ=255)
                                                                                                 tGATGGATTCATGACAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGAAATCTACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCTCCACACAGGGTGaa                                                 >  2:334912/1‑148 (MQ=255)
                                                                                                   aTGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAgg                                               <  1:2854261/148‑1 (MQ=255)
                                                                                                    tGGATTCATGCCAACAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGa                                              <  2:725967/148‑1 (MQ=255)
                                                                                                    tGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGa                                              <  2:1214643/148‑1 (MQ=255)
                                                                                                     ggATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGAc                                             <  1:1195510/148‑1 (MQ=255)
                                                                                                     ggATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGAc                                             <  1:1854502/148‑1 (MQ=255)
                                                                                                        ttCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAgg                                          <  1:2370312/148‑1 (MQ=255)
                                                                                                        ttCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAgg                                          <  1:2168132/148‑1 (MQ=255)
                                                                                                         tCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGgt                                         <  1:2571347/148‑1 (MQ=255)
                                                                                                           aTGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGgtgt                                       <  1:278211/148‑1 (MQ=255)
                                                                                                             gCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTa                                     >  2:1493803/1‑148 (MQ=255)
                                                                                                             gCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTa                                     <  1:833520/148‑1 (MQ=255)
                                                                                                              ccAAAAGCACAGGAGAGCATTATGGGGCAACAACAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTAc                                    <  1:344536/148‑1 (MQ=255)
                                                                                                              ccAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTAc                                    <  2:2523167/148‑1 (MQ=255)
                                                                                                                aaaaGCACAGGAGAGCATTATGGGGCAACAAAAGCTGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTg                                  <  2:719694/148‑1 (MQ=255)
                                                                                                                 aaaGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTgc                                 >  2:1095675/1‑148 (MQ=255)
                                                                                                                 aaaGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTgc                                 >  1:2958757/1‑148 (MQ=255)
                                                                                                                  aaGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTgcg                                >  2:280475/1‑148 (MQ=255)
                                                                                                                  aaGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTgcg                                >  2:567323/1‑148 (MQ=255)
                                                                                                                   aGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTgcgc                               >  1:1002505/1‑148 (MQ=255)
                                                                                                                     cacaGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGcac                             >  1:576287/1‑148 (MQ=255)
                                                                                                                     cacaGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGcac                             <  1:1510218/148‑1 (MQ=255)
                                                                                                                     cacaGGAGAGCATGATGGGGCAACAACAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGcac                             <  2:3199744/148‑1 (MQ=255)
                                                                                                                      acaGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGcaca                            <  1:407176/148‑1 (MQ=255)
                                                                                                                      acaGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGcaca                            <  1:513642/148‑1 (MQ=255)
                                                                                                                       caGGAGAGCATTATGGGGCAACAAAAGCTGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAg                           >  1:754403/1‑148 (MQ=255)
                                                                                                                        agtatatcATTATGGTGCAACAAAAGCAGCGTAATCGACGTTGTGTTCTGGCCTCGAGTACACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACAACGTGTGAAGGACAGGTGTTACTGCGCACAGt                          <  2:3269292/141‑1 (MQ=255)
                                                                                                                        aGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGt                          <  2:1142495/148‑1 (MQ=255)
                                                                                                                        aGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGt                          <  2:1009005/148‑1 (MQ=255)
                                                                                                                        aGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGt                          <  1:3083624/148‑1 (MQ=255)
                                                                                                                         ggAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGtt                         <  1:376185/148‑1 (MQ=255)
                                                                                                                           agagCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTa                       <  2:878269/148‑1 (MQ=255)
                                                                                                                            gagCATTATGGGGCAACAAAAGCTGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTAt                      <  1:1712359/148‑1 (MQ=255)
                                                                                                                             agCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTAtt                     >  1:925456/1‑148 (MQ=255)
                                                                                                                                  tATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTcc                >  2:1443203/1‑148 (MQ=255)
                                                                                                                                  tATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTcc                >  2:1443205/1‑148 (MQ=255)
                                                                                                                                    tGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCt              <  2:310602/148‑1 (MQ=255)
                                                                                                                                     ggggCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTg             <  1:1634022/148‑1 (MQ=255)
                                                                                                                                         caacaaAAGCTGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGAcc         >  1:1335406/1‑148 (MQ=255)
                                                                                                                                          aacaaAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGAccc        <  2:373503/148‑1 (MQ=255)
                                                                                                                                          aacaaAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGAccc        <  1:1844056/148‑1 (MQ=255)
                                                                                                                                              aaaGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCTTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGtat    >  1:2780788/1‑148 (MQ=255)
                                                                                                                                               aaGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACAAGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGGGTATTTGTCCCTGGACCCGtata   >  2:2984471/1‑148 (MQ=255)
                                                                                                                                                agcagcGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTatat  >  1:683826/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ACAGCCCGAAGCAACGGATTCCTCGATGAAATTATTGCGTGTTGCAGAGGGGATGTGACGGCTGCAAACAAAACCAATCCCGTGAGCCTGCTACTATTGATGGATTCATGCCAAAAGCACAGGAGAGCATTATGGGGCAACAAAAGCAGCGTAATCGACGTTGGGTTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATAT  >  CP009273/1513153‑1513444

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: