Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,514,893 A→C 8.2% Q32P (CAG→CCG)  urtB → urea ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,514,8930AC 8.2% 683.2 / 7.6 281Q32P (CAG→CCG) urtBurea ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (131/104);  new base C (10/19);  total (152/129)
Fisher's exact test for biased strand distribution p-value = 4.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 285 total aligned reads displayed.

TGTAGGAGCGGGCTCACCCGCGAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGATCGACAACGCCCTGG  >  NC_002947/5514744‑5515042
                                                                                                                                                     |                                                                                                                                                     
 gtAGGAGCGGGCTCACCCGCGAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAg                                                                                                                                                      <  1:5233451/150‑1 (MQ=255)
  tAGGAGCGGGCTCACCCGCGAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAgg                                                                                                                                                     <  1:1310257/150‑1 (MQ=255)
  tAGGAGCGGGCTCACCCGCGAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAgg                                                                                                                                                     <  1:1809467/150‑1 (MQ=255)
      aGCGGGCTCACCCGCGAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAggccgg                                                                                                                                                 >  2:4243316/1‑150 (MQ=255)
         gggCTCACCCGCGAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCt                                                                                                                                              <  2:388080/150‑1 (MQ=255)
           gCTCACCCGCGAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGc                                                                                                                                            <  1:2851600/150‑1 (MQ=255)
           gCTCACCCGCGAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCCGCCTGc                                                                                                                                            >  2:5017583/1‑150 (MQ=255)
            cTCACCCGCGAAGAATCCAACACCGATGTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCt                                                                                                                                           <  2:306279/150‑1 (MQ=255)
             tCACCCGCGAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCCGGCCGGCCTGcgc                                                                                                                                          >  1:4725932/1‑148 (MQ=255)
             tCACCCGCGAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTc                                                                                                                                          >  2:4557054/1‑150 (MQ=255)
             tCACCCGCGAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTc                                                                                                                                          <  2:1929459/150‑1 (MQ=255)
             tCACCCGCGAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCGGCTc                                                                                                                                          >  1:950918/1‑150 (MQ=255)
                 ccGCGAAGAATCCAACCCCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCGGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCggcg                                                                                                                                      <  2:3691164/150‑4 (MQ=255)
                 ccGCGAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCGGGCCGGGCTGCGCGCAg                                                                                                                                      >  1:1913791/1‑150 (MQ=255)
                 ccGCGAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAg                                                                                                                                      <  1:3613703/150‑1 (MQ=255)
                  cgcgAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAgg                                                                                                                                     >  2:4677834/1‑150 (MQ=255)
                   gcgAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCGGCGCGAAGGc                                                                                                                                    >  1:908063/1‑150 (MQ=255)
                      aataaTCCAACACCGATTTCCAGTATTGCCCGAATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGAAACCCACGCCAGCGAGGGCGAATTCTTTCTAACCGCAAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTgg                                                                                                                                 <  2:713316/147‑1 (MQ=255)
                         aaTCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCGGGCCGGCCGGCGCGAAGGCGGGgca                                                                                                                              >  1:3601904/1‑150 (MQ=255)
                           tCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGgcagc                                                                                                                            >  2:5403315/1‑150 (MQ=255)
                           tCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGgcagc                                                                                                                            <  2:1278275/150‑1 (MQ=255)
                           tCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGgcagc                                                                                                                            <  2:4890427/150‑1 (MQ=255)
                              aaCACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGcaca                                                                                                                         >  2:4846242/1‑150 (MQ=255)
                              aaCACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCGGCGCGCAGGCGGGGCCGccca                                                                                                                         >  1:5347562/1‑150 (MQ=255)
                                cacCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAAc                                                                                                                       >  1:5614330/1‑150 (MQ=255)
                                  ccGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAAccc                                                                                                                     >  2:273855/1‑150 (MQ=255)
                                    gttttCCTGGAGTGCCCGCATGCTCTGACTCCTGCTCACCCTCCTGCTGATAAGCCCCCTGGCAACCCACGCCAGCGTGGGCTATTTCTTTCTCCCCGCCACGCCCGCCGAGCCGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGa                                                                                                                   <  1:3301255/148‑1 (MQ=255)
                                    gATTGCCAGTATTGCCCGCATGCTCAGACTCCGGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCCGCGAGGGCTAATTCTTTCTCACCGCCACGCCCGCCGCGCAGGCCGGCCGGCGCGAAGGCTGGGCAGCACAACCCGa                                                                                                                   <  1:4523739/150‑1 (MQ=255)
                                     aTTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGAt                                                                                                                  <  2:3115914/150‑1 (MQ=255)
                                                tGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCt                                                                                                       <  2:555533/150‑1 (MQ=255)
                                                   ccGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCGGGCCCGCCGGCTCGCAGGCTGGGCAGCACAAACCCAGGCCGCGCGCCTGcc                                                                                                    >  1:1692571/1‑150 (MQ=255)
                                                   ccGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCAGGGCAGCACAACCCCATGCCGCGCGCCTGcc                                                                                                    >  1:2734822/1‑150 (MQ=255)
                                                    cGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCCATGCCGCGCGCCAGCCg                                                                                                   >  1:1460521/1‑150 (MQ=255)
                                                     gCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCgc                                                                                                  <  2:4729283/150‑1 (MQ=255)
                                                         gCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCCAAGCCGCGCGCCAGccccggcg                                                                                              >  1:3850729/1‑144 (MQ=255)
                                                           tCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCGGGCCGGCCGGCGCGAAGGCGGGGCCGCACCACCCGCAGCCGCCCGCCGGCCCCTGCggg                                                                                            >  1:2778383/1‑150 (MQ=255)
                                                           tCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTgg                                                                                            >  2:3014857/1‑150 (MQ=255)
                                                             aGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGACCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGaa                                                                                          >  1:3646848/1‑150 (MQ=255)
                                                                cTCCTGATCACCCCCCTGCTGCTAACCCCCCTGCCAACCCACGCCAGCGAGGGCGAATTCTTCCTCACCCCCACGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAAc                                                                                       <  1:4050740/150‑1 (MQ=255)
                                                                 tCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCCGGCCGGCCGGCACGAAGGCTGGGCCGCCCAACCCCAAGCCGCGCGCCGGCCGCTCCAGGAAAAcc                                                                                      >  1:3148588/1‑150 (MQ=255)
                                                                 tCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCGCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCAGCCGCGGCTGGAAAAcc                                                                                      >  1:1850194/1‑150 (MQ=255)
                                                                 tCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCGGCTCGAAGGCTGGGCAGCACAACCCCAAGCCGCGCGCCCGCCCCAGCGGGAAAAcc                                                                                      >  1:4290778/1‑150 (MQ=255)
                                                                   cTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCCATGCCGCGCGCCAGCCGCTGCTGGGAAAccag                                                                                    >  1:3438016/1‑148 (MQ=255)
                                                                    tGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCGGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCGGCCGCAGCTGGgaaaacagc                                                                                   >  1:1126442/1‑145 (MQ=255)
                                                                    ggCTCACCCTCCTGCTGCTAAGCCCCCTGGCACCCCCCGCCAGCGAGGGCGATTTCCTTCCCACCGCCACGCCCCCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTgc                                                                                   <  1:2537903/149‑1 (MQ=255)
                                                                      cTCACCCTCCATCTGCTAAACCCCCTGCCACCCCACCCCCGCGAGGGCGAATTCTTCCACCCCGCCACGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTgcgc                                                                                 <  1:757008/150‑1 (MQ=255)
                                                                          cccTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAg                                                                             <  2:4084833/150‑1 (MQ=255)
                                                                              cctGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCTTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCg                                                                         <  2:4701926/150‑1 (MQ=255)
                                                                                  ctgctAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATc                                                                     <  2:3756893/150‑1 (MQ=255)
                                                                                      tAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCAGCTCGAAGGCTGGGCAGCACAACCCCAAGCCGCGCGCCTGCCGCTGCTGGAAAACCAGCGCCAAGGCcgccacgcca                                                                 >  1:5003418/1‑150 (MQ=255)
                                                                                       aaTCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGccac                                                                <  2:2214365/150‑1 (MQ=255)
                                                                                       aaTCCCCCTCGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGccac                                                                <  2:207788/150‑1 (MQ=255)
                                                                                       aaTCCCCCGGGCACCCCCCGCCCGCGTGGGCGAATTCTTTCCCCCCGCCACGCCCGCCGAGCAGGCGGGCCGGCGCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGccac                                                                <  1:870847/150‑1 (MQ=255)
                                                                                          cccccTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCaa                                                             <  2:443650/150‑1 (MQ=255)
                                                                                            cccTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCCGGCCGGCCGGCTCGACGGCGGGGCCGCCCCACCCGCAGCCCCGCGCCGGCCGCGGCGGGAAAACCGGCGCCCAGGCCGCCGCGccccccaac                                                           >  1:2059675/1‑146 (MQ=255)
                                                                                            cccTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCGCGACGGCGGGGCCGCACAACCCGCGCCCCCGCGCCGGCCCCGGCTGGCCAACCTCCCCCCAGGCCGCAGCGccccccaag                                                           >  1:4891149/1‑146 (MQ=255)
                                                                                            cccTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCGGCGCGAAGGCTGGGGAGCACAACCCCATGCCGCGCGCCTGCCGCGGCTGGAAAACCAGCGCCAAGGCCGCCACCCCACCaaag                                                           >  1:3758813/1‑148 (MQ=255)
                                                                                                   aaCCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAgcg                                                    >  2:3276731/1‑150 (MQ=255)
                                                                                                      ccACGCCTGCGTGGGCGAATTCTTTCACACCGCCACGCCCGCCGAGCCGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCaa                                                 <  1:3375629/150‑1 (MQ=255)
                                                                                                        aCGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAgg                                               >  2:3125852/1‑150 (MQ=255)
                                                                                                             aGCGCGGGCTAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTgcgc                                          <  1:5509700/150‑1 (MQ=255)
                                                                                                             aGCGAGGGCGATTTCTTTCTCCCCGCCACGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTgcgc                                          <  1:1017472/150‑1 (MQ=255)
                                                                                                             aGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTgcgc                                          <  2:5614330/150‑1 (MQ=255)
                                                                                                             aGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTgcgc                                          >  2:5429074/1‑150 (MQ=255)
                                                                                                                gtgggCGAATTCTTTCTCCCCGCCCAGCCCGCCGCGCCGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACACGCGCAAGGTGCGCCTg                                       <  1:3366835/148‑1 (MQ=255)
                                                                                                                    gCGAATTCTTTCTCACAGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAAACCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGaaca                                   >  2:2786955/1‑150 (MQ=255)
                                                                                                                     cGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGCTGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGTCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGaccac                                  >  2:1992538/1‑149 (MQ=255)
                                                                                                                      gatttctttctCCCCGCCACGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGaacaac                                 <  1:4629577/147‑1 (MQ=255)
                                                                                                                       aaTTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAATGCGGGGCAGCACAACCCGAAGCCGCGCGCCTGCCGCTGATGGAAAACATGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAATGTGCGCCTGAACAAcc                                >  1:446520/1‑150 (MQ=255)
                                                                                                                       aaTTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAAcc                                <  2:3374969/150‑1 (MQ=255)
                                                                                                                       aaTTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAAcc                                <  2:1496849/150‑1 (MQ=255)
                                                                                                                          tctttctCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAGACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGACCAACCGcc                             >  2:4224454/1‑150 (MQ=255)
                                                                                                                          tctttctCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGcc                             >  2:4224608/1‑150 (MQ=255)
                                                                                                                            tttctCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTg                           >  2:1595121/1‑150 (MQ=255)
                                                                                                                            tttctCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTg                           >  2:3855428/1‑150 (MQ=255)
                                                                                                                            tttctCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCAGCGCCAAGGCCGCCTCGCCACCAAAGACAAGCGCAAGGTGCGCCTGAACAACCGCCTg                           >  1:443643/1‑150 (MQ=255)
                                                                                                                            tttctCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAACGCGGGGCCGCACAACCCGATGCCGCCCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCAGCGCCACCCATGACAAGCGCAAGGTGCGCCGGAACAAACGccgg                           >  1:3026980/1‑148 (MQ=255)
                                                                                                                             ttctCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCCCATCGCCACCCATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTgc                          >  2:2588745/1‑150 (MQ=255)
                                                                                                                                tcACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGcgcg                       >  2:1376721/1‑150 (MQ=255)
                                                                                                                                tcACCGCCAAGCCCGCCGAGCAGGCCGGCCAGCTCGAAGGCTGGGCAGCACAACCCCATGCCGCGCGCCAGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAAAGACAAGCGCAAGGTGCGCCTGAACAACCCCCTGcgcg                       >  1:1579155/1‑150 (MQ=255)
                                                                                                                                  cccgccccgccCGCCGCGCCGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGc                     <  1:4243321/149‑1 (MQ=255)
                                                                                                                                  aCCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGc                     >  2:4670348/1‑150 (MQ=255)
                                                                                                                                  aCCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGc                     >  2:2169095/1‑150 (MQ=255)
                                                                                                                                  aCCCCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGc                     <  1:968634/150‑1 (MQ=255)
                                                                                                                                    cGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCt                   <  2:5167599/150‑1 (MQ=255)
                                                                                                                                      ccACGCCCGCCGCGCCGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGa                 <  1:5545112/150‑1 (MQ=255)
                                                                                                                                      ccAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGa                 <  2:792576/150‑1 (MQ=255)
                                                                                                                                      ccAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGa                 <  2:839505/150‑1 (MQ=255)
                                                                                                                                      ccAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGa                 <  2:3780229/150‑1 (MQ=255)
                                                                                                                                      ccAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGa                 <  2:2373562/150‑1 (MQ=255)
                                                                                                                                        aaGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGATc               <  2:3601904/150‑1 (MQ=255)
                                                                                                                                        aaGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGATc               <  2:2979698/150‑1 (MQ=255)
                                                                                                                                           cccGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCAGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAAAGACAACCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGATCGAc            >  1:4195820/1‑150 (MQ=255)
                                                                                                                                           cccGCCGAGCAGGCCGGCCAGCTCGAAGGCTGGGCAGCACAACACCATGCCCCGCGCCAGCCGCTGCAGGAAAAACAGCGCCAAGGCCGCAACGCCACCAAACACAAAACCAAGGTGCGCCAGAAAAACCCCCAGCGCGGCCAGATCGAc            >  1:1069831/1‑150 (MQ=255)
                                                                                                                                             cgccgAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGATCGACaa          <  2:2468564/150‑1 (MQ=255)
                                                                                                                                             cgccgAGCAGGCCGGCCTGCTCGAAGGCTAGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGATCGACaa          <  2:2405685/150‑1 (MQ=255)
                                                                                                                                              ttcgAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGATCGACAAc         >  2:1906061/3‑150 (MQ=255)
                                                                                                                                              ttcgAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGATCGACAAc         >  2:1692943/3‑150 (MQ=255)
                                                                                                                                                    cAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCCATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGATCGACAACGCCCTg   <  2:5299207/150‑1 (MQ=255)
                                                                                                                                                     aGGCCGGCCTTCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTgcgc                                          <  1:2436725/110‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTAGGAGCGGGCTCACCCGCGAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGATCGACAACGCCCTGG  >  NC_002947/5514744‑5515042

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