Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 450,766 G→T 100% A60S (GCG→TCG)  tig → peptidyl‑prolyl cis/trans isomerase (trigger factor)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273450,7660GT95.0% 589.4 / 20.5 199A60S (GCG→TCG) tigpeptidyl‑prolyl cis/trans isomerase (trigger factor)
Reads supporting (aligned to +/- strand):  ref base G (9/1);  new base T (77/112);  total (86/113)
Fisher's exact test for biased strand distribution p-value = 2.60e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 200 total aligned reads displayed.

TGGCCGCCGTGTAACGATTACTATCGCTGCTGACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCGCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGACTTCACTTACTCTGTAG  >  CP009273/450621‑450910
                                                                                                                                                 |                                                                                                                                                
   ccgccgTGTAACGATTACTATCGCTGCTGACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCt                                                                                                                                             <  2:2683213/148‑1 (MQ=255)
   ccgccgTGTAACGATTACTATCGCTGCTGACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCt                                                                                                                                             <  1:468049/148‑1 (MQ=255)
    cgccgTGTAACGATTACTATCGCTGCTGACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTg                                                                                                                                            >  1:3062771/1‑148 (MQ=255)
         tgtAACGATTACTATCGCTGCTGACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGc                                                                                                                                       <  1:561974/148‑1 (MQ=255)
         tgtAACGATTACTATCGCTGCTGACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGc                                                                                                                                       <  1:2347443/148‑1 (MQ=255)
           tAACGATTACTATCGCTGCTGACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCa                                                                                                                                     <  2:959949/148‑1 (MQ=255)
           tAACGATTACTATCGCTGCTGACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCa                                                                                                                                     <  1:1041652/148‑1 (MQ=255)
              cGATTACTATCGCTGCTGACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGa                                                                                                                                  >  2:2041602/1‑148 (MQ=255)
               gATTACTATCGCTGCTGACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGAc                                                                                                                                 >  1:2022558/1‑148 (MQ=255)
                aTTACTATCGCTGCTGACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACg                                                                                                                                >  2:346533/1‑148 (MQ=255)
                 ttACTATCGCTGCTGACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGt                                                                                                                               >  2:2912052/1‑148 (MQ=255)
                  tACTATCGCTGCTGACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGtt                                                                                                                              >  2:2553326/1‑148 (MQ=255)
                  tACTATCGCTGCTGACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCGCGTCTGTACGCCAGGACGtt                                                                                                                              >  1:2470526/1‑148 (MQ=255)
                         gctgctGACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTg                                                                                                                       >  2:3478238/1‑148 (MQ=255)
                           tgctgACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGAc                                                                                                                     <  1:3287376/148‑1 (MQ=255)
                            gctgACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGAcc                                                                                                                    <  1:360787/148‑1 (MQ=255)
                             ctgACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGCTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCt                                                                                                                   >  1:2259445/1‑148 (MQ=255)
                              tgACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCtg                                                                                                                  <  1:1912755/148‑1 (MQ=255)
                               gACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCtga                                                                                                                 >  2:1637331/1‑148 (MQ=255)
                               gACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCtga                                                                                                                 >  2:1637220/1‑148 (MQ=255)
                                 cAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCGCGTCTGTACGCCAGGACGTTCTGGGTGACCtgatg                                                                                                               >  1:3217668/1‑148 (MQ=255)
                                  aGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCtgatga                                                                                                              <  2:3153628/148‑1 (MQ=255)
                                  aGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCtgatga                                                                                                              <  2:118583/148‑1 (MQ=255)
                                  aGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCtgatga                                                                                                              <  1:1651216/148‑1 (MQ=255)
                                     aTCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGcc                                                                                                           <  2:478070/148‑1 (MQ=255)
                                      tCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCg                                                                                                          >  1:1687158/1‑148 (MQ=255)
                                       cGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGt                                                                                                         >  2:1793793/1‑148 (MQ=255)
                                           aCCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAATTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACt                                                                                                     <  1:767440/148‑1 (MQ=255)
                                            ccGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACtt                                                                                                    >  2:2889724/1‑148 (MQ=255)
                                              gCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCa                                                                                                  <  2:3378224/148‑1 (MQ=255)
                                               cTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCAt                                                                                                 <  1:1969802/148‑1 (MQ=255)
                                                tGTTAATAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCACAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCAtt                                                                                                <  1:3259713/148‑1 (MQ=255)
                                                tGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCAtt                                                                                                <  2:3233031/148‑1 (MQ=255)
                                                 gTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGATCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTg                                                                                               >  2:3359165/1‑148 (MQ=255)
                                                        gCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCcat                                                                                        <  1:2278627/148‑1 (MQ=255)
                                                        gCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGAcat                                                                                        >  1:301864/1‑148 (MQ=255)
                                                           aGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCcatcat                                                                                     >  2:1604379/1‑148 (MQ=255)
                                                                   aaCGTTGCGAAAAAAGTACGTAGTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGaaa                                                                             <  2:2856378/148‑1 (MQ=255)
                                                                     cGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGaaaaa                                                                           >  1:1997953/1‑148 (MQ=255)
                                                                       ttGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAaat                                                                         <  2:3074321/148‑1 (MQ=255)
                                                                       ttGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAaat                                                                         <  2:1304474/148‑1 (MQ=255)
                                                                       ttGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAaat                                                                         >  2:657563/1‑148 (MQ=255)
                                                                        tGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAaatc                                                                        <  2:3431350/148‑1 (MQ=255)
                                                                           gAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAaatcaat                                                                     <  1:1638409/148‑1 (MQ=255)
                                                                           gAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCGCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAaatcaat                                                                     <  1:2801665/148‑1 (MQ=255)
                                                                            aaaaaaGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAaatcaatc                                                                    <  1:3478238/148‑1 (MQ=255)
                                                                            aaaaaaGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAaatcaatc                                                                    <  1:500178/148‑1 (MQ=255)
                                                                              aaaaGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCg                                                                  >  2:1800942/1‑148 (MQ=255)
                                                                               aaaGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCgg                                                                 <  2:854008/148‑1 (MQ=255)
                                                                                aaGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGc                                                                <  1:1011644/148‑1 (MQ=255)
                                                                                 aGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCt                                                               >  2:3553577/1‑148 (MQ=255)
                                                                                   tACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTgg                                                             >  1:2158444/1‑148 (MQ=255)
                                                                                        aTTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCAc                                                        <  1:3046624/148‑1 (MQ=255)
                                                                                          tGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCg                                                      <  1:1174992/148‑1 (MQ=255)
                                                                                           gACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGa                                                     <  2:2054742/148‑1 (MQ=255)
                                                                                           gACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGa                                                     <  1:1494381/148‑1 (MQ=255)
                                                                                           gACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGa                                                     >  1:2361630/1‑148 (MQ=255)
                                                                                            aCGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGAc                                                    <  1:3402427/148‑1 (MQ=255)
                                                                                            aCGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGAc                                                    >  1:2350218/1‑148 (MQ=255)
                                                                                               gCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTa                                                 <  2:748491/148‑1 (MQ=255)
                                                                                                 ttCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCGCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATg                                               >  1:140894/1‑148 (MQ=255)
                                                                                                  tCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGt                                              >  2:2553992/1‑148 (MQ=255)
                                                                                                         aggcaaagTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGc                                       >  1:2329119/1‑148 (MQ=255)
                                                                                                         aggcaaagTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGc                                       >  2:2271792/1‑148 (MQ=255)
                                                                                                          ggcaaagTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCg                                      <  1:892327/148‑1 (MQ=255)
                                                                                                            caaagTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGaa                                    <  1:1815359/148‑1 (MQ=255)
                                                                                                            caaagTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGaa                                    <  2:954480/148‑1 (MQ=255)
                                                                                                             aaagTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAAt                                   <  1:2455101/148‑1 (MQ=255)
                                                                                                              aagTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATa                                  <  1:3494443/148‑1 (MQ=255)
                                                                                                               agTGCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATAc                                 <  1:933379/148‑1 (MQ=255)
                                                                                                               agTGCCAATGAATATCGTTGCTCAGCGTTATGGCGCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATAc                                 >  2:50934/1‑148 (MQ=255)
                                                                                                                  gCCAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAg                              <  1:1637316/148‑1 (MQ=255)
                                                                                                                   ccAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGc                             <  2:2022558/148‑1 (MQ=255)
                                                                                                                    cAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCt                            <  1:461683/148‑1 (MQ=255)
                                                                                                                    cAATGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCt                            <  1:649795/148‑1 (MQ=255)
                                                                                                                     aaTGAATATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTg                           >  2:1469686/1‑148 (MQ=255)
                                                                                                                         aaTATCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTg                       >  1:1745923/1‑148 (MQ=255)
                                                                                                                          atatCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGa                      <  2:1667509/148‑1 (MQ=255)
                                                                                                                          atatCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGa                      >  1:1072660/1‑148 (MQ=255)
                                                                                                                             tCGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGa                   >  2:3500027/1‑148 (MQ=255)
                                                                                                                              cGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGAc                  <  1:2732016/148‑1 (MQ=255)
                                                                                                                              cGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGAc                  >  1:3418883/1‑148 (MQ=255)
                                                                                                                              cGTTGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGAc                  <  1:2317602/148‑1 (MQ=255)
                                                                                                                                ttGCTCTGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACGTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTTGCGCACCGACTTATGTTCCGGGCGAATACCAGCTGGGTGAAGACtt                <  2:2950852/148‑1 (MQ=255)
                                                                                                                                ttGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGACtt                >  1:2063487/1‑148 (MQ=255)
                                                                                                                                 tGCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGACTTc               <  1:208988/148‑1 (MQ=255)
                                                                                                                                  gCTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGACTTCa              >  1:131620/1‑148 (MQ=255)
                                                                                                                                   cTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGACTTCAc             <  1:1582289/148‑1 (MQ=255)
                                                                                                                                   cTCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGACTTCAc             >  2:866747/1‑148 (MQ=255)
                                                                                                                                    tCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGACTTCACt            <  1:2768524/148‑1 (MQ=255)
                                                                                                                                    tCAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGACTTCACt            <  1:1399348/148‑1 (MQ=255)
                                                                                                                                     cAGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGACTTCACtt           <  2:2896772/148‑1 (MQ=255)
                                                                                                                                      aGCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGACTTCACTTa          <  1:3460342/148‑1 (MQ=255)
                                                                                                                                       gCGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGACTTCACTTAc         >  2:232743/1‑148 (MQ=255)
                                                                                                                                        cGTTATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGACTTCACTTAct        >  2:1969083/1‑148 (MQ=255)
                                                                                                                                           tATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGACTTCACTTACTCTg     <  1:1768820/148‑1 (MQ=255)
                                                                                                                                           tATGGCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGACTTCACTTACTCTg     <  2:2262714/148‑1 (MQ=255)
                                                                                                                                              ggCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGACTTCACTTACTCTGTag  >  2:3050482/1‑148 (MQ=255)
                                                                                                                                              ggCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGACTTCACTTACTCTGTag  <  2:2748815/148‑1 (MQ=255)
                                                                                                                                              ggCTCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGACTTCACTTACTCTGTag  <  2:1318037/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TGGCCGCCGTGTAACGATTACTATCGCTGCTGACAGCATCGAGACCGCTGTTAAAAGCGAGCTGGTCAACGTTGCGAAAAAAGTACGTATTGACGGCTTCCGCAAAGGCAAAGTGCCAATGAATATCGTTGCTCAGCGTTATGGCGCGTCTGTACGCCAGGACGTTCTGGGTGACCTGATGAGCCGTAACTTCATTGACGCCATCATTAAAGAAAAAATCAATCCGGCTGGCGCACCGACTTATGTTCCGGGCGAATACAAGCTGGGTGAAGACTTCACTTACTCTGTAG  >  CP009273/450621‑450910

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: