Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,995,335 C→A 10.5% V544V (GTC→GTA xdhA → xanthine dehydrogenase, molybdenum binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,995,3350CA10.5% 601.6 / 32.5 181V544V (GTC→GTAxdhAxanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base C (81/81);  new base A (7/12);  total (88/93)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 184 total aligned reads displayed.

CGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACGCGGCGGGCAGCTCTCTGCTGAA  >  CP009273/2995191‑2995482
                                                                                                                                                |                                                                                                                                                   
cGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATg                                                                                                                                                  <  2:745514/148‑1 (MQ=255)
  ttCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCt                                                                                                                                                >  2:1996617/1‑148 (MQ=255)
  ttCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCt                                                                                                                                                >  2:1996597/1‑148 (MQ=255)
      cgTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACAt                                                                                                                                            <  2:911775/148‑1 (MQ=255)
       gTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACAtc                                                                                                                                           <  2:1402422/148‑1 (MQ=255)
        ttATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACAtca                                                                                                                                          >  1:266345/1‑148 (MQ=255)
         tATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACAtcag                                                                                                                                         >  2:313461/1‑148 (MQ=255)
          aTTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACAtcagt                                                                                                                                        <  1:2417637/148‑1 (MQ=255)
          aTTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTAATGCTACAtcagt                                                                                                                                        <  2:3049685/148‑1 (MQ=255)
           tttCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACAtcagtc                                                                                                                                       >  1:1607104/1‑148 (MQ=255)
           tttCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACGATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACAtcagtc                                                                                                                                       >  1:1606820/1‑148 (MQ=255)
            ttCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACAtcagtca                                                                                                                                      >  1:897100/1‑148 (MQ=255)
              cAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGc                                                                                                                                    >  2:547214/1‑148 (MQ=255)
                aCTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTAATGCTACATCAGTCAGCGa                                                                                                                                  <  1:721843/148‑1 (MQ=255)
                   cAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGa                                                                                                                               <  2:647583/148‑1 (MQ=255)
                   cAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGa                                                                                                                               >  2:789506/1‑148 (MQ=255)
                    aaGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGaa                                                                                                                              <  1:212677/148‑1 (MQ=255)
                     agagACCGACGTTACGCCGTTCGATCCCGGCGCATGTGCCTCAAGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAACTCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAAt                                                                                                                             <  2:2786189/148‑1 (MQ=255)
                        tACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTg                                                                                                                          >  1:1902791/1‑148 (MQ=255)
                        tACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTg                                                                                                                          <  1:3141225/148‑1 (MQ=255)
                                   cGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATaaa                                                                                                               <  1:2204656/148‑1 (MQ=255)
                                     ccGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAg                                                                                                             >  2:3255322/1‑148 (MQ=255)
                                        ttCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGcc                                                                                                          <  1:2896656/148‑1 (MQ=255)
                                         tCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCa                                                                                                         >  2:1743356/1‑148 (MQ=255)
                                         tCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCa                                                                                                         >  2:1744060/1‑148 (MQ=255)
                                            aTCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCatat                                                                                                      <  2:414781/148‑1 (MQ=255)
                                             tCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATc                                                                                                     >  1:365758/1‑148 (MQ=255)
                                             tCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATc                                                                                                     >  1:1854382/1‑148 (MQ=255)
                                             tCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTAATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATc                                                                                                     >  1:1246383/1‑148 (MQ=255)
                                              cccGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCg                                                                                                    >  2:59489/1‑148 (MQ=255)
                                                cGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTg                                                                                                  <  2:492772/148‑1 (MQ=255)
                                                cGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTg                                                                                                  >  1:310177/1‑148 (MQ=255)
                                                   cgcATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTg                                                                                               <  2:344960/148‑1 (MQ=255)
                                                   cgcATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTAATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTg                                                                                               >  1:444584/1‑148 (MQ=255)
                                                    gcATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTgg                                                                                              <  1:1313756/148‑1 (MQ=255)
                                                     cATTTGCCTCACGCCTGATCTCTGTTGCCGCGCCTGCGCGGCGGAGTGCGGCTCTATTATTATTAGAGTAACTCATGGCTCACGCCGCTGTCATGCTACATCCGTCAGCGATGAATCTGACACTGATAAAAGGCCATATCGTGCTGGt                                                                                             <  2:3128848/148‑1 (MQ=255)
                                                     cATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGt                                                                                             <  2:1487647/148‑1 (MQ=255)
                                                     cATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGt                                                                                             <  2:1136817/148‑1 (MQ=255)
                                                       tttGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTg                                                                                           >  2:2039916/1‑148 (MQ=255)
                                                             tCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATAAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGAc                                                                                     >  2:1016184/1‑148 (MQ=255)
                                                              cACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGAcc                                                                                    >  1:56371/1‑148 (MQ=255)
                                                                 gCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGa                                                                                 <  1:892828/148‑1 (MQ=255)
                                                                     gagCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAagag                                                                             >  1:637690/1‑148 (MQ=255)
                                                                     gagCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAagag                                                                             >  1:637712/1‑148 (MQ=255)
                                                                       gCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGcc                                                                           >  1:3080274/1‑148 (MQ=255)
                                                                       gCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGcc                                                                           <  1:2277281/148‑1 (MQ=255)
                                                                            gTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTaa                                                                      <  2:54531/148‑1 (MQ=255)
                                                                               gCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGt                                                                   <  2:3022455/148‑1 (MQ=255)
                                                                                 cgcgCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCg                                                                 <  2:903486/148‑1 (MQ=255)
                                                                                  gcgcCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGt                                                                <  1:2017817/148‑1 (MQ=255)
                                                                                      ctgcgctgcgCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTaaa                                                            >  1:1535569/1‑148 (MQ=255)
                                                                                      ctgcgctgcgCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTaaa                                                            >  1:2981660/1‑148 (MQ=255)
                                                                                      ctgcgctgcgCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTAATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTaaa                                                            >  2:3060190/1‑148 (MQ=255)
                                                                                       tgcgctgcgcAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTaaaa                                                           <  2:2832037/148‑1 (MQ=255)
                                                                                       tgcgctgcgcAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTaaaa                                                           >  2:3163378/1‑148 (MQ=255)
                                                                                       tgcgctgcgcAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTaaaa                                                           >  1:1584636/1‑148 (MQ=255)
                                                                                              cgcAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTAATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTgg                                                    <  1:1657612/148‑1 (MQ=255)
                                                                                                cAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCg                                                  <  1:1743356/148‑1 (MQ=255)
                                                                                                cAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCg                                                  <  1:1744060/148‑1 (MQ=255)
                                                                                                 aGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGa                                                 <  1:2054082/148‑1 (MQ=255)
                                                                                                  gTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGAt                                                <  2:74000/148‑1 (MQ=255)
                                                                                                   tGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTTGCGATg                                               >  1:1550133/1‑148 (MQ=255)
                                                                                                    gCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATgg                                              <  2:45876/148‑1 (MQ=255)
                                                                                                    gCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATgg                                              <  2:501267/148‑1 (MQ=255)
                                                                                                    gCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATgg                                              <  2:2365731/148‑1 (MQ=255)
                                                                                                       gCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACg                                           >  2:1580731/1‑148 (MQ=255)
                                                                                                        cACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCTATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGc                                          <  1:1020757/148‑1 (MQ=255)
                                                                                                        cACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGc                                          <  1:1996617/148‑1 (MQ=255)
                                                                                                        cACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGc                                          <  1:1944792/148‑1 (MQ=255)
                                                                                                        cACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGc                                          >  2:844810/1‑148 (MQ=255)
                                                                                                        cACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGc                                          <  1:323264/148‑1 (MQ=255)
                                                                                                        cACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGc                                          <  1:2666089/148‑1 (MQ=255)
                                                                                                        cACTAGTATTAACAGAGAAAATCCTCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGc                                          <  1:1020665/148‑1 (MQ=255)
                                                                                                         aCTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCt                                         >  2:1061583/1‑148 (MQ=255)
                                                                                                          cTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTAATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTAGTTAAAAGATTTGGCGATGGACGCtt                                        >  1:3285177/1‑148 (MQ=255)
                                                                                                           tattatTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCttt                                       >  2:1390356/1‑148 (MQ=255)
                                                                                                            attattAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTc                                      >  1:2716440/1‑148 (MQ=255)
                                                                                                              tattaAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTa                                    >  1:3250245/1‑148 (MQ=255)
                                                                                                                ttaAAAGAGAAAATCATCGCTCACGCCGCAGTAATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTacc                                  >  2:1032500/1‑148 (MQ=255)
                                                                                                                  aaaaGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTaccac                                >  2:716822/1‑148 (MQ=255)
                                                                                                                     agagAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCt                             >  1:1059875/1‑148 (MQ=255)
                                                                                                                      gagaAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTTGCGATGGTCGCTTTCTACCACCCTg                            >  1:1141661/1‑148 (MQ=255)
                                                                                                                       agaAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGa                           >  1:2328249/1‑148 (MQ=255)
                                                                                                                       agaAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGa                           >  2:2061765/1‑148 (MQ=255)
                                                                                                                       agaAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGa                           >  2:401938/1‑148 (MQ=255)
                                                                                                                       agaAAATCATCGCTCACGCCGCAGTAATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGa                           <  2:2255197/148‑1 (MQ=255)
                                                                                                                        gaAAATCATCGCTCACGCCGCAGTCATGCTACAACAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGaa                          <  1:938923/148‑1 (MQ=255)
                                                                                                                          aaaTCATCGCTCACGCCGCAGTAATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAAcg                        <  1:1400010/148‑1 (MQ=255)
                                                                                                                             tcatcGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACgcgg                     >  2:2949332/1‑148 (MQ=255)
                                                                                                                              catcGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACgcggc                    >  2:2605277/1‑148 (MQ=255)
                                                                                                                                tcGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACGcggcgg                  >  2:1496632/1‑148 (MQ=255)
                                                                                                                                 cGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACGCGGCggg                 <  1:716822/148‑1 (MQ=255)
                                                                                                                                    tCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACGCGGCGGGCAg              >  2:907262/1‑148 (MQ=255)
                                                                                                                                      aCGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCAGAACGCGGCGGGCAGct            >  1:1049984/1‑148 (MQ=255)
                                                                                                                                       cgccgcAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAATGGCCATATCGTGCTGGTTGAACTACCGGAAGAGCCGTTTATTTCGTTAAAATATTTTGCGAAGGACGCTTTATACCACACTGAAAGAGGCGGGCAGatc           >  2:616686/1‑148 (MQ=255)
                                                                                                                                       cgccgcAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACGCGGCGGGCAGctc           >  1:1321108/1‑148 (MQ=255)
                                                                                                                                         ccgcAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACGCGGCGGGCAGctctc         <  2:2981660/148‑1 (MQ=255)
                                                                                                                                               tCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTTGCGCTTGACGCGTTGTACCACCCTTAACGCGGCGGGCAGCTCTctggtg  >  2:1883827/1‑147 (MQ=255)
                                                                                                                                                cATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACGCGGCGGGCCGCTCTCTGCTGaa  >  1:2825843/1‑148 (MQ=255)
                                                                                                                                                cATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACGCGGCGGGCAGCTCTCTGCTGaa  >  2:2689742/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACGCGGCGGGCAGCTCTCTGCTGAA  >  CP009273/2995191‑2995482

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