Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,017,285 G→T 89.3% D871Y (GAC→TAC)  xdhD → putative hypoxanthine oxidase, molybdopterin‑binding/Fe‑S binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,017,2850GT89.3% 823.8 / 87.7 286D871Y (GAC→TAC) xdhDputative hypoxanthine oxidase, molybdopterin‑binding/Fe‑S binding
Reads supporting (aligned to +/- strand):  ref base G (21/10);  new base T (150/105);  total (171/115)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 288 total aligned reads displayed.

CAACACGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTA  >  CP009273/3017137‑3017434
                                                                                                                                                    |                                                                                                                                                     
cAACACGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGa                                                                                                                                                      >  2:969222/1‑150 (MQ=255)
  acacGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACg                                                                                                                                                    <  1:643312/150‑1 (MQ=255)
   cacGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGc                                                                                                                                                   >  1:2842480/1‑150 (MQ=255)
   cacGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGc                                                                                                                                                   >  2:1755039/1‑150 (MQ=255)
    acGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGcc                                                                                                                                                  <  2:1516570/150‑1 (MQ=255)
      gcgTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGa                                                                                                                                                <  2:1051108/150‑1 (MQ=255)
        gTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAg                                                                                                                                              <  2:3100005/150‑1 (MQ=255)
           cGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTc                                                                                                                                           <  2:2302954/150‑1 (MQ=255)
            gggTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCa                                                                                                                                          <  1:3399049/150‑1 (MQ=255)
             ggTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCAc                                                                                                                                         <  2:4296263/150‑1 (MQ=255)
               tGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCAccc                                                                                                                                       >  1:995331/1‑150 (MQ=255)
                 aaaTCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGt                                                                                                                                     >  2:1195837/1‑150 (MQ=255)
                 aaaTCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGt                                                                                                                                     >  1:3714518/1‑150 (MQ=255)
                  aaTCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGtt                                                                                                                                    >  1:2522237/1‑150 (MQ=255)
                     ccGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAAc                                                                                                                                 <  2:267629/150‑1 (MQ=255)
                      cGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAAcg                                                                                                                                >  1:2806764/1‑150 (MQ=255)
                       gCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAAcgc                                                                                                                               <  1:4514903/150‑1 (MQ=255)
                        ccTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAAcgcg                                                                                                                              >  1:4292831/1‑150 (MQ=255)
                        ccTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAAcgcg                                                                                                                              <  2:4112021/150‑1 (MQ=255)
                         cTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGt                                                                                                                             >  1:2552647/1‑150 (MQ=255)
                            gATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGAt                                                                                                                          >  2:2773400/1‑150 (MQ=255)
                            gATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGAt                                                                                                                          >  2:3258938/1‑150 (MQ=255)
                            gATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGAt                                                                                                                          >  2:3258093/1‑150 (MQ=255)
                            gATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGAt                                                                                                                          <  1:3222063/150‑1 (MQ=255)
                              tAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGAttt                                                                                                                        >  1:904065/1‑150 (MQ=255)
                              tAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGAttt                                                                                                                        <  2:3110235/150‑1 (MQ=255)
                                     tACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGt                                                                                                                 >  1:1281542/1‑150 (MQ=255)
                                     tACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGt                                                                                                                 >  1:1430117/1‑150 (MQ=255)
                                        gCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTAc                                                                                                              <  1:844296/150‑1 (MQ=255)
                                         ccTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACg                                                                                                             >  2:2638489/1‑150 (MQ=255)
                                          cTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACgg                                                                                                            >  2:2728464/1‑150 (MQ=255)
                                           ttGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGgc                                                                                                           <  2:256229/150‑1 (MQ=255)
                                            tGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGgcg                                                                                                          <  1:2954692/150‑1 (MQ=255)
                                             gCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGgcgc                                                                                                         >  2:895812/1‑150 (MQ=255)
                                               tGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCAc                                                                                                       >  2:4138883/1‑150 (MQ=255)
                                                ggACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCAcc                                                                                                      >  1:4412300/1‑150 (MQ=255)
                                                ggACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCAcc                                                                                                      >  2:3436757/1‑150 (MQ=255)
                                                    tGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACTCGTGATTTACGCAGTTACGGCGCACCGaac                                                                                                  >  1:4201940/1‑149 (MQ=255)
                                                    tGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGaaa                                                                                                  >  2:1294156/1‑150 (MQ=255)
                                                     gCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGaaaa                                                                                                 >  1:1166348/1‑150 (MQ=255)
                                                        gTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTg                                                                                              <  1:709490/150‑1 (MQ=255)
                                                        gTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTg                                                                                              <  1:2903839/150‑1 (MQ=255)
                                                          acacCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGt                                                                                            <  1:3582151/150‑1 (MQ=255)
                                                          acacCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGt                                                                                            <  1:2358549/150‑1 (MQ=255)
                                                           cacCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTg                                                                                           >  1:2083270/1‑150 (MQ=255)
                                                           cacCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTg                                                                                           <  2:1121467/150‑1 (MQ=255)
                                                            acCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGa                                                                                          <  1:164667/150‑1 (MQ=255)
                                                            acCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGa                                                                                          <  2:3895666/150‑1 (MQ=255)
                                                              cGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACa                                                                                        <  2:466628/150‑1 (MQ=255)
                                                              cGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACa                                                                                        >  2:2318277/1‑150 (MQ=255)
                                                                gTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACAtt                                                                                      <  2:1436581/150‑1 (MQ=255)
                                                                  cAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTcc                                                                                    >  2:1879028/1‑150 (MQ=255)
                                                                   aaTCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCcg                                                                                   >  1:1535337/1‑150 (MQ=255)
                                                                       cAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTg                                                                               >  2:17390/1‑150 (MQ=255)
                                                                        agagTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGa                                                                              >  1:3781678/1‑150 (MQ=255)
                                                                            ttAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTc                                                                          <  2:1260279/150‑1 (MQ=255)
                                                                             tAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTcc                                                                         <  1:2755242/150‑1 (MQ=255)
                                                                              aGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCcg                                                                        >  1:1953946/1‑150 (MQ=255)
                                                                                 gTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCgcgc                                                                     >  2:4392688/1‑150 (MQ=255)
                                                                                          aaTCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGt                                                            <  2:1635459/150‑1 (MQ=255)
                                                                                           agcgaCGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTg                                                           <  2:2476459/146‑1 (MQ=255)
                                                                                             cTACGGTGCCACCCTGCGTGCTATCGGCCACCGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGcc                                                         >  2:3861886/1‑150 (MQ=255)
                                                                                             cTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGcc                                                         <  2:1997499/150‑1 (MQ=255)
                                                                                                cGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAg                                                      >  2:220997/1‑150 (MQ=255)
                                                                                                cGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAg                                                      >  2:221112/1‑150 (MQ=255)
                                                                                                 ggTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGc                                                     >  2:676693/1‑150 (MQ=255)
                                                                                                  gTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGcg                                                    >  1:608118/1‑150 (MQ=255)
                                                                                                        cccTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATa                                              >  1:314833/1‑150 (MQ=255)
                                                                                                         ccTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATaa                                             >  2:3732137/1‑150 (MQ=255)
                                                                                                         ccTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATaa                                             >  1:3983289/1‑150 (MQ=255)
                                                                                                          cTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGAGaaa                                            >  2:2377335/1‑150 (MQ=255)
                                                                                                             cGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTc                                         >  1:3116101/1‑150 (MQ=255)
                                                                                                              gTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCg                                        >  1:2282542/1‑150 (MQ=255)
                                                                                                              gTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCg                                        >  2:3740588/1‑150 (MQ=255)
                                                                                                               tGCTATCGGCCAAAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCgg                                       >  2:1351294/1‑150 (MQ=255)
                                                                                                                  tATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGccc                                    <  2:3630115/150‑1 (MQ=255)
                                                                                                                  tATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGccc                                    >  2:2820702/1‑150 (MQ=255)
                                                                                                                  tATCGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGccc                                    <  2:3003275/150‑1 (MQ=255)
                                                                                                                     cGGCCACAGTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGtt                                 <  2:599558/150‑1 (MQ=255)
                                                                                                                     cGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGtt                                 >  2:3564051/1‑150 (MQ=255)
                                                                                                                             gTATGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGa                         >  2:805741/1‑150 (MQ=255)
                                                                                                                               aTGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGaaa                       <  2:3249546/150‑1 (MQ=255)
                                                                                                                               aTGAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGaaa                       >  2:4050236/1‑150 (MQ=255)
                                                                                                                                 gAGCGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATc                     <  1:1195837/150‑1 (MQ=255)
                                                                                                                                 gAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATc                     <  2:330474/150‑1 (MQ=255)
                                                                                                                                    cGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAt                  <  1:2668648/150‑1 (MQ=255)
                                                                                                                                    cGAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAt                  >  2:43307/1‑150 (MQ=255)
                                                                                                                                     gAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAACTCGAtc                 >  2:3052358/1‑150 (MQ=255)
                                                                                                                                     gAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtc                 >  1:306411/1‑150 (MQ=255)
                                                                                                                                     gAAGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtc                 >  2:179663/1‑150 (MQ=255)
                                                                                                                                      aaGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtct                <  2:3794761/150‑1 (MQ=255)
                                                                                                                                      aaGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtct                >  1:1814926/1‑150 (MQ=255)
                                                                                                                                      aaGAGATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtct                >  1:4091228/1‑150 (MQ=255)
                                                                                                                                       agagATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtctc               <  2:4258355/150‑1 (MQ=255)
                                                                                                                                       agagATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtctc               >  1:2016223/1‑150 (MQ=255)
                                                                                                                                       agagATCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtctc               >  1:748096/1‑150 (MQ=255)
                                                                                                                                           aTCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGaa           >  1:2768462/1‑150 (MQ=255)
                                                                                                                                            tCATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGaaa          >  2:3083764/1‑150 (MQ=255)
                                                                                                                                             cATTTATTACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAAt         >  2:1829345/1‑150 (MQ=255)
                                                                                                                                                ttattaCGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCgg      >  2:634582/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAACACGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTA  >  CP009273/3017137‑3017434

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: