Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,850,839 G→T 53.6% P309Q (CCG→CAG)  yidJ ← sulfatase/phosphatase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,850,8390GT53.6% 2.1 / inf 332P309Q (CCG→CAG) yidJsulfatase/phosphatase superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (71/83);  new base T (79/99);  total (150/182)
Fisher's exact test for biased strand distribution p-value = 8.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 334 total aligned reads displayed.

TGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTC  >  CP009273/3850691‑3850987
                                                                                                                                                    |                                                                                                                                                    
  gCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGa                                                                                                                                                    >  1:2778845/1‑149 (MQ=255)
    tctTTCACGGCAAGGATATTGTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTTTCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAAt                                                                                                                                                  <  1:2188871/149‑1 (MQ=255)
     ctTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATg                                                                                                                                                 >  1:3770446/1‑149 (MQ=255)
      tttCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGc                                                                                                                                                >  1:5156046/1‑149 (MQ=255)
      tttCACGGCAAGGATATTGTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGc                                                                                                                                                <  2:876220/149‑1 (MQ=255)
         cACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCggg                                                                                                                                             <  2:4118202/149‑1 (MQ=255)
         cACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCggg                                                                                                                                             <  1:3714673/149‑1 (MQ=255)
          aCGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGt                                                                                                                                            <  1:951024/149‑1 (MQ=255)
            ggCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGa                                                                                                                                          <  1:841473/149‑1 (MQ=255)
             gCAAGTATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGAt                                                                                                                                         <  2:4704477/149‑1 (MQ=255)
             gCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGAt                                                                                                                                         <  2:2841877/149‑1 (MQ=255)
             gCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGAt                                                                                                                                         >  1:1408544/1‑149 (MQ=255)
              cAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGACCGGATGATCAGCTGAATGCGGGTGATg                                                                                                                                        <  2:1906737/149‑1 (MQ=255)
              cAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATg                                                                                                                                        <  1:2733255/149‑1 (MQ=255)
                aGGATATTTTCCCCCGGCAGAATCGCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATgtc                                                                                                                                      >  1:3034481/1‑149 (MQ=255)
                    tatTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCa                                                                                                                                  <  1:1888458/149‑1 (MQ=255)
                      ttttCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATa                                                                                                                                <  1:3553367/149‑1 (MQ=255)
                       tttCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATAc                                                                                                                               <  2:967497/149‑1 (MQ=255)
                       tttCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATAc                                                                                                                               <  1:163799/149‑1 (MQ=255)
                           ccccGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATcg                                                                                                                           >  2:3240512/1‑149 (MQ=255)
                           ccccGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATcg                                                                                                                           >  1:2400655/1‑149 (MQ=255)
                            cccGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAAAGCGAGAGATGACGTCAAACAAcga                                                                                                                          >  2:4147494/1‑148 (MQ=255)
                                gCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCgccgcc                                                                                                                      >  1:3440266/1‑149 (MQ=255)
                                 cAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGccc                                                                                                                     <  2:165842/149‑1 (MQ=255)
                                 cAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGccc                                                                                                                     <  2:1475042/149‑1 (MQ=255)
                                  aGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGcccc                                                                                                                    <  1:2756872/149‑1 (MQ=255)
                                   gAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCt                                                                                                                   <  1:1041347/149‑1 (MQ=255)
                                   gAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCt                                                                                                                   <  2:4952765/149‑1 (MQ=255)
                                    aaTCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCtt                                                                                                                  >  1:2471291/1‑149 (MQ=255)
                                     aTCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCttt                                                                                                                 >  2:2909949/1‑149 (MQ=255)
                                     aTCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCttt                                                                                                                 <  1:2877813/149‑1 (MQ=255)
                                      tctcTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTa                                                                                                                <  2:2158332/149‑1 (MQ=255)
                                       ctctGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTAc                                                                                                               >  2:1168016/1‑149 (MQ=255)
                                        tctGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACt                                                                                                              <  2:4407485/149‑1 (MQ=255)
                                         ctGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTg                                                                                                             >  2:2520412/1‑149 (MQ=255)
                                          tGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGa                                                                                                            >  2:3485156/1‑149 (MQ=255)
                                           ggTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGAt                                                                                                           <  1:909594/149‑1 (MQ=255)
                                            gTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATc                                                                                                          <  2:2240442/149‑1 (MQ=255)
                                               tttCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGc                                                                                                       <  1:111400/149‑1 (MQ=255)
                                               tttCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGc                                                                                                       <  1:2139194/149‑1 (MQ=255)
                                               tttCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGc                                                                                                       <  2:333433/149‑1 (MQ=255)
                                                    atatCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTAtg                                                                                                  <  2:4019226/149‑1 (MQ=255)
                                                     tatCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTAtgt                                                                                                 >  1:4537971/1‑149 (MQ=255)
                                                     tatCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTAtgt                                                                                                 >  2:3671211/1‑149 (MQ=255)
                                                              gcgcCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCcatc                                                                                        >  2:2378162/1‑149 (MQ=255)
                                                               cgcCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCcatca                                                                                       <  2:1267157/149‑1 (MQ=255)
                                                                 ccATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCAtt                                                                                     >  2:4488675/1‑149 (MQ=255)
                                                                   atcatTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTc                                                                                   <  1:2255452/149‑1 (MQ=255)
                                                                    tcatTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCg                                                                                  >  1:4733021/1‑149 (MQ=255)
                                                                      atTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGcc                                                                                <  2:5220790/149‑1 (MQ=255)
                                                                        tGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGt                                                                              <  1:3126960/149‑1 (MQ=255)
                                                                         gTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTg                                                                             <  1:5053080/149‑1 (MQ=255)
                                                                         gTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTg                                                                             <  1:693762/149‑1 (MQ=255)
                                                                         gTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTg                                                                             >  2:4644896/1‑149 (MQ=255)
                                                                             gCAGTAAATCGATATGACTGATTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCg                                                                         <  2:2222854/149‑1 (MQ=255)
                                                                             gCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCg                                                                         <  2:1556900/149‑1 (MQ=255)
                                                                                gTACATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAg                                                                      <  1:4012722/149‑1 (MQ=255)
                                                                                 tAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAgg                                                                     >  1:1686642/1‑149 (MQ=255)
                                                                                  aaaTCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGt                                                                    <  1:539853/149‑1 (MQ=255)
                                                                                    aTCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGtat                                                                  <  1:2499937/149‑1 (MQ=255)
                                                                                     tCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGtata                                                                 <  2:3093244/149‑1 (MQ=255)
                                                                                     tCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGtata                                                                 >  1:2081819/1‑149 (MQ=255)
                                                                                     tCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCGTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGtata                                                                 <  2:2088239/149‑1 (MQ=255)
                                                                                       gATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTataaa                                                               <  2:3832608/149‑1 (MQ=255)
                                                                                         tatGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTataaata                                                             >  2:804260/1‑149 (MQ=255)
                                                                                           tgactgacTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAAc                                                           >  1:878004/1‑149 (MQ=255)
                                                                                                 actgGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGt                                                     >  2:4472320/1‑149 (MQ=255)
                                                                                                  ctgGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTa                                                    >  1:459525/1‑149 (MQ=255)
                                                                                                      cGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTAtttt                                                <  1:4979500/149‑1 (MQ=255)
                                                                                                          tCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACg                                            >  1:62475/1‑149 (MQ=255)
                                                                                                           cGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGt                                           <  1:2324239/149‑1 (MQ=255)
                                                                                                           cGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGt                                           <  1:3200600/149‑1 (MQ=255)
                                                                                                           cGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGt                                           <  1:4957333/149‑1 (MQ=255)
                                                                                                            gACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGtt                                          >  2:672047/1‑149 (MQ=255)
                                                                                                            gACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGtt                                          >  1:3294430/1‑149 (MQ=255)
                                                                                                            gACCTGGCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGtt                                          <  2:5081527/149‑1 (MQ=255)
                                                                                                             aCCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTg                                         >  1:3695445/1‑149 (MQ=255)
                                                                                                             aCCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTg                                         >  2:2682722/1‑149 (MQ=255)
                                                                                                                tGTCGCCGCTCCCCGTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGctc                                      <  2:3749749/149‑1 (MQ=255)
                                                                                                                  tCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTg                                    >  1:3607365/1‑149 (MQ=255)
                                                                                                                   cgccgcTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTgg                                   >  1:26642/1‑149 (MQ=255)
                                                                                                                   cgccgcTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTgg                                   >  1:3335930/1‑149 (MQ=255)
                                                                                                                    gccgcTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGc                                  >  2:1296631/1‑149 (MQ=255)
                                                                                                                       gcTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGTCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGtt                               >  2:4819463/1‑149 (MQ=255)
                                                                                                                       gcTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGtt                               <  2:3672185/149‑1 (MQ=255)
                                                                                                                         tCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTaa                             >  1:2643650/1‑149 (MQ=255)
                                                                                                                         tCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTaa                             >  1:1095061/1‑149 (MQ=255)
                                                                                                                         tCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTaa                             <  2:1392753/149‑1 (MQ=255)
                                                                                                                         tCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTaa                             >  1:3820415/1‑149 (MQ=255)
                                                                                                                           cccGTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAgg                           <  2:1166739/149‑1 (MQ=255)
                                                                                                                             cTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCa                         <  1:781167/149‑1 (MQ=255)
                                                                                                                                      aaCGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGAcc                >  2:2045169/1‑149 (MQ=255)
                                                                                                                                      aaCGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTGACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGAcc                <  2:3929664/149‑1 (MQ=255)
                                                                                                                                        cGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCg              >  1:4892943/1‑149 (MQ=255)
                                                                                                                                             gatCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGa         <  1:4826600/149‑1 (MQ=255)
                                                                                                                                             gatCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGa         <  1:4436181/149‑1 (MQ=255)
                                                                                                                                               tCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGAtt       >  2:4540805/1‑149 (MQ=255)
                                                                                                                                                cAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTGTGGGGTTAAGGCATTGATGACCCGTCCGAttt      <  1:4132503/149‑1 (MQ=255)
                                                                                                                                                 aGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTg     <  2:771389/149‑1 (MQ=255)
                                                                                                                                                    ggAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGtc  >  1:4039615/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTC  >  CP009273/3850691‑3850987

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: