Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,850,839 G→T 43.0% 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,8390GT43.0% 244.5 / inf 279P309Q (CCG→CAG) yidJsulfatase/phosphatase superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (79/80);  new base T (53/67);  total (132/147)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 282 total aligned reads displayed.

GGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTC  >  CP009273/3850692‑3850987
                                                                                                                                                   |                                                                                                                                                    
 gCTCTTTCACGGCAAGGATATTGTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGa                                                                                                                                                    <  2:2716469/149‑1 (MQ=255)
            gCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGAt                                                                                                                                         <  1:3225323/149‑1 (MQ=255)
              aaGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATgt                                                                                                                                       <  2:2345578/149‑1 (MQ=255)
              aaGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATgt                                                                                                                                       <  2:2389744/149‑1 (MQ=255)
              aaGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATgt                                                                                                                                       <  1:902192/149‑1 (MQ=255)
                ggATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATgtcg                                                                                                                                     >  2:897364/1‑149 (MQ=255)
                   tatTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCa                                                                                                                                  <  1:3454780/149‑1 (MQ=255)
                      tttCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATAc                                                                                                                               >  1:1599794/1‑149 (MQ=255)
                        tCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACAt                                                                                                                             <  2:1335829/149‑1 (MQ=255)
                         cccccGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATc                                                                                                                            <  2:3325230/149‑1 (MQ=255)
                         cccccGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATc                                                                                                                            >  1:2650840/1‑149 (MQ=255)
                         cccccGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATc                                                                                                                            <  2:1448196/149‑1 (MQ=255)
                           cccGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATcgc                                                                                                                          >  2:1802333/1‑149 (MQ=255)
                               gCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCgccgcc                                                                                                                      >  1:1090746/1‑149 (MQ=255)
                                cAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGccc                                                                                                                     <  2:4300905/149‑1 (MQ=255)
                                 aGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGcccc                                                                                                                    <  2:2487394/149‑1 (MQ=255)
                                 aGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGcccc                                                                                                                    <  1:3903053/149‑1 (MQ=255)
                                  gAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCt                                                                                                                   <  2:1058945/149‑1 (MQ=255)
                                  gAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCt                                                                                                                   <  1:3095597/149‑1 (MQ=255)
                                  gAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCACCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCt                                                                                                                   >  1:2091964/1‑149 (MQ=255)
                                  gAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCACCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCt                                                                                                                   >  1:2092092/1‑149 (MQ=255)
                                  gAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCACCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCt                                                                                                                   >  1:1454950/1‑149 (MQ=255)
                                      ctctGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTAc                                                                                                               <  1:4670533/149‑1 (MQ=255)
                                       tctGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACt                                                                                                              <  1:42559/149‑1 (MQ=255)
                                           gTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATc                                                                                                          >  1:2963664/1‑149 (MQ=255)
                                              tttCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGc                                                                                                       <  2:485815/149‑1 (MQ=255)
                                               ttCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCt                                                                                                      >  1:1401177/1‑149 (MQ=255)
                                                   atatCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTAtg                                                                                                  <  2:4501360/149‑1 (MQ=255)
                                                   atatCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTAtg                                                                                                  <  2:2001209/149‑1 (MQ=255)
                                                    tatCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTAtgt                                                                                                 >  2:291891/1‑149 (MQ=255)
                                                     atCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTAtgtg                                                                                                >  2:245915/1‑149 (MQ=255)
                                                      tCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTgc                                                                                               <  1:3485101/149‑1 (MQ=255)
                                                         gCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGcc                                                                                            <  2:2468381/149‑1 (MQ=255)
                                                         gCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGcc                                                                                            <  2:4048842/149‑1 (MQ=255)
                                                           cAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCca                                                                                          >  2:667770/1‑149 (MQ=255)
                                                           cAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCca                                                                                          <  1:3590553/149‑1 (MQ=255)
                                                             gcgcCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCcatc                                                                                        <  2:4226627/149‑1 (MQ=255)
                                                              cgcCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCcatca                                                                                       <  1:1941153/149‑1 (MQ=255)
                                                              cgcCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCcatca                                                                                       <  1:4614148/149‑1 (MQ=255)
                                                               gcCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCcatcat                                                                                      <  1:4356019/149‑1 (MQ=255)
                                                               gcCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCcatcat                                                                                      <  2:722883/149‑1 (MQ=255)
                                                                ccATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCAtt                                                                                     >  2:3297693/1‑149 (MQ=255)
                                                                ccATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCAtt                                                                                     >  1:187487/1‑149 (MQ=255)
                                                                ccATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCAtt                                                                                     <  1:2598563/149‑1 (MQ=255)
                                                                ccATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCAtt                                                                                     >  1:230233/1‑149 (MQ=255)
                                                                ccATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCAtt                                                                                     >  1:2112758/1‑149 (MQ=255)
                                                                 catcatTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCAttt                                                                                    >  1:2478758/1‑149 (MQ=255)
                                                                 catcatTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCAttt                                                                                    >  2:4298664/1‑149 (MQ=255)
                                                                   tcatTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCg                                                                                  >  1:4083253/1‑149 (MQ=255)
                                                                       tGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGt                                                                              <  1:2669136/149‑1 (MQ=255)
                                                                       tGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGt                                                                              <  1:3349472/149‑1 (MQ=255)
                                                                        gTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTGTGTGCGCCCATCATTTCGCCGTg                                                                             >  2:747648/1‑149 (MQ=255)
                                                                             cAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCgg                                                                        <  1:102876/149‑1 (MQ=255)
                                                                               gTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAg                                                                      >  2:396137/1‑149 (MQ=255)
                                                                               gTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAg                                                                      >  1:1631082/1‑149 (MQ=255)
                                                                               gTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAg                                                                      <  1:3889437/149‑1 (MQ=255)
                                                                                tAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAgg                                                                     >  2:3256463/1‑149 (MQ=255)
                                                                                 aaaTCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGt                                                                    <  2:4572840/149‑1 (MQ=255)
                                                                                  aaTCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGta                                                                   >  2:398235/1‑149 (MQ=255)
                                                                                  aaTCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGta                                                                   >  1:853515/1‑149 (MQ=255)
                                                                                  aaTCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGta                                                                   >  1:4293390/1‑149 (MQ=255)
                                                                                  aaTCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGACGTGATCGGAGGta                                                                   >  1:1770008/1‑149 (MQ=255)
                                                                                   aTCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGtat                                                                  <  2:728563/149‑1 (MQ=255)
                                                                                   aTCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGtat                                                                  >  1:2153291/1‑149 (MQ=255)
                                                                                         atGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTataaataa                                                            <  2:1124535/149‑1 (MQ=255)
                                                                                           gactgactGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAAcc                                                          >  2:2670317/1‑149 (MQ=255)
                                                                                             ctgactgGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCa                                                        <  1:1829805/149‑1 (MQ=255)
                                                                                             ctgactgGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCa                                                        <  2:3384667/149‑1 (MQ=255)
                                                                                              tgactgGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCAc                                                       <  1:655667/149‑1 (MQ=255)
                                                                                                   ggCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTAtt                                                  >  2:2931734/1‑149 (MQ=255)
                                                                                                   ggCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTAtt                                                  >  2:4127326/1‑149 (MQ=255)
                                                                                                   ggCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTAtt                                                  >  2:3317137/1‑149 (MQ=255)
                                                                                                    gCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTAttt                                                 >  1:1013953/1‑149 (MQ=255)
                                                                                                      gTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTc                                               <  1:1251286/149‑1 (MQ=255)
                                                                                                       tATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCa                                              >  2:2509938/1‑149 (MQ=255)
                                                                                                        aTCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCAc                                             <  2:1837480/149‑1 (MQ=255)
                                                                                                           gcccTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGACGGCGGGTGATGTCGTCCGACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATACATAACCCACGTATTTTCACGtt                                          <  2:1747838/147‑1 (MQ=255)
                                                                                                           gACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGtt                                          >  2:3103098/1‑149 (MQ=255)
                                                                                                           gACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGtt                                          <  1:21969/149‑1 (MQ=255)
                                                                                                           gACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGtt                                          >  2:3727042/1‑149 (MQ=255)
                                                                                                           gACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGtt                                          >  2:1209079/1‑149 (MQ=255)
                                                                                                             ccTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGc                                        >  1:2421653/1‑149 (MQ=255)
                                                                                                              cTGTCGCCGCTCCCCTTGCGGCAAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGct                                       >  2:3830349/1‑149 (MQ=255)
                                                                                                                gTCGCCGCTCCCCTTGCGGCGCACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGCTCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATCACCCACGTATTTTCACGTTGctct                                     <  2:1297794/149‑1 (MQ=255)
                                                                                                                gTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGctct                                     <  2:906833/149‑1 (MQ=255)
                                                                                                                 tCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTg                                    <  1:3719284/149‑1 (MQ=255)
                                                                                                                    ccgcTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCg                                 <  1:2046916/149‑1 (MQ=255)
                                                                                                                    ccgcTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCg                                 >  1:3446097/1‑149 (MQ=255)
                                                                                                                    ccgcTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCg                                 >  1:1152547/1‑149 (MQ=255)
                                                                                                                     cgcTCCCCTTGCTGCGAACGGATGATCAGCGGAATGCGTGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCTTCATTTCGCCGTGATCGGAGGTATAAATATCCCACGTATTTTCACGTTGCTCTGGctt                                <  1:470752/149‑3 (MQ=255)
                                                                                                                        tCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTaa                             >  2:950554/1‑149 (MQ=255)
                                                                                                                          cccTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAgg                           >  2:3434390/1‑149 (MQ=255)
                                                                                                                           ccTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGc                          <  2:2112700/149‑1 (MQ=255)
                                                                                                                             ttGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCAt                        >  1:755871/1‑149 (MQ=255)
                                                                                                                              tGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCAtt                       <  1:1900019/149‑1 (MQ=255)
                                                                                                                                cggcgAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATtga                     <  1:153272/149‑1 (MQ=255)
                                                                                                                                     aaCGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGAcc                >  1:4156686/1‑149 (MQ=255)
                                                                                                                                           tgatCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCg          >  2:2169369/1‑149 (MQ=255)
                                                                                                                                            gatCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGa         <  1:578255/149‑1 (MQ=255)
                                                                                                                                                aGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTGATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTg     >  1:2919759/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTC  >  CP009273/3850692‑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: