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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273524,4960TG97.0% 726.9 / 3.7 197pseudogene (401/491 nt)ybbDpseudogene
Reads supporting (aligned to +/- strand):  ref base T (4/2);  new base G (104/86);  total (109/88)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 201 total aligned reads displayed.

TTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTGAGGAAA  >  CP009273/524349‑524642
                                                                                                                                                    |                                                                                                                                                   
ttCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGgggt                                                                                                                                                   <  1:4556206/150‑1 (MQ=255)
ttCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGgggt                                                                                                                                                   <  1:1155453/150‑1 (MQ=255)
ttCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGgggt                                                                                                                                                   >  1:1874528/1‑150 (MQ=255)
 tCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTa                                                                                                                                                  <  2:1181810/150‑1 (MQ=255)
  cccTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTaa                                                                                                                                                 >  1:4088854/1‑150 (MQ=255)
  cccTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTaa                                                                                                                                                 >  2:255455/1‑150 (MQ=255)
     tAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAAcc                                                                                                                                              <  2:3012668/150‑1 (MQ=255)
       aCAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCtt                                                                                                                                           >  2:4738083/1‑150 (MQ=255)
       aCAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACAtt                                                                                                                                           >  1:3200651/1‑150 (MQ=255)
        cAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCttt                                                                                                                                          <  2:406179/150‑1 (MQ=255)
        cAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCttt                                                                                                                                          >  1:52832/1‑150 (MQ=255)
        cAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCttt                                                                                                                                          <  2:468583/150‑1 (MQ=255)
             gatgatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCtttttttt                                                                                                                                     >  1:3048600/1‑150 (MQ=255)
               tgatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTttgt                                                                                                                                   >  2:3662964/1‑150 (MQ=255)
                 atgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTttgttg                                                                                                                                 >  2:859632/1‑150 (MQ=255)
                    ggggTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAGAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTg                                                                                                                              <  2:1156187/150‑1 (MQ=255)
                      ggTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGaa                                                                                                                            <  2:3562723/150‑1 (MQ=255)
                      ggTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGaa                                                                                                                            >  1:2400459/1‑150 (MQ=255)
                         aGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAAtt                                                                                                                         <  2:3255598/150‑1 (MQ=255)
                         aGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAAtt                                                                                                                         <  1:2916741/150‑1 (MQ=255)
                           gTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTaa                                                                                                                       <  1:1377723/150‑1 (MQ=255)
                            tttAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAAt                                                                                                                      <  1:3746598/150‑1 (MQ=255)
                             tnAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATg                                                                                                                     >  1:1931417/1‑150 (MQ=255)
                               aaGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGtt                                                                                                                   >  1:79223/1‑150 (MQ=255)
                                aGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTa                                                                                                                  >  2:4015608/1‑150 (MQ=255)
                                aGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTa                                                                                                                  >  2:1437677/1‑150 (MQ=255)
                                   aaaCATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGc                                                                                                               <  2:200484/150‑1 (MQ=255)
                                    aaCATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCa                                                                                                              >  2:253185/1‑150 (MQ=255)
                                        ttGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCattat                                                                                                          >  1:1011627/1‑150 (MQ=255)
                                          gACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTAttt                                                                                                        >  1:4246056/1‑150 (MQ=255)
                                            cAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTta                                                                                                      >  2:2364869/1‑150 (MQ=255)
                                                gAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAAc                                                                                                  <  2:3418390/150‑1 (MQ=255)
                                                      cAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAatat                                                                                            >  1:1781019/1‑150 (MQ=255)
                                                      cAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAatat                                                                                            >  1:4100459/1‑150 (MQ=255)
                                                       aGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATAtt                                                                                           >  2:4571141/1‑150 (MQ=255)
                                                          gAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGa                                                                                        >  2:4081732/1‑150 (MQ=255)
                                                          gAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAg                                                                                       >  2:276885/1‑150 (MQ=255)
                                                            aCAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGc                                                                                      >  2:1902997/1‑150 (MQ=255)
                                                              aaGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTg                                                                                    <  1:256604/150‑1 (MQ=255)
                                                                gTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGa                                                                                  >  1:659762/1‑150 (MQ=255)
                                                                 tGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGAt                                                                                 <  2:1516078/150‑1 (MQ=255)
                                                                 tGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGAt                                                                                 <  2:2133113/150‑1 (MQ=255)
                                                                   gCGAAAATCGTATCAAAGAATGATACAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTATTAGTTGGTGAAATTAATGTTACGCATAATGTTATAACAAATATATGAGCTGGATTg                                                                               >  1:1550054/1‑150 (MQ=255)
                                                                      aaaaTCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCt                                                                            >  1:2066539/1‑150 (MQ=255)
                                                                       aaaTCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTg                                                                           >  2:382833/1‑150 (MQ=255)
                                                                       aaaTCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTg                                                                           <  2:1954837/150‑1 (MQ=255)
                                                                       aaaTCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTg                                                                           <  1:859632/150‑1 (MQ=255)
                                                                        aaTCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTg                                                                           >  1:1344405/1‑150 (MQ=255)
                                                                           cGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGc                                                                       <  2:1954986/150‑1 (MQ=255)
                                                                            gTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCt                                                                      <  1:1429618/150‑1 (MQ=255)
                                                                             tATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCtt                                                                     <  2:2850982/150‑1 (MQ=255)
                                                                              aTCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTtg                                                                    <  1:158304/150‑1 (MQ=255)
                                                                                   aGAATGATCCAGATACAAAAAAAGTGTGGTGTATATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCag                                                               <  2:4325868/150‑1 (MQ=255)
                                                                                   aGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCag                                                               >  2:4688635/1‑150 (MQ=255)
                                                                                   aGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCag                                                               >  1:107029/1‑150 (MQ=255)
                                                                                     aaTGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCagag                                                             >  1:1085974/1‑150 (MQ=255)
                                                                                     aaTGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCagag                                                             >  1:4041571/1‑150 (MQ=255)
                                                                                      aTGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGt                                                            >  1:3903648/1‑150 (MQ=255)
                                                                                       tGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAAATAATGTTACGCATTATTTTATAACAAATTTTGGAGCGGGATTGGCTGATGCTTGTGCAGAGTa                                                           >  2:4877380/1‑150 (MQ=255)
                                                                                         aTCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAAt                                                         >  1:887043/1‑150 (MQ=255)
                                                                                         aTCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAAt                                                         >  1:3567311/1‑150 (MQ=255)
                                                                                          tCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGAGGGGGTAACCCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAAtt                                                        <  2:1542219/150‑1 (MQ=255)
                                                                                           ccAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTg                                                       <  2:4173092/150‑1 (MQ=255)
                                                                                             aGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGca                                                     <  2:4434710/150‑2 (MQ=255)
                                                                                             aGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCt                                                     >  2:455240/1‑150 (MQ=255)
                                                                                                tACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGa                                                  <  1:961212/150‑1 (MQ=255)
                                                                                                tACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGa                                                  >  1:1704050/1‑150 (MQ=255)
                                                                                                   aaaaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAAtt                                               >  2:3719284/1‑150 (MQ=255)
                                                                                                    aaaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGTGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTa                                              >  1:1712076/1‑150 (MQ=255)
                                                                                                      caaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTCATTGCTTGAATTAag                                            >  2:1413196/2‑150 (MQ=255)
                                                                                                      aaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAag                                            >  1:1242901/1‑150 (MQ=255)
                                                                                                         aGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGt                                         >  2:1816476/1‑150 (MQ=255)
                                                                                                             tggtgTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTAt                                     <  1:3793583/150‑1 (MQ=255)
                                                                                                             tggtgTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTAt                                     <  1:3402347/150‑1 (MQ=255)
                                                                                                               gtgtAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATcc                                   <  2:4623759/150‑1 (MQ=255)
                                                                                                               gtgtAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATcc                                   >  2:3180396/1‑150 (MQ=255)
                                                                                                                 gtAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCa                                 <  1:5088174/150‑1 (MQ=255)
                                                                                                                 gtAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCa                                 <  1:3652328/150‑1 (MQ=255)
                                                                                                                   aaaTATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCata                               >  2:4542150/1‑150 (MQ=255)
                                                                                                                   aaaTATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCata                               >  1:4413584/1‑150 (MQ=255)
                                                                                                                   aaaTATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCata                               >  1:2286280/1‑150 (MQ=255)
                                                                                                                    aaTATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCatat                              >  2:4482926/1‑150 (MQ=255)
                                                                                                                    aaTATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCatat                              <  2:447006/150‑1 (MQ=255)
                                                                                                                      tatGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCataaaa                            >  2:114421/1‑147 (MQ=255)
                                                                                                                      tatGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCg                            <  2:695358/150‑1 (MQ=255)
                                                                                                                           tAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAgt                        >  1:173985/1‑150 (MQ=255)
                                                                                                                            aaGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAgtcg                      <  2:3602865/150‑1 (MQ=255)
                                                                                                                                tACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCaa                  <  2:2908369/150‑1 (MQ=255)
                                                                                                                                   cAGGGCAATGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACt                >  2:1265412/1‑150 (MQ=255)
                                                                                                                                    aGGTCAAGGTGATGGGGTAAACATTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTc              <  2:5067427/150‑1 (MQ=255)
                                                                                                                                        cAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAg          <  2:5116790/150‑1 (MQ=255)
                                                                                                                                        cAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAg          >  1:4443352/1‑150 (MQ=255)
                                                                                                                                        cAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAg          <  2:322669/150‑1 (MQ=255)
                                                                                                                                         aaGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGt         >  1:3221604/1‑150 (MQ=255)
                                                                                                                                          aGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTg        <  1:2688810/150‑1 (MQ=255)
                                                                                                                                          aGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTg        <  1:3664317/150‑1 (MQ=255)
                                                                                                                                            ggggATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTGa       >  2:2468730/1‑150 (MQ=255)
                                                                                                                                               gATGGGGTAAACCTTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTGAgg     <  1:5014231/150‑1 (MQ=255)
                                                                                                                                               gATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTGAGGa    >  2:1782895/1‑150 (MQ=255)
                                                                                                                                                 tGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTGAGGaaa  <  2:3204505/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTGAGGAAA  >  CP009273/524349‑524642

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: