Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,523,802 T→G 40.4% 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,8020TG40.4% 60.4 / 304.7 218pseudogene (1644/2037 nt)rhsEpseudogene, Rhs family protein
Reads supporting (aligned to +/- strand):  ref base T (70/60);  new base G (44/44);  total (114/104)
Fisher's exact test for biased strand distribution p-value = 5.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 223 total aligned reads displayed.

TATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGATCA  >  CP009273/1523655‑1523949
                                                                                                                                                   |                                                                                                                                                   
tatCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCCATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCt                                                                                                                                                     >  1:3241340/1‑148 (MQ=255)
  tCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTgg                                                                                                                                                   <  1:1595135/148‑1 (MQ=255)
   cACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGga                                                                                                                                                  <  2:2029962/148‑1 (MQ=255)
    aCCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGgat                                                                                                                                                 <  2:1354785/148‑1 (MQ=255)
    aCCCCGGACCCGAGTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATGGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGgat                                                                                                                                                 >  1:695585/1‑148 (MQ=255)
      cccGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGgattg                                                                                                                                               <  1:1702681/148‑1 (MQ=255)
      cccGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGgattg                                                                                                                                               <  1:845813/148‑1 (MQ=255)
      cccGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGgattg                                                                                                                                               <  1:1613340/148‑1 (MQ=255)
      cccGGACCCGATTGGGTTGAGAGGTGGATGGAACATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGgattg                                                                                                                                               <  1:27888/148‑1 (MQ=255)
       ccGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGgattga                                                                                                                                              <  2:603273/148‑1 (MQ=255)
         ggACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGgattgatt                                                                                                                                            <  2:2760373/148‑1 (MQ=255)
          gACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGattgattg                                                                                                                                           <  2:1917656/148‑1 (MQ=255)
          gACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGattgattg                                                                                                                                           <  2:2442271/148‑1 (MQ=255)
           aCCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATtgc                                                                                                                                          >  1:1042809/1‑148 (MQ=255)
                    ggTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAg                                                                                                                                 <  2:2465660/148‑1 (MQ=255)
                    ggTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAg                                                                                                                                 <  1:1885543/148‑1 (MQ=255)
                    ggTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAg                                                                                                                                 <  1:186350/148‑1 (MQ=255)
                     gTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTTCCTGGATTGATTGCTGCAAACAGc                                                                                                                                <  1:740049/148‑1 (MQ=255)
                      ttGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCa                                                                                                                               <  1:1094725/148‑1 (MQ=255)
                       tGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCAt                                                                                                                              >  2:92519/1‑148 (MQ=255)
                          gagGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCACTGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTAc                                                                                                                           >  2:3291250/1‑148 (MQ=255)
                             gTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTaa                                                                                                                        <  1:1270051/148‑1 (MQ=255)
                                gatggaATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAg                                                                                                                     <  1:2282178/148‑1 (MQ=255)
                                   ggaATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGc                                                                                                                  >  2:3250424/1‑148 (MQ=255)
                                     aaTATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCtt                                                                                                                <  1:246641/148‑1 (MQ=255)
                                      atatGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTa                                                                                                               >  2:2401653/1‑148 (MQ=255)
                                      atatGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTa                                                                                                               <  1:2210494/148‑1 (MQ=255)
                                       tatGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTAc                                                                                                              >  1:1121771/1‑148 (MQ=255)
                                       tatGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTAc                                                                                                              >  1:1337335/1‑148 (MQ=255)
                                         tGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCa                                                                                                            >  2:2105889/1‑148 (MQ=255)
                                              cAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCg                                                                                                       >  2:2645921/1‑148 (MQ=255)
                                              cAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCg                                                                                                       >  1:1449859/1‑148 (MQ=255)
                                                gTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGt                                                                                                     >  2:1631426/1‑148 (MQ=255)
                                                 tATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGtt                                                                                                    <  1:665733/148‑1 (MQ=255)
                                                  aTCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTa                                                                                                   <  1:265786/148‑1 (MQ=255)
                                                   tCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGAGTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTAt                                                                                                  >  2:2225859/1‑148 (MQ=255)
                                                    ccGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATg                                                                                                 >  2:453624/1‑148 (MQ=255)
                                                     cGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGa                                                                                                <  1:171420/148‑1 (MQ=255)
                                                        tGAATCCCATACAAGTGAGAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGatat                                                                                             <  1:1532236/148‑1 (MQ=255)
                                                          aaTCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATgg                                                                                           <  1:612881/148‑1 (MQ=255)
                                                             cccATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATg                                                                                        >  1:390013/1‑148 (MQ=255)
                                                              ccagtcAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATTACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCAGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGc                                                                                       <  2:798983/143‑1 (MQ=255)
                                                                aTACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTa                                                                                     <  2:436267/148‑1 (MQ=255)
                                                                   cAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTAttg                                                                                  <  1:2579384/148‑1 (MQ=255)
                                                                   cAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTAttg                                                                                  <  1:2503640/148‑1 (MQ=255)
                                                                    aaGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTAttgt                                                                                 <  2:671718/148‑1 (MQ=255)
                                                                     aGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTAttgtt                                                                                >  1:849094/1‑148 (MQ=255)
                                                                      gTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTAttgttg                                                                               >  2:2441295/1‑148 (MQ=255)
                                                                       tGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTgg                                                                              >  1:3093315/1‑148 (MQ=255)
                                                                       tGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTgg                                                                              >  2:3079039/1‑148 (MQ=255)
                                                                          tAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCgg                                                                           >  1:1179799/1‑148 (MQ=255)
                                                                          tAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCgg                                                                           <  2:1276648/148‑1 (MQ=255)
                                                                          tAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCgg                                                                           >  2:2631801/1‑148 (MQ=255)
                                                                             aCCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGc                                                                        >  1:2304002/1‑148 (MQ=255)
                                                                                 aaTGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCAt                                                                    >  1:1850742/1‑148 (MQ=255)
                                                                                  aTGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATa                                                                   <  1:2960757/148‑1 (MQ=255)
                                                                                  aTGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATa                                                                   >  2:1897633/1‑148 (MQ=255)
                                                                                  aTGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATa                                                                   <  2:1041721/148‑1 (MQ=255)
                                                                                     gggTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATg                                                                >  2:1229935/1‑148 (MQ=255)
                                                                                     gggTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATg                                                                <  1:1936031/148‑1 (MQ=255)
                                                                                      ggTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATgg                                                               <  2:1170184/148‑1 (MQ=255)
                                                                                      ggTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATgg                                                               >  1:1625693/1‑148 (MQ=255)
                                                                                        ttAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATgggg                                                             <  2:2829808/148‑1 (MQ=255)
                                                                                        ttAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATgggg                                                             <  2:2831958/148‑1 (MQ=255)
                                                                                         tataTGCGATTGAGAATTTGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGCCGTGATTGCTGCAAACAGCATTACTAACCGTGCTTACCCGTGCGGTTATGAGTTGGATGCTCGTGTTGGCGGAGCTCATAAGGGGgc                                                            >  2:1842213/4‑148 (MQ=255)
                                                                                           gATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGgccg                                                          >  2:183635/1‑148 (MQ=255)
                                                                                             tGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGgccgcc                                                        <  1:2800499/148‑1 (MQ=255)
                                                                                                     agaATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAAt                                                >  1:476924/1‑148 (MQ=255)
                                                                                                       aaTATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATga                                              <  2:1329182/148‑1 (MQ=255)
                                                                                                        atatGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATgag                                             <  1:980150/148‑1 (MQ=255)
                                                                                                         tatGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATgaga                                            >  1:3263340/1‑148 (MQ=255)
                                                                                                          atGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGAc                                           >  2:2131798/1‑148 (MQ=255)
                                                                                                          atGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGAc                                           <  2:125003/148‑1 (MQ=255)
                                                                                                           tGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACa                                          <  2:2918828/148‑1 (MQ=255)
                                                                                                           tGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCAGAATGGGGCCGCCGATGCAATGAGACa                                          >  1:1723318/1‑148 (MQ=255)
                                                                                                             aCATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACAtt                                        <  1:433123/148‑1 (MQ=255)
                                                                                                               aTCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGt                                      <  2:768903/148‑1 (MQ=255)
                                                                                                               aTCAGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGt                                      >  2:2833688/1‑148 (MQ=255)
                                                                                                                tCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGtt                                     >  2:534857/1‑148 (MQ=255)
                                                                                                                  aGGTGGACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTAc                                   <  2:1576305/148‑1 (MQ=255)
                                                                                                                   gggggACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACt                                  >  2:2670061/4‑148 (MQ=255)
                                                                                                                    gTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCTCCGATGCAATGAGACATTGTTACtt                                 >  1:2015535/1‑148 (MQ=255)
                                                                                                                       gACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGAt                              <  2:2070803/148‑1 (MQ=255)
                                                                                                                       gACTAATTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGAt                              >  2:774262/1‑148 (MQ=255)
                                                                                                                        aCTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGAtg                             >  1:2113603/1‑148 (MQ=255)
                                                                                                                            aTTTATGTCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTc                         >  2:1403174/1‑148 (MQ=255)
                                                                                                                            aTTTATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTc                         >  1:1989673/1‑148 (MQ=255)
                                                                                                                               tATGCCGTATCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGaa                      <  2:1790223/148‑1 (MQ=255)
                                                                                                                                ctGCCGTATCTGGTGTACCGGGAGTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAAt                     <  2:1710572/147‑1 (MQ=255)
                                                                                                                                       aTCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAg              >  1:3304632/1‑148 (MQ=255)
                                                                                                                                       aTCTGGTGTACCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAg              >  1:1859583/1‑148 (MQ=255)
                                                                                                                                         cTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGAc            >  1:1052369/1‑148 (MQ=255)
                                                                                                                                          tGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACa           >  1:1692536/1‑148 (MQ=255)
                                                                                                                                           ggTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACAt          >  2:657654/1‑148 (MQ=255)
                                                                                                                                             tgtACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACAttt        <  2:1679466/148‑1 (MQ=255)
                                                                                                                                                aCCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGa     <  2:1448735/148‑1 (MQ=255)
                                                                                                                                                aCCGGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGa     <  2:1449859/148‑1 (MQ=255)
                                                                                                                                                  cTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGATc   >  1:742167/1‑148 (MQ=255)
                                                                                                                                                  cTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGATc   >  1:1087410/1‑148 (MQ=255)
                                                                                                                                                   tGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGATca  <  1:1229935/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGATCA  >  CP009273/1523655‑1523949

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: