Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,820,888 C→A 88.6% C98F (TGC→TTC)  astE ← succinylglutamate desuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,820,8880CA88.6% 410.4 / 48.5 166C98F (TGC→TTC) astEsuccinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base C (9/10);  new base A (71/76);  total (80/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 170 total aligned reads displayed.

TAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTG  >  CP009273/1820743‑1821035
                                                                                                                                                 |                                                                                                                                                   
tAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGaaa                                                                                                                                                   <  2:1275860/148‑1 (MQ=255)
 aGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAAt                                                                                                                                                  <  2:2539733/148‑1 (MQ=255)
 aGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAAt                                                                                                                                                  <  2:1719655/148‑1 (MQ=255)
  gATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATa                                                                                                                                                 <  1:3006547/148‑1 (MQ=255)
    tCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAAc                                                                                                                                               <  1:1737778/148‑1 (MQ=255)
    tCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAAc                                                                                                                                               >  2:23467/1‑148 (MQ=255)
     cAATGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACg                                                                                                                                              >  1:1078695/1‑148 (MQ=255)
     cAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGANAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACg                                                                                                                                              <  1:1683242/148‑1 (MQ=255)
     cAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACGGTTCCTGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACg                                                                                                                                              <  1:1152266/148‑1 (MQ=255)
      aaGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGt                                                                                                                                             <  1:2309758/148‑1 (MQ=255)
      aaGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGt                                                                                                                                             <  1:951878/148‑1 (MQ=255)
      aaGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGt                                                                                                                                             <  1:1318668/148‑1 (MQ=255)
       aGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGtt                                                                                                                                            >  2:1376755/1‑148 (MQ=255)
       aGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGtt                                                                                                                                            >  1:3153219/1‑148 (MQ=255)
        ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGttt                                                                                                                                           >  1:1754112/1‑148 (MQ=255)
        ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGttt                                                                                                                                           >  1:1754113/1‑148 (MQ=255)
          tGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTcc                                                                                                                                         <  2:1917901/148‑1 (MQ=255)
           gCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTccc                                                                                                                                        <  2:419603/148‑1 (MQ=255)
              aGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTg                                                                                                                                     >  1:1497988/1‑148 (MQ=255)
              aGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTg                                                                                                                                     <  1:2463741/148‑1 (MQ=255)
               gcgcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCACTAACGTTTCCCTTGc                                                                                                                                    >  2:2776639/1‑148 (MQ=255)
               gcgcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGc                                                                                                                                    <  2:1400037/148‑1 (MQ=255)
                 gcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCtt                                                                                                                                  <  2:699622/148‑1 (MQ=255)
                 gcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCtt                                                                                                                                  <  2:699550/148‑1 (MQ=255)
                 gcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCtt                                                                                                                                  <  2:2435173/148‑1 (MQ=255)
                  cacaGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTc                                                                                                                                 <  2:1718137/148‑1 (MQ=255)
                       aTTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAgcgc                                                                                                                            >  1:942940/1‑148 (MQ=255)
                       aTTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAgcgc                                                                                                                            >  1:1656230/1‑148 (MQ=255)
                        ttCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCa                                                                                                                           >  2:2798976/1‑148 (MQ=255)
                        ttCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCa                                                                                                                           <  2:466770/148‑1 (MQ=255)
                        ttCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCa                                                                                                                           >  1:2725453/1‑148 (MQ=255)
                        ttCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCa                                                                                                                           <  2:466557/148‑1 (MQ=255)
                        ttCTTTGCCCGGGTAATAAAAATCTTCCAGGCTCTGTTCCAGTTCGCGCGCCCTACAGGTTTCTCCGCTTTCAGCTATTAGCTGCCAACGACCGCCATACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCa                                                                                                                           <  2:466156/148‑1 (MQ=255)
                                ccTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGAtt                                                                                                                   <  1:632742/148‑1 (MQ=255)
                                ccTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGAtt                                                                                                                   <  2:203124/148‑1 (MQ=255)
                                 cTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTc                                                                                                                  <  1:347723/148‑1 (MQ=255)
                                  tGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTcc                                                                                                                 >  1:223060/1‑148 (MQ=255)
                                      cATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAg                                                                                                             >  2:2319033/1‑148 (MQ=255)
                                       aTAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAgg                                                                                                            <  2:2725068/148‑1 (MQ=255)
                                       aTAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAgg                                                                                                            <  1:1582328/148‑1 (MQ=255)
                                        tAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGa                                                                                                           <  1:1792700/148‑1 (MQ=255)
                                             aTCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCAcc                                                                                                      <  2:1728609/148‑1 (MQ=255)
                                                   cAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAAc                                                                                                >  2:831562/1‑148 (MQ=255)
                                                    aGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAccc                                                                                               >  2:2065414/1‑148 (MQ=255)
                                                    aGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAAcc                                                                                               >  2:1523815/1‑148 (MQ=255)
                                                      gCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGc                                                                                             >  2:3289754/1‑148 (MQ=255)
                                                       cACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGcc                                                                                            >  1:95/1‑148 (MQ=255)
                                                       cACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGcc                                                                                            >  2:1067289/1‑148 (MQ=255)
                                                        aCTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCa                                                                                           <  1:1032427/148‑1 (MQ=255)
                                                        aCTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCa                                                                                           <  2:347448/148‑1 (MQ=255)
                                                        aCTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCa                                                                                           <  1:1032402/148‑1 (MQ=255)
                                                             tCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTa                                                                                      >  1:2592875/1‑148 (MQ=255)
                                                              ccAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTaa                                                                                     <  2:1624246/148‑1 (MQ=255)
                                                                aGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACg                                                                                   >  2:1074697/1‑148 (MQ=255)
                                                                  ttCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACggg                                                                                 >  1:833158/1‑148 (MQ=255)
                                                                  ttCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACggg                                                                                 >  2:1534555/1‑148 (MQ=255)
                                                                  ttCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACggg                                                                                 >  1:861411/1‑148 (MQ=255)
                                                                   tCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGa                                                                                <  1:506097/148‑1 (MQ=255)
                                                                   tCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGa                                                                                >  2:1028710/1‑148 (MQ=255)
                                                                    cgcgcgCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGAt                                                                               >  2:317006/1‑148 (MQ=255)
                                                                     gcgcgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGAtc                                                                              <  1:2458303/148‑1 (MQ=255)
                                                                     gcgcgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGAtc                                                                              <  1:2457947/148‑1 (MQ=255)
                                                                      cgcgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGAtct                                                                             <  1:58214/148‑1 (MQ=255)
                                                                      cgcgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGAtct                                                                             <  2:2707369/148‑1 (MQ=255)
                                                                        cgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCg                                                                           >  2:2203410/1‑148 (MQ=255)
                                                                              aCAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTga                                                                     >  2:2817132/1‑148 (MQ=255)
                                                                               cAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTgag                                                                    <  2:1239030/148‑1 (MQ=255)
                                                                                  gTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTACCGGGATCTCGCCGTGAGata                                                                 >  2:1840127/1‑148 (MQ=255)
                                                                                  gTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGata                                                                 >  2:3276223/1‑148 (MQ=255)
                                                                                         cGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGcc                                                          >  2:1457696/1‑148 (MQ=255)
                                                                                           cTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCaa                                                        >  2:3212152/1‑148 (MQ=255)
                                                                                            tttCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAg                                                       >  2:1044826/1‑148 (MQ=255)
                                                                                             ttCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGc                                                      >  1:3130841/1‑148 (MQ=255)
                                                                                               cAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCaa                                                    <  2:1576748/148‑1 (MQ=255)
                                                                                               cAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCaa                                                    >  1:2155467/1‑148 (MQ=255)
                                                                                                  cAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAAcgc                                                 >  2:523475/1‑148 (MQ=255)
                                                                                                    aaTAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGt                                               <  2:147249/148‑1 (MQ=255)
                                                                                                        gCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAg                                           >  1:1885558/1‑148 (MQ=255)
                                                                                                         cTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGc                                          >  2:1798149/1‑148 (MQ=255)
                                                                                                          tGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCa                                         >  2:981152/1‑148 (MQ=255)
                                                                                                          tGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCa                                         >  2:2635306/1‑148 (MQ=255)
                                                                                                           gCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCAt                                        <  2:294670/148‑1 (MQ=255)
                                                                                                           gCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCAt                                        <  1:1074697/148‑1 (MQ=255)
                                                                                                           gCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCAt                                        <  1:2376612/148‑1 (MQ=255)
                                                                                                              aaCGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCAtctc                                     <  1:1961506/148‑1 (MQ=255)
                                                                                                                  aCCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCcaca                                 <  2:158005/148‑1 (MQ=255)
                                                                                                                  aCCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCcaca                                 <  2:158119/148‑1 (MQ=255)
                                                                                                                         aaCATTCGATTCATATCGCTATGGAAATTTGGGTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCg                          >  2:3155839/1‑148 (MQ=255)
                                                                                                                         aaCATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCg                          >  2:896203/1‑148 (MQ=255)
                                                                                                                          aCATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGt                         >  2:3078655/1‑148 (MQ=255)
                                                                                                                          aCATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGt                         >  2:3078437/1‑148 (MQ=255)
                                                                                                                                  ttCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAc                 >  1:2674593/1‑148 (MQ=255)
                                                                                                                                  ttCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAc                 >  1:671453/1‑148 (MQ=255)
                                                                                                                                     atatCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCgt              <  1:2666774/148‑1 (MQ=255)
                                                                                                                                     atatCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCgt              <  2:121607/148‑1 (MQ=255)
                                                                                                                                       atCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCgtgt            <  1:1044826/148‑1 (MQ=255)
                                                                                                                                         cGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTAt          <  1:1556198/148‑1 (MQ=255)
                                                                                                                                         cGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAAGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTAt          >  2:745143/1‑148 (MQ=255)
                                                                                                                                             aTGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTccc      >  2:2614031/1‑148 (MQ=255)
                                                                                                                                                 cAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTg  >  2:3224902/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTG  >  CP009273/1820743‑1821035

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: