Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,736,722 C→A 28.5% G210G (GGC→GGA yiaK → 2,3‑diketo‑L‑gulonate reductase, NADH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,736,7220CA28.5% 591.9 / 193.4 285G210G (GGC→GGAyiaK2,3‑diketo‑L‑gulonate reductase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base C (98/105);  new base A (40/41);  total (138/147)
Fisher's exact test for biased strand distribution p-value = 8.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 286 total aligned reads displayed.

GCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTT  >  CP009273/3736576‑3736870
                                                                                                                                                  |                                                                                                                                                    
gCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGt                                                                                                                                                    <  2:3236291/149‑1 (MQ=255)
gCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGt                                                                                                                                                    <  2:370704/149‑1 (MQ=255)
 ccATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGtt                                                                                                                                                   >  1:2347840/1‑149 (MQ=255)
   aTTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTAt                                                                                                                                                 >  1:1729365/1‑149 (MQ=255)
    ttccttccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATc                                                                                                                                                <  2:2730431/149‑1 (MQ=255)
       cttccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCgag                                                                                                                                             <  2:4672613/149‑1 (MQ=255)
        ttccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCgaga                                                                                                                                            >  1:790753/1‑149 (MQ=255)
          ccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGagaag                                                                                                                                          <  1:2969499/149‑1 (MQ=255)
                  aTCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGTTGATGAGGGCAATTTGACCAAAGAACCAGGCGTTATCGAGAAGAATcgccg                                                                                                                                  >  1:3216045/1‑149 (MQ=255)
                  aTCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATcgccg                                                                                                                                  >  1:2089179/1‑149 (MQ=255)
                   tCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATcgccgc                                                                                                                                 >  2:1090457/1‑149 (MQ=255)
                   tCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATcgccgc                                                                                                                                 <  1:116795/149‑1 (MQ=255)
                    cACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCa                                                                                                                                <  1:2214644/149‑1 (MQ=255)
                     aCCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCAt                                                                                                                               >  2:2548568/1‑149 (MQ=255)
                      ccATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGGTTTGATGATGAGGGCAATTTTAACAAAGAAACTGGCGTTATAGAGAAGAATAGCCGCAtt                                                                                                                              >  1:3227890/1‑149 (MQ=255)
                      ccATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCAtt                                                                                                                              >  2:1388035/1‑149 (MQ=255)
                       cATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCAttt                                                                                                                             >  2:1875591/1‑149 (MQ=255)
                       cATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCAttt                                                                                                                             <  2:806281/149‑1 (MQ=255)
                       cATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCAttt                                                                                                                             <  1:1784791/149‑1 (MQ=255)
                         tGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTg                                                                                                                           >  2:1560479/1‑149 (MQ=255)
                         tGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTg                                                                                                                           >  2:160400/1‑149 (MQ=255)
                             cGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGa                                                                                                                       <  1:3686944/149‑1 (MQ=255)
                                tatGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATgg                                                                                                                    >  2:1821918/1‑149 (MQ=255)
                                 atgtcgatgtcgATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATggg                                                                                                                   <  1:1494222/149‑1 (MQ=255)
                                   gtcgatgtcgatGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCt                                                                                                                 <  2:2769120/149‑1 (MQ=255)
                                       atgtcgatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAATAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTg                                                                                                             >  1:2428908/1‑149 (MQ=255)
                                        tgtcgatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACTAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTgg                                                                                                            <  1:4415542/149‑1 (MQ=255)
                                             atgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAgg                                                                                                       >  2:3220650/1‑149 (MQ=255)
                                                ttCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTc                                                                                                    <  1:1939181/149‑1 (MQ=255)
                                                 tctcTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCt                                                                                                   <  2:4176874/149‑1 (MQ=255)
                                                 tctcTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCt                                                                                                   <  2:722777/149‑1 (MQ=255)
                                                 tctcTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCt                                                                                                   <  2:61657/149‑1 (MQ=255)
                                                  ctctTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTg                                                                                                  >  2:3817452/1‑149 (MQ=255)
                                                        cGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGt                                                                                            <  2:3296546/149‑1 (MQ=255)
                                                            aTGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGAt                                                                                        <  1:1061235/149‑1 (MQ=255)
                                                             tGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGAtt                                                                                       <  2:2220469/149‑1 (MQ=255)
                                                             tGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGAtt                                                                                       >  2:2754849/1‑149 (MQ=255)
                                                               ttAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATtgt                                                                                     <  1:3157218/149‑1 (MQ=255)
                                                                tAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATtgtg                                                                                    >  2:207289/1‑149 (MQ=255)
                                                                tAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATtgtg                                                                                    >  1:1177542/1‑149 (MQ=255)
                                                                   aaGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTg                                                                                 >  1:1601575/1‑149 (MQ=255)
                                                                     gTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGa                                                                               >  2:567298/1‑149 (MQ=255)
                                                                      ttAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGGTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGat                                                                              >  1:849095/1‑149 (MQ=255)
                                                                       tAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGata                                                                             <  1:667592/149‑1 (MQ=255)
                                                                       tAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGata                                                                             <  1:3059037/149‑1 (MQ=255)
                                                                        aaCCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGatat                                                                            >  2:1437278/1‑149 (MQ=255)
                                                                        aaCCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGatat                                                                            >  1:2330132/1‑149 (MQ=255)
                                                                         aCCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATg                                                                           >  1:2407538/1‑149 (MQ=255)
                                                                         aCCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATg                                                                           <  2:1419911/149‑1 (MQ=255)
                                                                         aCCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATAGAGAAGAATCGCCGCATTTTGCCGATGGGCTACAGGAAAGGTTATGGCATGTCGATTGTGCTGGAGATg                                                                           >  1:14054/1‑149 (MQ=255)
                                                                          ccGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGa                                                                          <  2:1521230/149‑1 (MQ=255)
                                                                          ccGTCTGGCAGGTCGTCAGCTCCCGGTCGAGGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGa                                                                          <  2:3762467/149‑1 (MQ=255)
                                                                           cGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGAt                                                                         <  2:3523680/149‑1 (MQ=255)
                                                                              cTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGc                                                                      >  2:4451885/1‑149 (MQ=255)
                                                                               tGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCt                                                                     <  2:2793584/149‑1 (MQ=255)
                                                                               tGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCt                                                                     >  1:342484/1‑149 (MQ=255)
                                                                                        gtcAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCtt                                                            >  2:4407201/1‑149 (MQ=255)
                                                                                           aGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTcc                                                         <  1:3407742/149‑1 (MQ=255)
                                                                                           aGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTcc                                                         <  2:2024856/149‑1 (MQ=255)
                                                                                            gCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCg                                                        >  1:801985/1‑149 (MQ=255)
                                                                                            gCTCCCGGTCGATGGTGGCTTTGATGAAGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCTCCGCATTTTGCCGATGGGCTACTGGAAATGTTCTGGCATGTCGATTGTGCTTGATATGATCGCTACTCTCATTTACg                                                        >  2:4557487/1‑149 (MQ=255)
                                                                                               cccGGTCGATGGTGGCTTTTATGATGAGGGCAATTTGACCAAATAACCTGGAGTTATCGAGAAGAATAGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTacgacg                                                     >  1:3324893/1‑149 (MQ=255)
                                                                                                ccGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACgg                                                    <  1:2921158/149‑1 (MQ=255)
                                                                                                 cGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGgc                                                   >  1:3850770/1‑149 (MQ=255)
                                                                                                 cGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGgc                                                   >  1:1461545/1‑149 (MQ=255)
                                                                                                  ggTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGgcg                                                  >  1:1004023/1‑149 (MQ=255)
                                                                                                  ggTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGgcg                                                  >  2:1704645/1‑149 (MQ=255)
                                                                                                  ggTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGgcg                                                  >  1:1003817/1‑149 (MQ=255)
                                                                                                    tCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCa                                                <  1:4473592/149‑1 (MQ=255)
                                                                                                      gATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATc                                              >  1:1885091/1‑149 (MQ=255)
                                                                                                        tggtggCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCg                                            >  2:1188747/1‑149 (MQ=255)
                                                                                                            ggCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGc                                        >  2:295235/1‑149 (MQ=255)
                                                                                                             gCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGcc                                       <  2:4025328/149‑1 (MQ=255)
                                                                                                               tttGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGa                                     >  1:1285129/1‑149 (MQ=255)
                                                                                                                ttGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGaa                                    <  1:4314523/149‑1 (MQ=255)
                                                                                                                     gatgaGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCAc                               <  2:3605134/149‑1 (MQ=255)
                                                                                                                       tgaGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCAccc                             >  2:344922/1‑149 (MQ=255)
                                                                                                                        gaGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCa                            <  1:1560479/149‑1 (MQ=255)
                                                                                                                         aGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAg                           >  2:4709555/1‑149 (MQ=255)
                                                                                                                         aGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATAGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAg                           <  1:2744539/149‑1 (MQ=255)
                                                                                                                          gggCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGACGTCACCCAgg                          >  1:3620470/1‑149 (MQ=255)
                                                                                                                          gggCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAgg                          >  1:404360/1‑149 (MQ=255)
                                                                                                                          gggCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAgg                          >  2:3922855/1‑149 (MQ=255)
                                                                                                                                tttGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGacaaca                    <  2:909706/149‑1 (MQ=255)
                                                                                                                                  tGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGc                  <  2:1907731/149‑1 (MQ=255)
                                                                                                                                    aCCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGcga                <  1:4371956/149‑1 (MQ=255)
                                                                                                                                    aCCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGcga                <  1:4736850/149‑1 (MQ=255)
                                                                                                                                    aCCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGcga                >  1:1374155/1‑149 (MQ=255)
                                                                                                                                    aCCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGcga                >  1:3763762/1‑149 (MQ=255)
                                                                                                                                    aCCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGcga                >  2:3496094/1‑149 (MQ=255)
                                                                                                                                       aaaGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGcgacga             <  1:3745380/149‑1 (MQ=255)
                                                                                                                                       aaaGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGcgacga             <  1:4549900/149‑1 (MQ=255)
                                                                                                                                       aaaGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGcgacga             >  1:2648783/1‑149 (MQ=255)
                                                                                                                                           aaCCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATAc         <  2:17576/149‑1 (MQ=255)
                                                                                                                                              cTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGc      <  2:4120697/149‑1 (MQ=255)
                                                                                                                                              cTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGc      <  1:1852666/149‑1 (MQ=255)
                                                                                                                                               tGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCa     >  1:338770/1‑149 (MQ=255)
                                                                                                                                                ggCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCAt    <  1:3969948/149‑1 (MQ=255)
                                                                                                                                                 gCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCAtt   >  2:2861539/1‑149 (MQ=255)
                                                                                                                                                 gCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCAtt   >  1:4760677/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTT  >  CP009273/3736576‑3736870

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: