Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,067,754 T→A 20.7% H189L (CAC→CTC)  rutB ← ureidoacrylate amidohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,067,7540TA20.7% 396.2 / 100.5 174H189L (CAC→CTC) rutBureidoacrylate amidohydrolase
Reads supporting (aligned to +/- strand):  ref base T (71/67);  new base A (20/16);  total (91/83)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 176 total aligned reads displayed.

TTGGCATCGTTAAACTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTAT  >  CP009273/1067609‑1067898
                                                                                                                                                 |                                                                                                                                                
 tGGCATCGTTAAACTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAg                                                                                                                                               <  2:1215885/148‑1 (MQ=255)
 tGGCATCGTTAAACTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAg                                                                                                                                               <  1:2437035/148‑1 (MQ=255)
 tGGCATCGTTAAACTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAg                                                                                                                                               <  1:1963963/148‑1 (MQ=255)
  ggCATCGTTAAACTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGt                                                                                                                                              <  1:3012114/148‑1 (MQ=255)
  ggCATCGTTAAACTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGt                                                                                                                                              <  2:1387292/148‑1 (MQ=255)
   gCATCGTTAAACTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTGTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGtt                                                                                                                                             >  2:3399899/1‑148 (MQ=255)
    cATCGTTAAACTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTg                                                                                                                                            >  1:1486478/1‑148 (MQ=255)
        gTTAAACTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTc                                                                                                                                        <  1:884966/148‑1 (MQ=255)
         ttAAACTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCt                                                                                                                                       >  2:1501707/1‑148 (MQ=255)
            aaCTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCa                                                                                                                                    >  2:3132162/1‑148 (MQ=255)
              cTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAg                                                                                                                                  >  1:1251057/1‑148 (MQ=255)
              cTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGAGGCTTTCTGCGCAAATTTCGGCCACGCCTGGTGAGTTGCGTCTTCAAg                                                                                                                                  >  2:336546/1‑148 (MQ=255)
              cTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTAGCTGACCCAGACAAAAAATGTTTCTATATGGAACAACGCGTCTTTCTGCGCAAATTTCTGCCCCTCCGGGAGAGAAGCGTCTTCAAg                                                                                                                                  >  1:3029369/1‑148 (MQ=255)
               tCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGc                                                                                                                                 >  2:66706/1‑148 (MQ=255)
               tCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGc                                                                                                                                 >  2:2677898/1‑148 (MQ=255)
               tCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGc                                                                                                                                 >  2:2343402/1‑148 (MQ=255)
                ccTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGca                                                                                                                                >  2:2859069/1‑148 (MQ=255)
                ccTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGca                                                                                                                                >  2:1952789/1‑148 (MQ=255)
                ccTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGca                                                                                                                                >  1:1133097/1‑148 (MQ=255)
                 cTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGcac                                                                                                                               <  1:3397656/148‑1 (MQ=255)
                 cTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGcac                                                                                                                               <  2:2225155/148‑1 (MQ=255)
                  ttAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGcacc                                                                                                                              >  1:4033/1‑148 (MQ=255)
                  ttAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGcacc                                                                                                                              >  1:3090415/1‑148 (MQ=255)
                  ttAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGcacc                                                                                                                              <  2:2165465/148‑1 (MQ=255)
                     aGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACg                                                                                                                           <  2:3375266/148‑1 (MQ=255)
                       cGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGcc                                                                                                                         <  2:3189012/148‑1 (MQ=255)
                        gATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCg                                                                                                                        <  2:2495118/148‑1 (MQ=255)
                        gATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCg                                                                                                                        <  2:2495142/148‑1 (MQ=255)
                             gAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATa                                                                                                                   >  1:2110211/1‑148 (MQ=255)
                              aGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATAc                                                                                                                  >  2:102539/1‑148 (MQ=255)
                                  aaGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCa                                                                                                              >  1:1222462/1‑148 (MQ=255)
                                   aGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAg                                                                                                             >  2:1998035/1‑148 (MQ=255)
                                           ggAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGc                                                                                                     <  1:1518142/148‑1 (MQ=255)
                                               aaaGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTc                                                                                                 >  1:1357873/1‑148 (MQ=255)
                                               aaaGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTc                                                                                                 >  2:362976/1‑148 (MQ=255)
                                                aaGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTcg                                                                                                >  1:1836428/1‑148 (MQ=255)
                                                  gcgcGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTcgcg                                                                                              >  1:3344825/1‑148 (MQ=255)
                                                   cgcgTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGt                                                                                             >  2:237914/1‑148 (MQ=255)
                                                      gTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGc                                                                                          >  1:3420758/1‑148 (MQ=255)
                                                                gTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGa                                                                                <  2:2921135/148‑1 (MQ=255)
                                                                gTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGa                                                                                <  2:1301848/148‑1 (MQ=255)
                                                                 tttCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGAc                                                                               >  2:963804/1‑148 (MQ=255)
                                                                     gACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAg                                                                           <  1:2434676/148‑1 (MQ=255)
                                                                     gACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAg                                                                           <  2:1726195/148‑1 (MQ=255)
                                                                     gACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAg                                                                           <  1:2859069/148‑1 (MQ=255)
                                                                       cGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGAc                                                                         >  1:753964/1‑148 (MQ=255)
                                                                        gTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACg                                                                        <  2:2312320/148‑1 (MQ=255)
                                                                               aCCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAg                                                                 >  2:2606967/1‑148 (MQ=255)
                                                                                   aGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCCGAAAAAAGCCGTCGCGTAGCGTCGATACGAAGCACACGTTTGTAGCGAt                                                             >  1:698408/1‑148 (MQ=255)
                                                                                     ccAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGc                                                           <  2:1222462/148‑1 (MQ=255)
                                                                                        aaaaaGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCgg                                                        >  1:2701985/1‑148 (MQ=255)
                                                                                         aaaaGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGt                                                       >  1:1746879/1‑148 (MQ=255)
                                                                                         aaaaGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGt                                                       <  2:1322356/148‑1 (MQ=255)
                                                                                         aaaaGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGt                                                       >  1:1746895/1‑148 (MQ=255)
                                                                                           aaGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGa                                                     >  2:2368831/1‑148 (MQ=255)
                                                                                           aaGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGa                                                     >  1:835250/1‑148 (MQ=255)
                                                                                           aaGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGa                                                     >  1:835078/1‑148 (MQ=255)
                                                                                             ggTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGaaa                                                   >  2:2416419/1‑148 (MQ=255)
                                                                                             ggTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGaaa                                                   <  2:2962173/148‑1 (MQ=255)
                                                                                             ggTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGaaa                                                   >  1:1048620/1‑148 (MQ=255)
                                                                                             ggTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGaaa                                                   >  1:1048594/1‑148 (MQ=255)
                                                                                             ggTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGaaa                                                   <  2:1956931/148‑1 (MQ=255)
                                                                                               tttCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAAc                                                 <  1:1537388/148‑1 (MQ=255)
                                                                                               tttCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAAc                                                 <  2:1353075/148‑1 (MQ=255)
                                                                                               tttCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAAc                                                 >  1:1073330/1‑148 (MQ=255)
                                                                                                 tCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCa                                               <  2:2274991/148‑1 (MQ=255)
                                                                                                  cGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAg                                              >  1:1414614/1‑148 (MQ=255)
                                                                                                  cGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAg                                              <  1:257874/148‑1 (MQ=255)
                                                                                                  cGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAg                                              <  1:408539/148‑1 (MQ=255)
                                                                                                    atatTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGAt                                            <  1:2497657/148‑1 (MQ=255)
                                                                                                    atatTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGAt                                            >  1:1565995/1‑148 (MQ=255)
                                                                                                     tatTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATg                                           >  1:1307066/1‑148 (MQ=255)
                                                                                                       ttGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGc                                         >  2:3135251/1‑148 (MQ=255)
                                                                                                        tGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCAGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCg                                        >  1:1916096/1‑148 (MQ=255)
                                                                                                          aacaacGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAg                                                                                                      >  2:814770/1‑84 (MQ=255)
                                                                                                          aacaacGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTa                                      >  1:676845/1‑148 (MQ=255)
                                                                                                          aacaacGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTa                                      >  1:2465862/1‑148 (MQ=255)
                                                                                                           acaacGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTAt                                     <  2:3001182/148‑1 (MQ=255)
                                                                                                           acaacGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTAt                                     >  2:554544/1‑148 (MQ=255)
                                                                                                                  ggCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCgg                              <  2:3277072/148‑1 (MQ=255)
                                                                                                                  ggCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCgg                              <  2:3277024/148‑1 (MQ=255)
                                                                                                                    cTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCt                            >  2:1606986/1‑148 (MQ=255)
                                                                                                                          gcgcAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCaa                      >  2:1868940/1‑148 (MQ=255)
                                                                                                                          gcgcAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCaa                      <  1:137319/148‑1 (MQ=255)
                                                                                                                           cgcAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCaaa                     <  1:1508794/148‑1 (MQ=255)
                                                                                                                           cgcAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCaaa                     >  1:410638/1‑148 (MQ=255)
                                                                                                                           cgcAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCaaa                     >  1:2907775/1‑148 (MQ=255)
                                                                                                                           cgcAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCaaa                     >  2:1063721/1‑148 (MQ=255)
                                                                                                                           cgcAAATTTCGGCCCCGCCTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCaaa                     <  1:2211460/148‑1 (MQ=255)
                                                                                                                             cAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAAt                   >  2:3180462/1‑148 (MQ=255)
                                                                                                                                atgtCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCt                <  1:2416419/145‑1 (MQ=255)
                                                                                                                                  ttCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGt              <  2:378878/148‑1 (MQ=255)
                                                                                                                                     ggCCCCGCCTGGAGAGTTGCGTCTTCGAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCa           >  1:2201253/1‑148 (MQ=255)
                                                                                                                                      gCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAg          <  2:33268/148‑1 (MQ=255)
                                                                                                                                       ccccGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAgc         >  2:2847983/1‑148 (MQ=255)
                                                                                                                                           gCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTATCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAgcggcg     <  2:1463717/148‑1 (MQ=255)
                                                                                                                                            ccTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGt    <  2:1933650/148‑1 (MQ=255)
                                                                                                                                             cTGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTa   <  2:2748364/148‑1 (MQ=255)
                                                                                                                                              tGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTAt  >  2:552293/1‑148 (MQ=255)
                                                                                                                                              tGGAGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTAt  <  1:1596501/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTGGCATCGTTAAACTCCTTAAGCGATATGAGCAAAGGACGTGGGAGAAAGCGCGTCGCAGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTAT  >  CP009273/1067609‑1067898

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: