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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273176,7830GT98.1% 981.2 / 7.4 310D354Y (GAT→TAT) dgtdeoxyguanosine triphosphate triphosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base T (143/161);  total (147/163)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 321 total aligned reads displayed.

GGTTCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAG  >  CP009273/176636‑176929
                                                                                                                                                   |                                                                                                                                                  
ggTTCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGACGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGAt                                                                                                                                                  >  2:3667736/1‑150 (MQ=255)
   tCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAAt                                                                                                                                               >  1:2667670/1‑150 (MQ=255)
        ctTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATctgcc                                                                                                                                          >  1:1422000/1‑150 (MQ=255)
         ttttCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATctgcct                                                                                                                                         <  1:684988/150‑1 (MQ=255)
         ttttCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATctgcct                                                                                                                                         <  2:2806823/150‑1 (MQ=255)
          tttCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATctgcctg                                                                                                                                        <  2:4387245/150‑1 (MQ=255)
             cGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGa                                                                                                                                     >  2:4700847/1‑150 (MQ=255)
                tggtggTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGAttt                                                                                                                                  >  2:1228328/1‑150 (MQ=255)
                 ggtggtTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGAtttt                                                                                                                                 <  1:2347630/150‑1 (MQ=255)
                   tggtTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCg                                                                                                                               >  2:355533/1‑150 (MQ=255)
                     gtTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGcc                                                                                                                             <  1:2524739/150‑1 (MQ=255)
                         aaaaTGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGaa                                                                                                                         >  1:3735219/1‑150 (MQ=255)
                             tGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGtt                                                                                                                     <  1:2821780/150‑1 (MQ=255)
                                cTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTaa                                                                                                                  >  1:981750/1‑150 (MQ=255)
                                 tGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAAt                                                                                                                 >  2:5083853/1‑150 (MQ=255)
                                 tGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAAAACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAAt                                                                                                                 <  2:2172260/150‑1 (MQ=255)
                                  gggAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATc                                                                                                                >  1:2754800/1‑150 (MQ=255)
                                   ggAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCa                                                                                                               >  2:4658657/1‑150 (MQ=255)
                                     aaaaaTCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATg                                                                                                             >  2:416981/1‑150 (MQ=255)
                                       aaaTCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCa                                                                                                           >  1:3363848/1‑150 (MQ=255)
                                       aaaTCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCa                                                                                                           <  2:4505082/150‑1 (MQ=255)
                                       aaaTCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCa                                                                                                           >  2:5270070/1‑150 (MQ=255)
                                         aTCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCatt                                                                                                         >  1:1755800/1‑150 (MQ=255)
                                           cACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGACTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCattat                                                                                                       <  1:3472622/150‑1 (MQ=255)
                                            aCGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCattatt                                                                                                      <  1:608163/150‑1 (MQ=255)
                                            aCGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTGTAATCATGCattatt                                                                                                      >  2:2878725/1‑150 (MQ=255)
                                              gCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTgg                                                                                                    <  1:2461149/150‑1 (MQ=255)
                                              gCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTgg                                                                                                    >  2:3548491/1‑150 (MQ=255)
                                               cTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGa                                                                                                   <  2:1586772/150‑1 (MQ=255)
                                               cTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGa                                                                                                   <  1:837995/150‑1 (MQ=255)
                                                tCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGaa                                                                                                  >  1:1678413/1‑150 (MQ=255)
                                                 cAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAg                                                                                                 >  1:4201739/1‑150 (MQ=255)
                                                   aaTAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGAt                                                                                               >  1:5246779/1‑150 (MQ=255)
                                                    aTAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATg                                                                                              <  2:621168/150‑1 (MQ=255)
                                                         ttAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGc                                                                                         >  2:3624377/1‑150 (MQ=255)
                                                           aaGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGa                                                                                       >  1:3831202/1‑150 (MQ=255)
                                                              ccGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATg                                                                                    <  2:888960/150‑1 (MQ=255)
                                                               cGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGc                                                                                   <  1:3576953/150‑1 (MQ=255)
                                                               cGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGc                                                                                   <  1:5004426/150‑1 (MQ=255)
                                                                 cAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAg                                                                                 <  2:1037483/150‑1 (MQ=255)
                                                                  aGTACGTAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGAATTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGc                                                                                <  2:796552/150‑1 (MQ=255)
                                                                  aGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGc                                                                                <  1:965303/150‑1 (MQ=255)
                                                                      cGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACTATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGAtc                                                                            <  1:2934261/150‑1 (MQ=255)
                                                                      cGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGAtc                                                                            <  1:1159661/150‑1 (MQ=255)
                                                                      cGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGAtc                                                                            >  2:4759768/1‑150 (MQ=255)
                                                                       ggAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGAtct                                                                           <  2:3925440/150‑1 (MQ=255)
                                                                       ggAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGACACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGAtct                                                                           <  2:5124903/150‑1 (MQ=255)
                                                                        gAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGAtctt                                                                          <  1:4761112/150‑1 (MQ=255)
                                                                        gAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGAtctt                                                                          <  2:1522980/150‑1 (MQ=255)
                                                                           gATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATcttctt                                                                       >  1:4860869/1‑150 (MQ=255)
                                                                            aTCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTa                                                                      >  2:324086/1‑150 (MQ=255)
                                                                             tCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTaa                                                                     >  1:1272641/1‑150 (MQ=255)
                                                                              cAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAg                                                                    >  1:394699/1‑150 (MQ=255)
                                                                              cAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAg                                                                    <  2:1678413/150‑1 (MQ=255)
                                                                              cAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAg                                                                    <  2:845763/150‑1 (MQ=255)
                                                                                gTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCt                                                                  >  1:1949749/1‑150 (MQ=255)
                                                                                 ttttttATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCta                                                                 >  2:3338460/1‑150 (MQ=255)
                                                                                        tGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATaaaa                                                          >  1:3384368/1‑150 (MQ=255)
                                                                                        tGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATaaaa                                                          >  1:445345/1‑150 (MQ=255)
                                                                                          tATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAAt                                                        >  1:3504546/1‑150 (MQ=255)
                                                                                           aTTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATg                                                       >  2:4094198/1‑150 (MQ=255)
                                                                                           aTTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATg                                                       >  1:366104/1‑150 (MQ=255)
                                                                                           aTTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATg                                                       <  1:2067762/150‑1 (MQ=255)
                                                                                           aTTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATg                                                       >  2:484589/1‑150 (MQ=255)
                                                                                                cGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCt                                                  >  1:3126759/1‑150 (MQ=255)
                                                                                                 gggTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTg                                                 <  2:784947/150‑1 (MQ=255)
                                                                                                  ggTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGt                                                >  2:2053715/1‑150 (MQ=255)
                                                                                                    tAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTaa                                              >  2:431012/1‑150 (MQ=255)
                                                                                                     aaaCACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTaaa                                             <  2:1308801/150‑1 (MQ=255)
                                                                                                     aaaCACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTaaa                                             >  1:3974629/1‑150 (MQ=255)
                                                                                                       acacCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAAc                                           <  1:2471289/150‑1 (MQ=255)
                                                                                                          cccTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACAtg                                        >  2:4111533/1‑150 (MQ=255)
                                                                                                           ccTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACAtgt                                       <  1:881646/150‑1 (MQ=255)
                                                                                                            cTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACAtgtg                                      <  2:424601/150‑1 (MQ=255)
                                                                                                              aaataaaCTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGtt                                    >  1:5281978/1‑150 (MQ=255)
                                                                                                               aataaaCTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGttt                                   <  2:3189709/150‑1 (MQ=255)
                                                                                                               aataaaCTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGttt                                   <  1:5155223/150‑1 (MQ=255)
                                                                                                               aataaaCTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGttt                                   >  1:3435498/1‑150 (MQ=255)
                                                                                                                   aaCTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGcc                               >  1:4309372/1‑150 (MQ=255)
                                                                                                                     cTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCAt                             <  2:2514097/150‑1 (MQ=255)
                                                                                                                     cTGGTACCCTACGCGGCACAACGATTTATTTATAAGCTGCCTGCGATTTTCGCCGGAACGTATAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCAt                             <  2:3378463/150‑1 (MQ=255)
                                                                                                                        gTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCa                          <  2:187913/150‑1 (MQ=255)
                                                                                                                           cccTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGAt                       <  2:3831202/150‑1 (MQ=255)
                                                                                                                           cccTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGAt                       >  1:3008040/1‑150 (MQ=255)
                                                                                                                           cccTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGAt                       >  1:4610893/1‑150 (MQ=255)
                                                                                                                           cccTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGAt                       <  2:3011473/150‑1 (MQ=255)
                                                                                                                           cccTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGAt                       <  2:290755/150‑1 (MQ=255)
                                                                                                                           cccTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGAt                       <  2:1919455/150‑1 (MQ=255)
                                                                                                                            ccTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATg                      <  1:437476/150‑1 (MQ=255)
                                                                                                                                 gcgGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAg                 <  1:2587384/150‑1 (MQ=255)
                                                                                                                                  cgGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGc                <  2:4612254/150‑1 (MQ=255)
                                                                                                                                   ggCACACCGAGTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCg               <  2:3685514/150‑1 (MQ=255)
                                                                                                                                        aaCGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTg          >  2:2734307/1‑150 (MQ=255)
                                                                                                                                            atttatttATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAAtt      >  1:3731409/1‑150 (MQ=255)
                                                                                                                                            atttatttATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAAtt      >  2:2330137/1‑150 (MQ=255)
                                                                                                                                             tttatttaTAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTg     <  2:4011442/150‑1 (MQ=255)
                                                                                                                                               tATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCa   <  2:4948948/150‑1 (MQ=255)
                                                                                                                                               tATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCa   <  2:1401481/150‑1 (MQ=255)
                                                                                                                                                aTTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAg  <  2:1422000/150‑1 (MQ=255)
                                                                                                                                                aTTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAg  <  1:5058644/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGTTCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAG  >  CP009273/176636‑176929

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: