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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273176,7830GT89.6% 601.5 / 64.2 219D354Y (GAT→TAT) dgtdeoxyguanosine triphosphate triphosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (12/11);  new base T (89/107);  total (101/118)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 231 total aligned reads displayed.

GGTTCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGC  >  CP009273/176636‑176932
                                                                                                                                                   |                                                                                                                                                     
ggTTCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATttatttat                                                                                                                                                     <  1:12388/150‑1 (MQ=255)
ggTTCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATttatttat                                                                                                                                                     <  1:1681672/150‑1 (MQ=255)
   tCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAAt                                                                                                                                                  >  1:2726790/1‑150 (MQ=255)
    cGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATc                                                                                                                                                 >  2:1700381/1‑150 (MQ=255)
    cGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATc                                                                                                                                                 >  1:2071294/1‑150 (MQ=255)
         ttttCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATctgcct                                                                                                                                            >  1:2284201/1‑150 (MQ=255)
             cGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGa                                                                                                                                        >  1:2043998/1‑150 (MQ=255)
                tggtggTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGAttt                                                                                                                                     <  1:4519741/150‑1 (MQ=255)
                tggtggTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGAttt                                                                                                                                     <  1:2359576/150‑1 (MQ=255)
                  gtggtTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTc                                                                                                                                   <  2:4189612/150‑1 (MQ=255)
                          aaaTGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAAc                                                                                                                           >  1:3610801/1‑150 (MQ=255)
                          aaaTGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAAc                                                                                                                           <  2:555326/150‑1 (MQ=255)
                            aTGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGt                                                                                                                         >  2:3921321/1‑150 (MQ=255)
                            aTGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGt                                                                                                                         >  2:3235362/1‑150 (MQ=255)
                            aTGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGt                                                                                                                         >  2:842685/1‑150 (MQ=255)
                              gCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGttt                                                                                                                       <  1:218775/150‑1 (MQ=255)
                              gCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGttt                                                                                                                       >  2:3028552/1‑150 (MQ=255)
                              gCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGttt                                                                                                                       >  2:4153832/1‑150 (MQ=255)
                              gCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGttt                                                                                                                       <  1:1779124/150‑1 (MQ=255)
                              gCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGttt                                                                                                                       <  2:1902615/150‑1 (MQ=255)
                               ccTGGGAAAAATCACGCTCAAAGAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTa                                                                                                                      <  2:1638740/150‑1 (MQ=255)
                                   ggAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCa                                                                                                                  >  1:1427968/1‑150 (MQ=255)
                                    gAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCAt                                                                                                                 <  2:1589889/150‑1 (MQ=255)
                                     aaaaaTCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATg                                                                                                                <  2:424774/150‑1 (MQ=255)
                                     aaaaaTCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATg                                                                                                                <  2:469732/150‑1 (MQ=255)
                                       aaaTCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCa                                                                                                              >  1:275715/1‑150 (MQ=255)
                                       aaaTCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCa                                                                                                              <  1:593968/150‑1 (MQ=255)
                                          tCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCatta                                                                                                           <  2:3019511/150‑1 (MQ=255)
                                           cACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCattat                                                                                                          <  2:3717992/150‑1 (MQ=255)
                                              gCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTgg                                                                                                       <  2:2528834/150‑1 (MQ=255)
                                              gCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTgg                                                                                                       >  1:2661515/1‑150 (MQ=255)
                                                 cAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTTAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAg                                                                                                    <  1:3421955/150‑1 (MQ=255)
                                                  aaaTAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGa                                                                                                   >  1:3359905/1‑150 (MQ=255)
                                                  aaaTAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGa                                                                                                   <  2:3817060/150‑1 (MQ=255)
                                                   aaTAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGAt                                                                                                  <  2:2995829/150‑1 (MQ=255)
                                                      aGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGcc                                                                                               <  1:1871851/150‑1 (MQ=255)
                                                      aGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGcc                                                                                               >  2:761417/1‑150 (MQ=255)
                                                        tttAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAg                                                                                             <  1:3560919/150‑1 (MQ=255)
                                                           aaGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGa                                                                                          <  1:1227065/150‑1 (MQ=255)
                                                            aGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGaa                                                                                         >  1:1463978/1‑150 (MQ=255)
                                                            aGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGAGTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCAGTCTTGGAAGATGCCAGCGac                                                                                         >  2:3517985/1‑149 (MQ=255)
                                                             gCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAAt                                                                                        <  2:2231402/150‑1 (MQ=255)
                                                               cGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGc                                                                                      >  2:217190/1‑150 (MQ=255)
                                                                gCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCa                                                                                     <  1:1012239/150‑1 (MQ=255)
                                                                     aCGGAAGATCAGTTTTTTATGTATTTACGGNTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGAt                                                                                >  1:4279551/1‑150 (MQ=255)
                                                                     aCGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGAt                                                                                <  2:2760319/150‑1 (MQ=255)
                                                                     aCGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGAt                                                                                <  1:3698880/150‑1 (MQ=255)
                                                                         aaGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGAtcttc                                                                            <  2:4447680/150‑1 (MQ=255)
                                                                           gATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATcttctt                                                                          >  1:4109005/1‑150 (MQ=255)
                                                                           gATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATATTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATcttctt                                                                          >  1:2958402/1‑150 (MQ=255)
                                                                            aTCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTa                                                                         <  1:2145144/150‑1 (MQ=255)
                                                                            aTCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTa                                                                         >  2:2296922/1‑150 (MQ=255)
                                                                            aTCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTa                                                                         >  2:3342622/1‑150 (MQ=255)
                                                                             tCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTaa                                                                        >  2:734688/1‑150 (MQ=255)
                                                                              cAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAg                                                                       >  2:83406/1‑150 (MQ=255)
                                                                                gTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCt                                                                     >  2:4018070/1‑150 (MQ=255)
                                                                                 ttttttATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCta                                                                    <  2:3054138/150‑1 (MQ=255)
                                                                                 ttttttATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTta                                                                    >  2:2736258/1‑150 (MQ=255)
                                                                                  tttttATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCtat                                                                   >  2:1272050/1‑150 (MQ=255)
                                                                                  tttttATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATCATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCtat                                                                   <  1:2110496/150‑1 (MQ=255)
                                                                                  tttttATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCtat                                                                   <  2:2810164/150‑1 (MQ=255)
                                                                                   ttttATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCtata                                                                  <  1:969907/150‑1 (MQ=255)
                                                                                       aTGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGACTGCAGCGATCTTCTTAAGCTATATaaa                                                              >  1:3123787/1‑150 (MQ=255)
                                                                                       aTGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGGTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATaaa                                                              >  2:2009569/1‑150 (MQ=255)
                                                                                         gTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATaaaaa                                                            >  1:4582885/1‑150 (MQ=255)
                                                                                         gTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATaaaaa                                                            >  1:556934/1‑150 (MQ=255)
                                                                                         gTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATaaaaa                                                            <  2:323890/150‑1 (MQ=255)
                                                                                         gTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATaaaaa                                                            >  1:2467953/1‑150 (MQ=255)
                                                                                          tATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAAt                                                           >  1:3934826/1‑150 (MQ=255)
                                                                                             ttACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTc                                                        <  2:1111188/150‑1 (MQ=255)
                                                                                             ttACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTc                                                        <  2:486278/150‑1 (MQ=255)
                                                                                             ttACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTc                                                        >  2:4051559/1‑150 (MQ=255)
                                                                                                cGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAAAGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATATTCTTAAGCTATATAAAAATGTCGCt                                                     <  1:2355681/150‑1 (MQ=255)
                                                                                                 gggTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTg                                                    <  2:2635272/150‑1 (MQ=255)
                                                                                                      aaCACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTaaaa                                               <  2:939862/150‑1 (MQ=255)
                                                                                                        cacCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACa                                             <  1:4590839/150‑1 (MQ=255)
                                                                                                         acCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGACCGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACAt                                            >  2:3593462/1‑150 (MQ=255)
                                                                                                         acCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACAt                                            >  1:2858086/1‑150 (MQ=255)
                                                                                                             taaataaaCTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATgtgt                                        <  2:286765/150‑1 (MQ=255)
                                                                                                                  aaaCTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGc                                   >  1:1972806/1‑150 (MQ=255)
                                                                                                                    aCTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCa                                 <  2:2453023/150‑1 (MQ=255)
                                                                                                                      tGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATc                               <  1:4339143/150‑1 (MQ=255)
                                                                                                                       ggTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATcc                              >  1:328161/1‑150 (MQ=255)
                                                                                                                        gTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCa                             <  2:2633714/150‑1 (MQ=255)
                                                                                                                           cccTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGAt                          <  2:2109291/150‑1 (MQ=255)
                                                                                                                           cccTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGAt                          >  2:1008766/1‑150 (MQ=255)
                                                                                                                              tACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTc                       <  2:1147666/150‑1 (MQ=255)
                                                                                                                                cgcgGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGa                     >  1:1844056/1‑150 (MQ=255)
                                                                                                                                cgcgGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGa                     <  1:4402245/150‑1 (MQ=255)
                                                                                                                                 gcgGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAg                    <  1:2126605/150‑1 (MQ=255)
                                                                                                                                    gCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCgg                 <  1:4152494/150‑1 (MQ=255)
                                                                                                                                       caACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTTTATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCtt              <  2:1890563/150‑1 (MQ=255)
                                                                                                                                        aaCGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTg             >  2:697177/1‑150 (MQ=255)
                                                                                                                                        aaCGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTg             <  2:3824038/150‑1 (MQ=255)
                                                                                                                                        aaCGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTg             <  2:3823972/150‑1 (MQ=255)
                                                                                                                                            atttatttATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAAtt         <  1:3973848/150‑1 (MQ=255)
                                                                                                                                             tttatttaTAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTg        <  1:3015667/150‑1 (MQ=255)
                                                                                                                                              ttatttatAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGc       <  2:1955482/150‑1 (MQ=255)
                                                                                                                                                   gATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGc  <  1:4224007/150‑1 (MQ=255)
                                                                                                                                                   gATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGc  <  1:4069980/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGTTCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGC  >  CP009273/176636‑176932

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: