Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 39,866 G→T 6.2% S174Y (TCT→TAT)  caiA ← crotonobetaine reductase subunit II, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927339,8660GT 6.2% 1102.8 / 36.3 293S174Y (TCT→TAT) caiAcrotonobetaine reductase subunit II, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base G (106/168);  new base T (7/11);  total (113/180)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 299 total aligned reads displayed.

CGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGT  >  CP009273/39717‑40014
                                                                                                                                                     |                                                                                                                                                    
       gTCTTTCTCGTCCAGGTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGcccc                                                                                                                                               <  2:2455521/150‑1 (MQ=255)
          tttCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATCTCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGACGCCCCGTc                                                                                                                                            >  2:1191855/1‑150 (MQ=255)
          tttCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTc                                                                                                                                            <  2:956735/150‑1 (MQ=255)
          tttCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTc                                                                                                                                            >  1:4988421/1‑150 (MQ=255)
          tttCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTc                                                                                                                                            <  2:40371/150‑1 (MQ=255)
              tcGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTcgcgc                                                                                                                                        >  2:3739028/1‑150 (MQ=255)
                  ccAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCa                                                                                                                                    <  2:1512741/150‑1 (MQ=255)
                  ccAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCa                                                                                                                                    <  2:1344927/150‑1 (MQ=255)
                   cAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAATCCCCGTCGCGCGCCAt                                                                                                                                   >  2:3225044/1‑150 (MQ=255)
                      ttCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGATAAGCCCCGTCGCGCGCCATcac                                                                                                                                <  1:1579609/150‑1 (MQ=255)
                       tCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATcccc                                                                                                                               >  2:4907869/1‑150 (MQ=255)
                         cACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAACCCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATcaccac                                                                                                                             <  1:2756695/150‑1 (MQ=255)
                         cACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATcaccac                                                                                                                             >  1:3356403/1‑150 (MQ=255)
                          aCGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACg                                                                                                                            <  2:3321836/150‑1 (MQ=255)
                            gtcgtcAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGAt                                                                                                                          >  1:88810/1‑150 (MQ=255)
                             tcgtcAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATg                                                                                                                         <  2:2551302/150‑1 (MQ=255)
                              cgtcAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGt                                                                                                                        >  1:344826/1‑150 (MQ=255)
                               gtcAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTa                                                                                                                       <  2:456652/150‑1 (MQ=255)
                               gtcAAAGGTGATTACACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTa                                                                                                                       >  1:3579821/1‑150 (MQ=255)
                                  aaaGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACgg                                                                                                                    >  1:1898485/1‑150 (MQ=255)
                                    aGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACgggg                                                                                                                  >  1:2592399/1‑150 (MQ=255)
                                    aGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACgggg                                                                                                                  >  1:3689325/1‑150 (MQ=255)
                                      gTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGgtg                                                                                                                >  2:3040045/1‑150 (MQ=255)
                                      gTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGgtg                                                                                                                <  2:4315502/150‑1 (MQ=255)
                                        gATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGATAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTa                                                                                                              >  1:767450/1‑150 (MQ=255)
                                        gATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGATAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTa                                                                                                              >  1:3392174/1‑150 (MQ=255)
                                         atgtctcagcagCTCTCCCTCCGCCGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCATCACACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCGATCACCACGATGTACGGGGTGTAg                                                                                                             <  2:3806326/144‑1 (MQ=255)
                                            tCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGgcg                                                                                                          <  1:1115311/150‑1 (MQ=255)
                                               cagcagCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTGGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCgctg                                                                                                       <  1:803844/150‑1 (MQ=255)
                                                     cTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGta                                                                                                 <  1:1803963/150‑1 (MQ=255)
                                                      tATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGtaa                                                                                                <  2:4055484/150‑1 (MQ=255)
                                                        tCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGtaata                                                                                              <  2:767137/150‑1 (MQ=255)
                                                        tCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGtaata                                                                                              <  2:1614626/150‑1 (MQ=255)
                                                         ccATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCCCCACGATGTACGGGGTGTAGGCGCTGCTGGtaataa                                                                                             <  2:1824599/150‑1 (MQ=255)
                                                         ccATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGtaataa                                                                                             <  2:1824748/150‑1 (MQ=255)
                                                          cATACGCAGACCGAGCTTGTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATaaa                                                                                            <  1:4281029/150‑1 (MQ=255)
                                                               gCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACt                                                                                       >  2:443336/1‑150 (MQ=255)
                                                                cAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACtt                                                                                      <  2:2951658/150‑1 (MQ=255)
                                                                 aGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTta                                                                                     <  1:4316663/150‑1 (MQ=255)
                                                                 aGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTta                                                                                     >  1:2894158/1‑150 (MQ=255)
                                                                 aGACCGAGCTTTGCCAGTGTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTta                                                                                     <  2:2686263/150‑1 (MQ=255)
                                                                   aCCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTtact                                                                                   >  2:2794222/1‑150 (MQ=255)
                                                                     cGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTtactac                                                                                 <  1:3032371/150‑1 (MQ=255)
                                                                          ttttCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGATAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTa                                                                            >  2:3967094/1‑150 (MQ=255)
                                                                               aaGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATa                                                                       <  1:2990894/150‑1 (MQ=255)
                                                                                  tttGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAAt                                                                    <  1:2288860/150‑1 (MQ=255)
                                                                                  tttGGTCACTGTGATGCCCGTTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAAt                                                                    <  2:3129630/150‑1 (MQ=255)
                                                                                     ggTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCtt                                                                 <  2:5238162/150‑1 (MQ=255)
                                                                                      gTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTa                                                                <  1:1576978/150‑1 (MQ=255)
                                                                                        cACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTAcc                                                              >  1:2005642/1‑150 (MQ=255)
                                                                                        cACTTGGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTAcc                                                              <  1:2885250/150‑1 (MQ=255)
                                                                                          cTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCAt                                                            <  1:551367/150‑1 (MQ=255)
                                                                                          cTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCAt                                                            >  1:615710/1‑150 (MQ=255)
                                                                                           ttggATGCCCGGTGTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCAtt                                                           <  1:2414851/147‑1 (MQ=255)
                                                                                            ttGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCAttt                                                          <  2:139680/150‑1 (MQ=255)
                                                                                            ttGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCAttt                                                          >  2:2359254/1‑150 (MQ=255)
                                                                                             tGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTc                                                         >  1:4536207/1‑150 (MQ=255)
                                                                                             tGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTc                                                         >  1:1353711/1‑150 (MQ=255)
                                                                                             tGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTc                                                         <  1:4032032/150‑1 (MQ=255)
                                                                                             tGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGGGTAGGCGCTGCTGGTAATAAAACGCTTACTCCCATTCAGATCTATCTTACCATTTc                                                         >  2:1492782/1‑150 (MQ=255)
                                                                                              gATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCt                                                        >  1:2304613/1‑150 (MQ=255)
                                                                                                 gCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACg                                                     >  1:2357706/1‑150 (MQ=255)
                                                                                                   ccGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACggg                                                   <  2:4403664/150‑1 (MQ=255)
                                                                                                   ccGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACggg                                                   <  2:4622836/150‑1 (MQ=255)
                                                                                                      gTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGtat                                                <  1:1035843/150‑1 (MQ=255)
                                                                                                      gTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGACATCACCACGATTTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCTTTTCTTCGGGtat                                                <  2:4642177/150‑1 (MQ=255)
                                                                                                         tGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAg                                             >  1:4614258/1‑150 (MQ=255)
                                                                                                         tGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAg                                             >  1:3404357/1‑150 (MQ=255)
                                                                                                           cTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAATCACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTc                                           >  2:12269/1‑150 (MQ=255)
                                                                                                            tCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCg                                          >  2:216129/1‑150 (MQ=255)
                                                                                                             cATATCAACAACCCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGt                                         <  1:5206374/150‑1 (MQ=255)
                                                                                                              atatCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCTCCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGtt                                        <  1:4643668/150‑1 (MQ=255)
                                                                                                              atatCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGtt                                        <  1:2020756/150‑1 (MQ=255)
                                                                                                                atCAACAAACCATTCGGTGTAGACTGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGtttt                                      <  1:4922203/150‑1 (MQ=255)
                                                                                                                atCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGtttt                                      <  1:4173760/150‑1 (MQ=255)
                                                                                                                atCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGACGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGtttt                                      <  2:719353/150‑1 (MQ=255)
                                                                                                                  caacaaACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCa                                    <  2:3431285/150‑1 (MQ=255)
                                                                                                                  caacaaACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCa                                    <  1:1222679/150‑1 (MQ=255)
                                                                                                                   aacaaACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAg                                   >  1:4804282/1‑150 (MQ=255)
                                                                                                                    acaaACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAgg                                  <  1:745962/150‑1 (MQ=255)
                                                                                                                      aaaCCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCt                                <  2:4064657/150‑1 (MQ=255)
                                                                                                                         ccATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTAcc                             <  1:4734432/150‑1 (MQ=255)
                                                                                                                          cATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTAccc                            <  1:2688632/150‑1 (MQ=255)
                                                                                                                             tCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACg                         <  2:4988421/150‑1 (MQ=255)
                                                                                                                             tCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACg                         >  2:4175132/1‑150 (MQ=255)
                                                                                                                               ggTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTc                       >  1:3528135/1‑150 (MQ=255)
                                                                                                                                    aGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGc                  >  1:4241501/1‑150 (MQ=255)
                                                                                                                                      aCAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGccc                >  1:382522/1‑150 (MQ=255)
                                                                                                                                      aCAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGccc                <  1:363472/150‑1 (MQ=255)
                                                                                                                                      aCAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGccc                >  2:3880254/1‑150 (MQ=255)
                                                                                                                                       cAGGTTTGTCCGGATAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCcg               <  2:2056374/150‑1 (MQ=255)
                                                                                                                                       cAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCcg               <  1:3739028/150‑1 (MQ=255)
                                                                                                                                        aGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCcgc              >  1:3271368/1‑150 (MQ=255)
                                                                                                                                        aGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTAAGGCTACCCACGTCGGAGCCcgc              >  1:4450445/1‑150 (MQ=255)
                                                                                                                                            ttGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGccc          <  2:3574683/150‑1 (MQ=255)
                                                                                                                                                 cGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTc     <  2:3820035/150‑1 (MQ=255)
                                                                                                                                                 cGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTc     <  2:1020880/150‑1 (MQ=255)
                                                                                                                                                  ggAGCAGCCCCGTCGCGCGCCATCACCACGAGGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCg    <  2:4294282/150‑1 (MQ=255)
                                                                                                                                                  ggAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTCCGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCg    <  2:248197/150‑1 (MQ=255)
                                                                                                                                                    agaAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGt  <  2:4722267/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGT  >  CP009273/39717‑40014

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: