Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,604,418 C→A 11.7% D174Y (GAC→TAC)  uxaB ← altronate oxidoreductase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,604,4180CA11.7% 1056.4 / 82.3 308D174Y (GAC→TAC) uxaBaltronate oxidoreductase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base C (130/142);  new base A (21/15);  total (151/157)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 323 total aligned reads displayed.

TTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCG  >  CP009273/1604269‑1604567
                                                                                                                                                     |                                                                                                                                                     
ttCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGtc                                                                                                                                                       >  1:3058881/1‑150 (MQ=255)
  cATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGtcaa                                                                                                                                                     >  2:1764369/1‑150 (MQ=255)
  cATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGtcaa                                                                                                                                                     >  2:2052290/1‑150 (MQ=255)
    tCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGtcaatc                                                                                                                                                   <  1:2673137/150‑1 (MQ=255)
     cgcgcgGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGtcaatca                                                                                                                                                  <  2:4690157/150‑1 (MQ=255)
            aTAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTcaca                                                                                                                                           <  2:2142274/150‑1 (MQ=255)
            aTAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTcaca                                                                                                                                           >  1:337827/1‑150 (MQ=255)
            aTAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTAAATCAACTcaca                                                                                                                                           >  1:1381659/1‑150 (MQ=255)
              aaCCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTAGCGCATAGCGCAGCACCAGGTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACg                                                                                                                                         >  1:350105/1‑150 (MQ=255)
                  ggTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGact                                                                                                                                     >  1:3145302/1‑148 (MQ=255)
                  ggTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGact                                                                                                                                     >  2:3527896/1‑148 (MQ=255)
                  ggTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAat                                                                                                                                     >  2:2473028/1‑150 (MQ=255)
                   gTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAatg                                                                                                                                    <  2:2753143/150‑1 (MQ=255)
                     aaCGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTAAATCAACTCACACGGAatgat                                                                                                                                  >  2:1491621/1‑150 (MQ=255)
                       cGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAAtgatga                                                                                                                                >  1:1292978/1‑150 (MQ=255)
                       cGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAAtgatga                                                                                                                                <  1:76682/150‑1 (MQ=255)
                        gATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATgatgat                                                                                                                               >  1:1771880/1‑150 (MQ=255)
                         aTACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATc                                                                                                                              >  1:4995847/1‑150 (MQ=255)
                            cGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTAAATCAACTCACACGGAATGATGATCCaa                                                                                                                           <  2:4998085/150‑1 (MQ=255)
                              gTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAAcc                                                                                                                         <  1:359651/150‑1 (MQ=255)
                                ccaccaGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCtt                                                                                                                       >  2:338726/1‑150 (MQ=255)
                                  accaGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTa                                                                                                                     <  2:2388652/150‑1 (MQ=255)
                                   ccaGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTAt                                                                                                                    >  1:3638402/1‑150 (MQ=255)
                                   ccaGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTAt                                                                                                                    <  1:473250/150‑1 (MQ=255)
                                      gCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCa                                                                                                                 <  1:125417/150‑1 (MQ=255)
                                       cGTAGAAAAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAg                                                                                                                >  2:737328/1‑150 (MQ=255)
                                        gTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAgc                                                                                                               <  2:2628612/150‑1 (MQ=255)
                                          aGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAgcgc                                                                                                             <  2:194037/150‑1 (MQ=255)
                                          aGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAgcgc                                                                                                             <  2:89270/150‑1 (MQ=255)
                                          aGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAgcgc                                                                                                             >  2:3174605/1‑150 (MQ=255)
                                          aGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAgcgc                                                                                                             >  1:4735378/1‑150 (MQ=255)
                                            aaCAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCAc                                                                                                           <  1:1470471/150‑1 (MQ=255)
                                            aaCAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCAc                                                                                                           <  1:1099195/150‑1 (MQ=255)
                                             aCAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCAGAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCAcc                                                                                                          <  2:1197723/150‑1 (MQ=255)
                                              cAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCg                                                                                                         >  2:3444904/1‑150 (MQ=255)
                                                gAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGtt                                                                                                       <  1:5041369/150‑1 (MQ=255)
                                                 aaGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTg                                                                                                      >  2:2719329/1‑150 (MQ=255)
                                                 aaGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTg                                                                                                      >  1:1040803/1‑150 (MQ=255)
                                                 aaGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTg                                                                                                      >  2:4836711/1‑150 (MQ=255)
                                                  aGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGa                                                                                                     >  1:540645/1‑150 (MQ=255)
                                                     tGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTATATCCAGCGCACCGTGGAAAt                                                                                                  >  2:4660490/1‑150 (MQ=255)
                                                      gTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATg                                                                                                 <  2:3256625/150‑1 (MQ=255)
                                                       ttAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATgg                                                                                                >  2:4529282/1‑150 (MQ=255)
                                                       ttAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATgg                                                                                                >  2:2934400/1‑150 (MQ=255)
                                                        tAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTAAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGc                                                                                               <  2:4358556/150‑1 (MQ=255)
                                                             tGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAg                                                                                          >  2:4568023/1‑150 (MQ=255)
                                                             tGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAg                                                                                          <  2:460613/150‑1 (MQ=255)
                                                               aTCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCg                                                                                        <  1:2341048/150‑1 (MQ=255)
                                                               aTCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCAGAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCg                                                                                        <  2:2212582/150‑1 (MQ=255)
                                                                tCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTGATCCAGCGCACCGTTGAAATGGCTAAAGCGt                                                                                       <  2:2420506/150‑1 (MQ=255)
                                                                  cAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTc                                                                                     <  2:540645/150‑1 (MQ=255)
                                                                  cAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTc                                                                                     >  1:2805730/1‑150 (MQ=255)
                                                                   aGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCg                                                                                    >  2:3576474/1‑150 (MQ=255)
                                                                     ccATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGaa                                                                                  <  2:1725754/150‑1 (MQ=255)
                                                                     ccATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGaa                                                                                  >  1:415023/1‑150 (MQ=255)
                                                                      cATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAAc                                                                                 <  2:3897308/150‑1 (MQ=255)
                                                                       aTTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAAca                                                                                <  2:1449286/150‑1 (MQ=255)
                                                                          gAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAAcagca                                                                             <  2:3975389/150‑1 (MQ=255)
                                                                          gAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAAcagca                                                                             <  2:2302332/150‑1 (MQ=255)
                                                                            aTAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGa                                                                           >  1:387393/1‑150 (MQ=255)
                                                                                aTGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACggg                                                                       <  2:1453620/150‑1 (MQ=255)
                                                                                   cTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCa                                                                    <  1:4647844/150‑1 (MQ=255)
                                                                                          cAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTAAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGc                                                             >  2:4140888/1‑150 (MQ=255)
                                                                                            gTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTGATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCg                                                           <  2:3125413/150‑1 (MQ=255)
                                                                                             tGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCgg                                                          >  1:4002377/1‑150 (MQ=255)
                                                                                             tGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCgg                                                          <  2:3058881/150‑1 (MQ=255)
                                                                                                ccACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATa                                                       >  2:1121668/1‑150 (MQ=255)
                                                                                                ccACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATa                                                       <  1:1347357/150‑1 (MQ=255)
                                                                                                ccACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATa                                                       >  1:4291250/1‑150 (MQ=255)
                                                                                                ccACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTAAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATa                                                       <  2:846463/150‑1 (MQ=255)
                                                                                                   ctctTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCt                                                    >  1:1642362/1‑150 (MQ=255)
                                                                                                    tctTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCtt                                                   >  2:4301867/1‑150 (MQ=255)
                                                                                                    tctTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCtt                                                   >  2:4754/1‑150 (MQ=255)
                                                                                                     ctTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTa                                                  >  2:774175/1‑150 (MQ=255)
                                                                                                      ttGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTAc                                                 >  1:3767647/1‑150 (MQ=255)
                                                                                                       tGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTAcc                                                >  1:1007330/1‑150 (MQ=255)
                                                                                                          gcATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCt                                             >  1:2019816/1‑150 (MQ=255)
                                                                                                           cagaGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTg                                            <  1:4097331/147‑1 (MQ=255)
                                                                                                            aTAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTgg                                           <  1:3837431/150‑1 (MQ=255)
                                                                                                             tAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCACCTCACACGGAATGATGATCCAACCTTTAGCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGgc                                          <  1:4660490/150‑1 (MQ=255)
                                                                                                              aGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGgcg                                         <  2:4228407/150‑1 (MQ=255)
                                                                                                                    gcaCCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCAt                                   <  1:4863047/150‑1 (MQ=255)
                                                                                                                       ccAGTTCACGCAACGCGTCGCCATTGTAGTAAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGa                                >  2:3439285/1‑150 (MQ=255)
                                                                                                                        cAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGaa                               >  2:904215/1‑150 (MQ=255)
                                                                                                                         aGTTCACGCAACGCGTCGCCATTGTAGTAAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAAt                              <  2:1381659/150‑1 (MQ=255)
                                                                                                                          gTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAAtt                             >  1:4360715/1‑150 (MQ=255)
                                                                                                                          gTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAAtt                             <  1:686395/150‑1 (MQ=255)
                                                                                                                           ttCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTGATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAAttt                            <  1:160065/150‑1 (MQ=255)
                                                                                                                            tCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTa                           <  1:2253089/150‑1 (MQ=255)
                                                                                                                             cACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTAt                          >  1:2378622/1‑150 (MQ=255)
                                                                                                                              aCGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATc                         <  2:972159/150‑1 (MQ=255)
                                                                                                                               cGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCg                        >  1:1203928/1‑150 (MQ=255)
                                                                                                                                   aCGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCg                    <  2:4442033/150‑1 (MQ=255)
                                                                                                                                         cGCCATTGTAGTAAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGt              >  1:3974452/1‑150 (MQ=255)
                                                                                                                                           ccATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAg            >  1:4939471/1‑150 (MQ=255)
                                                                                                                                            cATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGc           >  1:3451452/1‑150 (MQ=255)
                                                                                                                                             aTTGTAGTAAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCt          >  2:1470189/1‑150 (MQ=255)
                                                                                                                                                 tAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATa      <  2:3701905/150‑1 (MQ=255)
                                                                                                                                                   gTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATAcc    >  1:3012419/1‑150 (MQ=255)
                                                                                                                                                     caatcaaCTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCg  >  2:4795042/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCG  >  CP009273/1604269‑1604567

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: