Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 524,496 T→G 8.6% pseudogene (401/491 nt) ybbD → pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273524,4960TG 8.6% 712.5 / 56.8 245pseudogene (401/491 nt)ybbDpseudogene
Reads supporting (aligned to +/- strand):  ref base T (117/107);  new base G (9/12);  total (126/119)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 251 total aligned reads displayed.

TTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTGAGGAAAA  >  CP009273/524348‑524643
                                                                                                                                                    |                                                                                                                                                    
 ttCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGgtgt                                                                                                                                                    >  1:2973998/1‑150 (MQ=255)
  tCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGACCAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTa                                                                                                                                                   >  1:1023896/1‑150 (MQ=255)
      tAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAAcc                                                                                                                                               >  1:3865337/1‑150 (MQ=255)
      tAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAAcc                                                                                                                                               >  1:3617173/1‑150 (MQ=255)
       aaCAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑t                                                                                                                                             >  2:1163115/1‑150 (MQ=255)
        actaaTGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑tt                                                                                                                                            <  1:2726672/147‑1 (MQ=255)
        aCAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑tt                                                                                                                                            >  2:4576708/1‑150 (MQ=255)
        aCAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑tt                                                                                                                                            >  2:1974478/1‑150 (MQ=255)
         cAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑ttt                                                                                                                                           <  1:4202009/150‑1 (MQ=255)
          aaaTGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑tttt                                                                                                                                          >  2:3883342/1‑150 (MQ=255)
          aaaTGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑tttt                                                                                                                                          <  2:3080437/150‑1 (MQ=255)
          aaaTGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACC‑tttt                                                                                                                                          <  2:2578238/150‑1 (MQ=255)
              gatgatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑tttttttt                                                                                                                                      <  1:2152491/150‑1 (MQ=255)
              gatgatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑tttttttt                                                                                                                                      <  1:3374165/150‑1 (MQ=255)
                  atgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTttgttg                                                                                                                                  >  2:3907062/1‑150 (MQ=255)
                  atgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TGTTTTttgttg                                                                                                                                  >  2:2882538/1‑150 (MQ=255)
                      gggTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGa                                                                                                                              <  2:4486623/150‑1 (MQ=255)
                      gggTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGa                                                                                                                              <  2:2963042/150‑1 (MQ=255)
                           ggTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTa                                                                                                                         >  1:4533330/1‑150 (MQ=255)
                            gTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTaa                                                                                                                        <  1:1475427/150‑1 (MQ=255)
                             tttAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAAt                                                                                                                       <  1:288130/150‑1 (MQ=255)
                              ttAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATg                                                                                                                      >  1:2706162/1‑150 (MQ=255)
                               tAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGt                                                                                                                     >  1:1871976/1‑150 (MQ=255)
                                aaGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGtt                                                                                                                    <  2:1717549/150‑1 (MQ=255)
                                   aaaaCATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACg                                                                                                                 >  1:143247/1‑150 (MQ=255)
                                   aaaaCATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACg                                                                                                                 <  2:2542486/150‑1 (MQ=255)
                                   aaaaCATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACg                                                                                                                 <  1:2415865/150‑1 (MQ=255)
                                   aaaaCATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACg                                                                                                                 <  2:3301502/150‑1 (MQ=255)
                                      aCATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCat                                                                                                              >  1:1397471/1‑150 (MQ=255)
                                       cATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCatt                                                                                                             >  1:789092/1‑150 (MQ=255)
                                       cATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCatt                                                                                                             >  1:39943/1‑150 (MQ=255)
                                        aTTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCatta                                                                                                            >  2:2480567/1‑150 (MQ=255)
                                           gACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTAttt                                                                                                         >  2:2277204/1‑150 (MQ=255)
                                             cAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTta                                                                                                       <  2:1900388/150‑1 (MQ=255)
                                                 gAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAAc                                                                                                   >  1:1623129/1‑150 (MQ=255)
                                                 gAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAAc                                                                                                   >  2:4537972/1‑150 (MQ=255)
                                                   gggCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACaa                                                                                                 <  2:311767/150‑1 (MQ=255)
                                                    ggCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACaaa                                                                                                <  1:3580601/150‑1 (MQ=255)
                                                     gCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAATTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAat                                                                                               >  1:3917739/1‑150 (MQ=255)
                                                     gCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAat                                                                                               >  2:4299816/1‑150 (MQ=255)
                                                      ccAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAata                                                                                              >  1:833924/1‑150 (MQ=255)
                                                         gggAACTAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAATAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATATCAAATATTg                                                                                           <  1:539638/150‑1 (MQ=255)
                                                         gggAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTg                                                                                           >  1:3105624/1‑150 (MQ=255)
                                                          ggAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTgg                                                                                          <  2:4404305/150‑1 (MQ=255)
                                                           gAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGa                                                                                         >  1:2042233/1‑150 (MQ=255)
                                                              cAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCt                                                                                      >  1:2172332/1‑150 (MQ=255)
                                                                 gTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGa                                                                                   >  2:1687413/1‑150 (MQ=255)
                                                                  tGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGAt                                                                                  >  2:61451/1‑150 (MQ=255)
                                                                  tGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGAt                                                                                  <  1:957870/150‑1 (MQ=255)
                                                                                   aaGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCa                                                                 >  1:4291721/1‑150 (MQ=255)
                                                                                   aaGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCa                                                                 <  1:4340938/150‑1 (MQ=255)
                                                                                   aaGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCa                                                                 >  1:642664/1‑150 (MQ=255)
                                                                                   aaGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCa                                                                 >  1:3420799/1‑150 (MQ=255)
                                                                                    aGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTTTTTGGTGAAATAAATTTAACGCATTATTTTAAAAAAAATATTGGAGCTGTATTGCAGGATGATTTTGcc                                                                 >  2:4234641/1‑149 (MQ=255)
                                                                                    aGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCag                                                                <  2:923222/150‑1 (MQ=255)
                                                                                     gAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCaga                                                               >  2:759569/1‑150 (MQ=255)
                                                                                      aaTGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCagag                                                              >  1:2538689/1‑150 (MQ=255)
                                                                                         gATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTaa                                                           <  2:3018749/150‑1 (MQ=255)
                                                                                         gATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTaa                                                           <  2:539857/150‑1 (MQ=255)
                                                                                           tCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAAtt                                                         <  2:4559314/150‑1 (MQ=255)
                                                                                             cAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGc                                                       >  2:2710140/1‑150 (MQ=255)
                                                                                               gATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCATGATGCTTGTGAAGAGTAATTTCtt                                                     >  2:1268584/1‑150 (MQ=255)
                                                                                                aTACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTg                                                    <  1:1231845/150‑1 (MQ=255)
                                                                                                       aaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAag                                             >  2:2807785/1‑150 (MQ=255)
                                                                                                        aaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAaga                                            >  1:110791/1‑150 (MQ=255)
                                                                                                            tgtgGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCt                                        >  1:4234021/1‑150 (MQ=255)
                                                                                                            tgtgGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCt                                        >  1:4467904/1‑150 (MQ=255)
                                                                                                                  gtAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCa                                  >  2:4152525/1‑150 (MQ=255)
                                                                                                                   tAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCat                                 <  2:2363683/150‑1 (MQ=255)
                                                                                                                    aaaTATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCata                                <  1:3187104/150‑1 (MQ=255)
                                                                                                                    aaaTATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCata                                <  2:522168/150‑1 (MQ=255)
                                                                                                                     aaTATGGTAAGATACCAGGGCATGGGGATGGTGTTAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTTTTACAAATATTGGATCTGGATTGCCGGATGCTTGTTCAGAGTAATTGCTTGAATTAAGAGTCTATCCCatat                               <  2:1238756/150‑1 (MQ=255)
                                                                                                                     aaTATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCatat                               <  2:3748336/150‑1 (MQ=255)
                                                                                                                      atatGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCCATCCCATATc                              <  2:1412444/150‑1 (MQ=255)
                                                                                                                        atGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGa                            <  1:331782/150‑1 (MQ=255)
                                                                                                                         tGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAATAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGaa                           >  2:2212269/1‑150 (MQ=255)
                                                                                                                          ggTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAg                          >  1:3956690/1‑150 (MQ=255)
                                                                                                                           gTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAgt                         <  1:2368603/150‑1 (MQ=255)
                                                                                                                            tAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAgtc                        <  1:3709424/150‑1 (MQ=255)
                                                                                                                            tAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAgtc                        <  1:1388935/150‑1 (MQ=255)
                                                                                                                            tAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAgtc                        <  1:3276969/150‑1 (MQ=255)
                                                                                                                              aGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAgtcgt                      >  1:2972628/1‑150 (MQ=255)
                                                                                                                                 tACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCaa                   >  2:1042326/1‑150 (MQ=255)
                                                                                                                                   ccAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACt                 >  1:2445737/1‑150 (MQ=255)
                                                                                                                                    cAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACtt                >  1:1143063/1‑150 (MQ=255)
                                                                                                                                     aGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTc               <  2:2433531/150‑1 (MQ=255)
                                                                                                                                     aGGGCAAGGGGATGGGGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTc               >  2:3330794/1‑150 (MQ=255)
                                                                                                                                       ggCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGt             <  2:143247/150‑1 (MQ=255)
                                                                                                                                       ggCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGt             <  2:2306989/150‑1 (MQ=255)
                                                                                                                                        gCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTa            <  1:2279343/150‑1 (MQ=255)
                                                                                                                                         cAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAg           >  1:2503067/1‑150 (MQ=255)
                                                                                                                                         cAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAg           <  1:729524/150‑1 (MQ=255)
                                                                                                                                           aGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAATAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTg         <  1:1000216/150‑1 (MQ=255)
                                                                                                                                           aGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTg         <  2:3634481/150‑1 (MQ=255)
                                                                                                                                           aGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTg         >  2:1354761/1‑150 (MQ=255)
                                                                                                                                           aGGGGATGGGGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTg         <  2:1247147/150‑1 (MQ=255)
                                                                                                                                            ggggATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTGa        <  1:2788422/150‑1 (MQ=255)
                                                                                                                                              ggATGGGGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTGTgg      >  1:4111424/1‑150 (MQ=255)
                                                                                                                                               gATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTGAGGa     >  2:2172779/1‑150 (MQ=255)
                                                                                                                                                  ggTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTGAGGaaaa  >  2:1081269/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTGAGGAAAA  >  CP009273/524348‑524643

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: ‑