Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 176,864 C→T 15.5% L381L (CTA→TTA)  dgt → deoxyguanosine triphosphate triphosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273176,8640CT15.5% 816.3 / 125.2 284L381L (CTA→TTA) dgtdeoxyguanosine triphosphate triphosphohydrolase
Reads supporting (aligned to +/- strand):  ref base C (106/134);  new base T (26/18);  total (132/152)
Fisher's exact test for biased strand distribution p-value = 7.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 287 total aligned reads displayed.

GTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAA  >  CP009273/176716‑177013
                                                                                                                                                    |                                                                                                                                                     
 ttttttATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCta                                                                                                                                                     >  2:950158/1‑150 (MQ=255)
    tttATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCtatat                                                                                                                                                  >  1:1097653/1‑150 (MQ=255)
    tttATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTtatat                                                                                                                                                  <  1:1743709/150‑1 (MQ=255)
    tttATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCtatat                                                                                                                                                  >  1:3801974/1‑150 (MQ=255)
      tATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATaa                                                                                                                                                <  1:2601248/150‑1 (MQ=255)
       aTGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATaaa                                                                                                                                               >  2:3406759/1‑150 (MQ=255)
          tATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCATCGATCTTCTTAAGCTATCTaagaaa                                                                                                                                            >  2:1929418/1‑149 (MQ=255)
            tttACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGt                                                                                                                                          >  1:1138070/1‑150 (MQ=255)
              tACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCg                                                                                                                                        <  2:2437132/150‑1 (MQ=255)
               aCGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGc                                                                                                                                       >  2:2740901/1‑150 (MQ=255)
                cGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCt                                                                                                                                      >  2:668326/1‑150 (MQ=255)
                cGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCt                                                                                                                                      >  1:2441629/1‑150 (MQ=255)
                cGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCt                                                                                                                                      <  1:2821208/150‑1 (MQ=255)
                cGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCt                                                                                                                                      <  1:4107077/150‑1 (MQ=255)
                     aaaCACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTaaa                                                                                                                                 >  1:203292/1‑150 (MQ=255)
                      aaCACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTaaaa                                                                                                                                >  1:1944772/1‑150 (MQ=255)
                        cacCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACa                                                                                                                              <  1:3171345/150‑1 (MQ=255)
                        cacCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACa                                                                                                                              >  2:1805684/1‑150 (MQ=255)
                        cacCCTAAATAAACTGGTACCCTACGCGGCAAATCGAATTATTGATAATCTGCCTGAGATATTCGCCGGAACTGTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACa                                                                                                                              >  2:2403605/1‑150 (MQ=255)
                             taaataaaCTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATgtgt                                                                                                                         <  1:1261495/150‑1 (MQ=255)
                             taaataaaCTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATgtgt                                                                                                                         >  1:444352/1‑150 (MQ=255)
                              aaataaaCTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGtt                                                                                                                        >  1:2173632/1‑150 (MQ=255)
                               aataaaCTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGttt                                                                                                                       <  1:5044846/150‑1 (MQ=255)
                                ataaaCTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTa                                                                                                                      <  2:3871779/150‑1 (MQ=255)
                                 taaaCTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAg                                                                                                                     >  2:1996031/1‑150 (MQ=255)
                                  aaaCTGGTACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGc                                                                                                                    >  2:787374/1‑150 (MQ=255)
                                    aCTGGTACCCTACGCGACACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCa                                                                                                                  <  1:4608065/150‑1 (MQ=255)
                                     cTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCAt                                                                                                                 >  2:1013390/1‑150 (MQ=255)
                                      tGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATc                                                                                                                <  1:347802/150‑1 (MQ=255)
                                       ggTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATcc                                                                                                               <  2:614742/150‑1 (MQ=255)
                                       ggTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATcc                                                                                                               <  2:446432/150‑1 (MQ=255)
                                       ggTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATcc                                                                                                               >  2:1181778/1‑150 (MQ=255)
                                        gTACCCTACGCGGCACATCGATTTATTTATAATCTGCCTGCGATTTTCGCCGGATCGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTTTATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATcta                                                                                                              <  2:2348136/150‑3 (MQ=255)
                                        gTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCa                                                                                                              <  2:43635/150‑1 (MQ=255)
                                         tACCCTACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAg                                                                                                             >  1:1515874/1‑150 (MQ=255)
                                            ccTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATg                                                                                                          >  2:2709316/1‑150 (MQ=255)
                                              tACGCGGCACAACGATTTATTTATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTc                                                                                                        >  2:3528935/1‑150 (MQ=255)
                                                 gcgGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAg                                                                                                     <  2:4307097/150‑1 (MQ=255)
                                                 gcgGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAg                                                                                                     <  2:2771582/150‑1 (MQ=255)
                                                    gCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCgg                                                                                                  <  2:1651145/150‑1 (MQ=255)
                                                    gCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCgg                                                                                                  >  1:2907750/1‑150 (MQ=255)
                                                    gCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCgg                                                                                                  <  2:1412274/150‑1 (MQ=255)
                                                       caACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCtt                                                                                               <  1:2292061/150‑1 (MQ=255)
                                                         aCGATTTATTGATCATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGa                                                                                             <  1:3604736/150‑1 (MQ=255)
                                                            atttatttATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAAtt                                                                                          >  2:4093882/1‑150 (MQ=255)
                                                             tttATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTg                                                                                         >  2:324636/1‑150 (MQ=255)
                                                                    aTAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCt                                                                                  >  1:3215333/1‑150 (MQ=255)
                                                                       aTCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATc                                                                               <  1:2412224/150‑1 (MQ=255)
                                                                            cctgcGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTc                                                                          >  1:1215253/1‑150 (MQ=255)
                                                                             ctgcGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCa                                                                         <  2:882461/150‑1 (MQ=255)
                                                                              tgcGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCAt                                                                        <  2:1234432/150‑1 (MQ=255)
                                                                              tgcGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCAt                                                                        <  1:4958920/150‑1 (MQ=255)
                                                                              tgcGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCAt                                                                        >  1:4337449/1‑150 (MQ=255)
                                                                               gcGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCAtt                                                                       >  1:4830599/1‑150 (MQ=255)
                                                                                  aTTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGc                                                                    >  2:4493788/1‑150 (MQ=255)
                                                                                     ttCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGTGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGa                                                                 <  2:2933872/150‑1 (MQ=255)
                                                                                       cGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGatt                                                               <  2:2232502/150‑1 (MQ=255)
                                                                                        gCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGatta                                                              <  2:282839/150‑1 (MQ=255)
                                                                                           ggAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATTTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGAttatta                                                           <  1:778643/150‑1 (MQ=255)
                                                                                           ggAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATACAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGAttatta                                                           <  2:1781011/150‑1 (MQ=255)
                                                                                              aCGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTagag                                                        <  2:231096/150‑1 (MQ=255)
                                                                                              aCGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTagag                                                        <  2:761573/150‑1 (MQ=255)
                                                                                              aCGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTagag                                                        >  2:3385670/1‑150 (MQ=255)
                                                                                              aCGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTagag                                                        <  2:3061723/150‑1 (MQ=255)
                                                                                                 tttAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGAtt                                                     <  2:5090277/150‑1 (MQ=255)
                                                                                                 tttAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGAtt                                                     <  2:5090125/150‑1 (MQ=255)
                                                                                                 tttAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGAtt                                                     >  1:639340/1‑150 (MQ=255)
                                                                                                     aTCATGCATGATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTGAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATc                                                 <  2:2778949/150‑1 (MQ=255)
                                                                                                        aTGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTc                                              >  1:2640561/1‑150 (MQ=255)
                                                                                                        aTGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTc                                              <  2:4778653/150‑1 (MQ=255)
                                                                                                        aTGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTc                                              >  2:4546967/1‑150 (MQ=255)
                                                                                                         tGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTcc                                             <  2:3701572/150‑1 (MQ=255)
                                                                                                          gCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCt                                            >  2:1210464/1‑150 (MQ=255)
                                                                                                            attattGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCttt                                          <  1:4742043/150‑1 (MQ=255)
                                                                                                              tattGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTttat                                        <  1:2010938/150‑1 (MQ=255)
                                                                                                              tattGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTttat                                        >  2:4853054/1‑150 (MQ=255)
                                                                                                                  tgAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAg                                    <  1:538859/149‑1 (MQ=255)
                                                                                                                     aGATGCCAGCGAATGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCt                                 <  1:3277942/150‑1 (MQ=255)
                                                                                                                       aTGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGt                               <  1:1472283/150‑1 (MQ=255)
                                                                                                                       aTGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGt                               <  2:4931516/150‑1 (MQ=255)
                                                                                                                          ccAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGt                            >  2:2491536/1‑150 (MQ=255)
                                                                                                                           cAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGtt                           >  2:3506592/1‑150 (MQ=255)
                                                                                                                            aGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTa                          >  2:1057871/1‑150 (MQ=255)
                                                                                                                              cGAATGCAGCGATATTCTTAAGCTATATAAAAATGTCGCTGTATTACTTGTGTTTAGCCATCCAGATGTCGAGCGGCTTGATTTGCAGGGCTTTCGGGTCATTAGCGGATTATTAGAGATTTTTCGTCCTTTATTAAGCCTGTCGTTATc                        <  1:4298148/150‑1 (MQ=255)
                                                                                                                                 aTGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGa                     >  2:3316494/1‑150 (MQ=255)
                                                                                                                                 aTGCAGCGATCTTCTTAAGTTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGa                     <  1:3757720/150‑1 (MQ=255)
                                                                                                                                 aTGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGa                     <  1:2016955/150‑1 (MQ=255)
                                                                                                                                   gCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACt                   <  1:4901684/150‑1 (MQ=255)
                                                                                                                                       cGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTAc               >  2:3702058/1‑150 (MQ=255)
                                                                                                                                       cGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGACATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTAc               >  1:715078/1‑150 (MQ=255)
                                                                                                                                       cGATCTACTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGACTGTCGTTATCAGACTTTAc               >  1:4193417/1‑150 (MQ=255)
                                                                                                                                        gATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACt              <  1:2103036/150‑1 (MQ=255)
                                                                                                                                          tcttctTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGa            >  2:369131/1‑150 (MQ=255)
                                                                                                                                             tcttAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACt         >  1:911527/1‑150 (MQ=255)
                                                                                                                                                tAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGt      <  1:4697864/150‑1 (MQ=255)
                                                                                                                                                 aaGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTa     >  1:3165120/1‑150 (MQ=255)
                                                                                                                                                 aaGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTa     <  2:2134042/150‑1 (MQ=255)
                                                                                                                                                 aaGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTa     <  1:2055464/150‑1 (MQ=255)
                                                                                                                                                    cTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGaa  >  2:1528368/1‑150 (MQ=255)
                                                                                                                                                    cTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGaa  <  2:2186741/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGCGGCACAACGATTTATTGATAATCTGCCTGCGATTTTCGCCGGAACGTTTAATCATGCATTATTGGAAGATGCCAGCGAATGCAGCGATCTTCTTAAGCTATATAAAAATGTCGCTGTAAAACATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAA  >  CP009273/176716‑177013

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