Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,523,802 T→G 100% pseudogene (1644/2037 nt) rhsE → pseudogene, Rhs family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,523,8020TG97.7% 846.8 / 1.0 212pseudogene (1644/2037 nt)rhsEpseudogene, Rhs family protein
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base G (119/88);  total (122/90)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 213 total aligned reads displayed.

TCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGA  >  CP009273/1523657‑1523946
                                                                                                                                                 |                                                                                                                                                
   ccccGGACCCGATTGGGGTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGAGTGAGAATATGACATCTGGTGGACTCATTTATGGCGTATCCGGTGTACCGGgatt                                                                                                                                             >  2:2817225/1‑148 (MQ=255)
    cccGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGgattg                                                                                                                                            <  2:3411013/148‑1 (MQ=255)
     ccGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGgattga                                                                                                                                           >  2:1432799/1‑148 (MQ=255)
      cGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGgattgat                                                                                                                                          <  2:587093/148‑1 (MQ=255)
       ggACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGgattgatt                                                                                                                                         >  2:1020329/1‑148 (MQ=255)
         aCCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATtgc                                                                                                                                       >  1:2408965/1‑148 (MQ=255)
           ccGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATtgctg                                                                                                                                     >  1:780822/1‑148 (MQ=255)
                  ggTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATGGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAg                                                                                                                              >  1:935189/1‑148 (MQ=255)
                   gTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGc                                                                                                                             <  1:2320957/148‑1 (MQ=255)
                   gTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGc                                                                                                                             <  1:408305/148‑1 (MQ=255)
                   gTTGAGAGGTGGAGGGAATATGTATCAGTATCCGTTGAATCCCATCCAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGc                                                                                                                             <  2:1257636/148‑1 (MQ=255)
                      gagaGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCAtt                                                                                                                          <  1:1755272/148‑1 (MQ=255)
                         agGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATGGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACt                                                                                                                       >  1:1625464/1‑148 (MQ=255)
                              gatggaATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAg                                                                                                                  >  2:3225923/1‑148 (MQ=255)
                                  gaATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCt                                                                                                              >  2:648389/1‑148 (MQ=255)
                                        gTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAg                                                                                                        >  1:3122290/1‑148 (MQ=255)
                                            cAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCg                                                                                                    >  2:2068989/1‑148 (MQ=255)
                                             aGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCgg                                                                                                   <  1:659579/148‑1 (MQ=255)
                                                aTCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTa                                                                                                >  2:2412040/1‑148 (MQ=255)
                                                  ccGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATg                                                                                              >  1:3377067/1‑148 (MQ=255)
                                                  ccGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATg                                                                                              >  2:1312125/1‑148 (MQ=255)
                                                  ccGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATg                                                                                              >  1:1708004/1‑148 (MQ=255)
                                                  ccGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATg                                                                                              >  2:3057100/1‑148 (MQ=255)
                                                    gTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGat                                                                                            >  1:124127/1‑148 (MQ=255)
                                                     ttGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGata                                                                                           <  1:2883287/148‑1 (MQ=255)
                                                        aaTCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATgg                                                                                        <  1:648389/148‑1 (MQ=255)
                                                        aaTCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATgg                                                                                        >  2:1471123/1‑148 (MQ=255)
                                                         aTCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGa                                                                                       >  1:2956272/1‑148 (MQ=255)
                                                         aTCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGa                                                                                       >  1:2688970/1‑148 (MQ=255)
                                                          tCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGAt                                                                                      <  1:2253304/148‑1 (MQ=255)
                                                          tCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGAt                                                                                      >  2:1049331/1‑148 (MQ=255)
                                                          tCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGAt                                                                                      >  2:1638391/1‑148 (MQ=255)
                                                          tCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGAt                                                                                      >  1:21768/1‑148 (MQ=255)
                                                           cccATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATg                                                                                     >  1:3358635/1‑148 (MQ=255)
                                                              aTACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTa                                                                                  <  1:3147483/148‑1 (MQ=255)
                                                               tACAAGTGCTAGACCCAATGGGGTTAGATGCGATTGAGAATATGCCATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTAt                                                                                 <  2:571688/148‑1 (MQ=255)
                                                                  aaGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTAttgt                                                                              <  2:1131377/148‑1 (MQ=255)
                                                                   aGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTAttgtt                                                                             <  2:653188/148‑1 (MQ=255)
                                                                   aGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTAttgtt                                                                             <  1:3218518/148‑1 (MQ=255)
                                                                   aGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTAttgtt                                                                             >  2:3389002/1‑148 (MQ=255)
                                                                    gTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTAttgttg                                                                            >  1:2687307/1‑148 (MQ=255)
                                                                      gATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGc                                                                          >  2:213598/1‑148 (MQ=255)
                                                                       aTAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCg                                                                         >  2:637878/1‑148 (MQ=255)
                                                                         aGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGa                                                                       <  2:3296552/148‑1 (MQ=255)
                                                                          gACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAg                                                                      <  2:2709300/148‑1 (MQ=255)
                                                                             ccAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTc                                                                   <  1:2340071/148‑1 (MQ=255)
                                                                             ccAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTc                                                                   >  1:2572680/1‑148 (MQ=255)
                                                                                aTGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATa                                                                >  1:1273567/1‑148 (MQ=255)
                                                                                aTGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATa                                                                >  2:1249416/1‑148 (MQ=255)
                                                                                aTGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGAAAACAGCATTACTAACATAGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGTCGGAGCTCaga                                                                >  1:3020917/1‑146 (MQ=255)
                                                                                    ggTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATgg                                                            <  2:2331560/148‑1 (MQ=255)
                                                                                      ttAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATgggg                                                          <  2:412908/148‑1 (MQ=255)
                                                                                          aTGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTAGCTTGTGTACCGGTAGTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGgccgc                                                      >  1:3447754/1‑148 (MQ=255)
                                                                                           tGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGgccgcc                                                     >  2:967533/1‑148 (MQ=255)
                                                                                             cGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGa                                                   <  2:419805/148‑1 (MQ=255)
                                                                                             cGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGa                                                   <  2:2331742/148‑1 (MQ=255)
                                                                                                ttGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGTCCGCCGATGc                                                >  1:247071/1‑148 (MQ=255)
                                                                                                ttGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGTCCGCCGATGc                                                >  1:617479/1‑148 (MQ=255)
                                                                                                ttGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGc                                                >  1:1258890/1‑148 (MQ=255)
                                                                                                ttGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGc                                                >  1:1381959/1‑148 (MQ=255)
                                                                                                ttGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGc                                                >  2:2915320/1‑148 (MQ=255)
                                                                                                 tGAGAATATGACCTCAGGTGGACTAATTTATGCCGTATCTGGTCTACCGGGAGTGATTGCTGCAAACATCATTACAAACAGTGCTTACCAGTTAGCTAATGATATGCATGCTATATTTGGCGGAGCACATAATGGGGCCGCCGATGCa                                               >  2:277775/1‑148 (MQ=255)
                                                                                                   agaATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAAt                                             >  2:1553560/1‑148 (MQ=255)
                                                                                                    gaATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATg                                            >  1:671288/1‑148 (MQ=255)
                                                                                                    gaATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATg                                            >  2:2505186/1‑148 (MQ=255)
                                                                                                      atatGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATgag                                          >  1:1410193/1‑148 (MQ=255)
                                                                                                         tGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACa                                       >  1:111982/1‑148 (MQ=255)
                                                                                                         tGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACa                                       <  1:3257459/148‑1 (MQ=255)
                                                                                                          gACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACAt                                      >  1:501726/1‑148 (MQ=255)
                                                                                                            cATCAGGTGGACTAATTTATGACGTATCTGGTGTACCGTGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTg                                    <  2:2140516/148‑1 (MQ=255)
                                                                                                             aTCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGt                                   <  1:948902/148‑1 (MQ=255)
                                                                                                             aTCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGt                                   >  1:2789100/1‑148 (MQ=255)
                                                                                                              tCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGtt                                  <  1:1471123/148‑1 (MQ=255)
                                                                                                                aGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTAc                                >  1:3446276/1‑148 (MQ=255)
                                                                                                                aGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTAc                                <  1:526126/148‑1 (MQ=255)
                                                                                                                  gTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACtt                              >  2:349292/1‑148 (MQ=255)
                                                                                                                   tGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTg                             >  2:1195467/1‑148 (MQ=255)
                                                                                                                    ggACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGa                            >  2:1935109/1‑148 (MQ=255)
                                                                                                                       cTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGAtgt                         >  1:360322/1‑148 (MQ=255)
                                                                                                                        tAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGAtgtg                        <  2:1997404/148‑1 (MQ=255)
                                                                                                                        tAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGAtgtg                        >  2:3210409/1‑148 (MQ=255)
                                                                                                                           tttATGCCGTATCTGGTGTCCCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCg                     >  1:817547/1‑148 (MQ=255)
                                                                                                                            ttATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGa                    <  1:1900646/148‑1 (MQ=255)
                                                                                                                             tATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGaa                   >  2:1373754/1‑148 (MQ=255)
                                                                                                                             tATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGaa                   <  2:1777214/148‑1 (MQ=255)
                                                                                                                             tATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGaa                   >  2:2254660/1‑148 (MQ=255)
                                                                                                                               tGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATg                 >  1:2355299/1‑148 (MQ=255)
                                                                                                                                gCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGa                <  2:2563673/148‑1 (MQ=255)
                                                                                                                                gCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGa                <  2:1991461/148‑1 (MQ=255)
                                                                                                                                 ccGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGAc               <  1:213598/148‑1 (MQ=255)
                                                                                                                                  cGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACt              >  1:1877990/1‑148 (MQ=255)
                                                                                                                                   gTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTa             >  2:3082118/1‑148 (MQ=255)
                                                                                                                                    tATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTaa            >  1:2526827/1‑148 (MQ=255)
                                                                                                                                       cTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGAc         >  1:1831180/1‑148 (MQ=255)
                                                                                                                                       cTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGAc         >  1:692551/1‑148 (MQ=255)
                                                                                                                                         ggTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACAt       <  2:780822/148‑1 (MQ=255)
                                                                                                                                         ggTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACAt       <  1:1697764/148‑1 (MQ=255)
                                                                                                                                           tgtACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACAgtt     >  2:104028/1‑148 (MQ=255)
                                                                                                                                             tACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTgg   >  2:1850639/1‑148 (MQ=255)
                                                                                                                                             tACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTgg   >  2:1071348/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGA  >  CP009273/1523657‑1523946

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: