Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 176,783 G→T 54.5% D354Y (GAT→TAT)  dgt → deoxyguanosine triphosphate triphosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273176,7830GT54.5% 10.3 / inf 231D354Y (GAT→TAT) dgtdeoxyguanosine triphosphate triphosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (46/59);  new base T (69/57);  total (115/116)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 238 total aligned reads displayed.

AAGGTTCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGG  >  CP009273/176634‑176931
                                                                                                                                                     |                                                                                                                                                    
aaGGTTCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTg                                                                                                                                                      <  2:1815012/150‑1 (MQ=255)
aaGGTTCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTg                                                                                                                                                      >  2:1768486/1‑150 (MQ=255)
  ggTTCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATttatttat                                                                                                                                                    <  1:1327569/150‑1 (MQ=255)
    ttCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATaa                                                                                                                                                  >  2:1795187/1‑150 (MQ=255)
     tCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAAt                                                                                                                                                 >  1:3036095/1‑150 (MQ=255)
        ctctTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATctg                                                                                                                                              >  2:3712798/1‑150 (MQ=255)
            tttCGCTGGTGGTTTAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCTCAGTTCGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATctgcctg                                                                                                                                          <  1:4018116/150‑1 (MQ=255)
              tCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCg                                                                                                                                        <  2:841035/150‑1 (MQ=255)
               cGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGa                                                                                                                                       >  2:633604/1‑150 (MQ=255)
                  tggtggTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGAttt                                                                                                                                    >  2:3018313/1‑150 (MQ=255)
                  tggtggTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGAttt                                                                                                                                    >  2:1291833/1‑150 (MQ=255)
                  tggtggTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGAttt                                                                                                                                    >  1:2041857/1‑150 (MQ=255)
                   ggtggtTTAAAATGCCTGGGATAAATCACGCTAAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGAtttt                                                                                                                                   <  1:914428/150‑1 (MQ=255)
                   ggtggtTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGAtttt                                                                                                                                   <  1:2942407/150‑1 (MQ=255)
                    gtggtTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTc                                                                                                                                  <  1:3546120/150‑1 (MQ=255)
                     tggtTTTAAATGCCTGGGAAAAATTACGCTCAAATATTTTTAGCCGCAGTACGGTAGATCATTTTTTTTTGTATTTACGGGTAAACACCCTAAATAAACTGGTAACCTACGGGGCAAAAAGTTTTATTTATAATATGCCTGCGAtttttc                                                                                                                                 <  1:361480/150‑3 (MQ=255)
                       gtTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGcc                                                                                                                               <  2:3335800/150‑1 (MQ=255)
                        ttGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCg                                                                                                                              <  2:393076/150‑1 (MQ=255)
                        ttGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCg                                                                                                                              >  1:3122156/1‑150 (MQ=255)
                          gAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGa                                                                                                                            >  2:421713/1‑150 (MQ=255)
                           aaaaTGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGaa                                                                                                                           <  2:3580229/150‑1 (MQ=255)
                           aaaaTGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGaa                                                                                                                           <  2:675736/150‑1 (MQ=255)
                            aaaTGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAAc                                                                                                                          <  2:1675911/150‑1 (MQ=255)
                            aaaTGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAAc                                                                                                                          >  2:1764487/1‑150 (MQ=255)
                            aaaTGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAAc                                                                                                                          <  1:2678207/150‑1 (MQ=255)
                               tGCCTGGGAAAAATCACGCTCAAATAGTTTAGGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGtt                                                                                                                       <  1:1196552/150‑1 (MQ=255)
                               tGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGtt                                                                                                                       >  1:1334435/1‑150 (MQ=255)
                               tGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGtt                                                                                                                       <  2:909213/150‑1 (MQ=255)
                               tGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGtt                                                                                                                       <  1:3042268/150‑1 (MQ=255)
                                gCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGttt                                                                                                                      <  2:3769466/150‑1 (MQ=255)
                                  cTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTaa                                                                                                                    <  2:2025763/150‑1 (MQ=255)
                                   tGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAAt                                                                                                                   <  2:3966675/150‑1 (MQ=255)
                                     ggAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCa                                                                                                                 <  2:2766519/150‑1 (MQ=255)
                                      gAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCAt                                                                                                                <  1:3822125/150‑1 (MQ=255)
                                        aaaaTCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGc                                                                                                              >  2:3514616/1‑150 (MQ=255)
                                        aaaaTCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGc                                                                                                              >  2:3690513/1‑150 (MQ=255)
                                           aTCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCatt                                                                                                           >  2:3758768/1‑150 (MQ=255)
                                              aCGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCattatt                                                                                                        <  1:2770918/150‑1 (MQ=255)
                                                 cTCAAATAGTTTAAGCCGCAGTACGGACGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTGATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGa                                                                                                     >  2:711092/1‑150 (MQ=255)
                                                  tCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGaa                                                                                                    >  1:699775/1‑150 (MQ=255)
                                                  tCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGaa                                                                                                    <  2:3992606/150‑1 (MQ=255)
                                                     aaTAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGAt                                                                                                 <  1:1126690/150‑1 (MQ=255)
                                                      aTAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATg                                                                                                <  1:1865717/150‑1 (MQ=255)
                                                         gTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCa                                                                                             <  1:2692234/150‑1 (MQ=255)
                                                         gTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCa                                                                                             >  2:3377933/1‑150 (MQ=255)
                                                            tAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCg                                                                                          <  1:490527/150‑1 (MQ=255)
                                                             aaGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGa                                                                                         <  1:3392862/150‑1 (MQ=255)
                                                             aaGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGa                                                                                         <  1:1772606/150‑1 (MQ=255)
                                                             aaGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGa                                                                                         <  2:591294/150‑1 (MQ=255)
                                                               gCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAAt                                                                                       <  2:2307704/150‑1 (MQ=255)
                                                               gCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAAt                                                                                       <  2:3374455/150‑1 (MQ=255)
                                                                  gCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCa                                                                                    <  2:2607676/150‑1 (MQ=255)
                                                                     gTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCg                                                                                 >  1:3206607/1‑150 (MQ=255)
                                                                      tACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGa                                                                                >  1:1822543/1‑150 (MQ=255)
                                                                         ggAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGAtct                                                                             <  1:2454482/150‑1 (MQ=255)
                                                                         ggAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGAtct                                                                             >  1:2032647/1‑150 (MQ=255)
                                                                          gAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGAGTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGAtctt                                                                            >  2:1865046/1‑150 (MQ=255)
                                                                           aaGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGAtcttc                                                                           >  2:2240939/1‑150 (MQ=255)
                                                                           aaGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGAtcttc                                                                           >  1:1513888/1‑150 (MQ=255)
                                                                               tCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTaa                                                                       >  1:1183799/1‑150 (MQ=255)
                                                                                  gTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCt                                                                    >  1:2353717/1‑150 (MQ=255)
                                                                                    tttttATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCtat                                                                  <  2:1653396/150‑1 (MQ=255)
                                                                                    tttttATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCtat                                                                  <  1:2539298/150‑1 (MQ=255)
                                                                                         aTGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATaaa                                                             >  1:3184994/1‑150 (MQ=255)
                                                                                         aTGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATaaa                                                             >  2:2398727/1‑150 (MQ=255)
                                                                                          tGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATaaaa                                                            <  1:600828/150‑1 (MQ=255)
                                                                                           gTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATaaaaa                                                           >  2:469713/1‑150 (MQ=255)
                                                                                            tATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAAt                                                          >  1:1571710/1‑150 (MQ=255)
                                                                                               ttACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGAGTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTc                                                       >  1:2889787/1‑150 (MQ=255)
                                                                                                  cGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCt                                                    <  1:2662865/150‑1 (MQ=255)
                                                                                                    ggTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGt                                                  >  1:1086277/1‑150 (MQ=255)
                                                                                                     gTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATAGAAAAATGTCGCTGTa                                                 >  1:662590/1‑150 (MQ=255)
                                                                                                       aaaCACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTaaa                                               <  2:2130770/150‑1 (MQ=255)
                                                                                                        aaCACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTaaaa                                              >  1:422364/1‑150 (MQ=255)
                                                                                                        aaCACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTaaaa                                              <  2:2910951/150‑1 (MQ=255)
                                                                                                         acacCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAAc                                             <  2:143253/150‑1 (MQ=255)
                                                                                                           acCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACAt                                           <  2:1130041/150‑1 (MQ=255)
                                                                                                           acCATAAATAATCTGGTACCCTTCGCGGCACAACTATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCTTTGATTATTGGAAGATGCCAGCGATTGCAGCGATCTTCTTAAGCTATTTAAAAATGTCGCTGTAAAACAt                                           <  2:407155/150‑1 (MQ=255)
                                                                                                            cccTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACAtg                                          <  2:1326363/150‑1 (MQ=255)
                                                                                                              cTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACAtgtg                                        >  2:2624277/1‑150 (MQ=255)
                                                                                                                  ataaaCTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTa                                    <  1:2382401/150‑1 (MQ=255)
                                                                                                                  ataaaCTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACATGTGTTTa                                    <  2:2674348/150‑1 (MQ=255)
                                                                                                                   taaaCTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAg                                   <  2:1738802/150‑1 (MQ=255)
                                                                                                                   taaaCTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAg                                   <  2:3286795/150‑1 (MQ=255)
                                                                                                                    aaaCTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGc                                  <  1:3215521/150‑1 (MQ=255)
                                                                                                                      aCTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCa                                >  1:1018637/1‑150 (MQ=255)
                                                                                                                           tACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAg                           <  1:2530415/150‑1 (MQ=255)
                                                                                                                             cccTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGAt                         <  2:3501292/150‑1 (MQ=255)
                                                                                                                             cccTACGCGGCACAACGATTTAATGATAATCTGCCTGCGATTTTCGCCGGAACTTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGAt                         <  2:2462961/150‑1 (MQ=255)
                                                                                                                                tACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTc                      >  2:1246277/1‑150 (MQ=255)
                                                                                                                                tACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTc                      >  2:3339458/1‑150 (MQ=255)
                                                                                                                                 aCGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCg                     <  2:2789365/150‑1 (MQ=255)
                                                                                                                                  cgcgGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGa                    >  2:1412963/1‑150 (MQ=255)
                                                                                                                                   gcgGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAg                   <  1:3970807/150‑1 (MQ=255)
                                                                                                                                    cgGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGc                  >  2:959907/1‑150 (MQ=255)
                                                                                                                                    cgGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGc                  <  2:140483/150‑1 (MQ=255)
                                                                                                                                     ggCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCg                 >  2:3658652/1‑150 (MQ=255)
                                                                                                                                          aaCGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTg            >  1:4050211/1‑150 (MQ=255)
                                                                                                                                                 tATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCa     >  1:1183968/1‑150 (MQ=255)
                                                                                                                                                    tGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAggg  >  1:1496815/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAGGTTCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGG  >  CP009273/176634‑176931

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: