Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 365,572 G→A 5.8% pseudogene (116/534 nt) mhpC → pseudogene, IS30 insertion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273365,5720GA 5.8% 1210.3 / 43.1 327pseudogene (116/534 nt)mhpCpseudogene, IS30 insertion
Reads supporting (aligned to +/- strand):  ref base G (176/132);  new base A (12/7);  total (188/139)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 331 total aligned reads displayed.

TTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTT  >  CP009273/365423‑365720
                                                                                                                                                     |                                                                                                                                                    
ttGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAgg                                                                                                                                                      >  1:4956362/1‑150 (MQ=255)
ttGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAgg                                                                                                                                                      <  2:1920207/150‑1 (MQ=255)
 tGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGc                                                                                                                                                     >  2:488344/1‑150 (MQ=255)
   gCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGcga                                                                                                                                                   <  2:181904/150‑1 (MQ=255)
   gCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGcga                                                                                                                                                   <  2:2855010/150‑1 (MQ=255)
     tCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAATATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGcgacg                                                                                                                                                 <  2:459668/150‑1 (MQ=255)
     tCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGcgacg                                                                                                                                                 <  1:276020/150‑1 (MQ=255)
      cGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGcgacga                                                                                                                                                >  2:8479/1‑150 (MQ=255)
      cGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGcgacga                                                                                                                                                >  1:4475533/1‑150 (MQ=255)
       gTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGaa                                                                                                                                               <  2:3883889/150‑1 (MQ=255)
        ttAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGaaa                                                                                                                                              <  2:1056269/150‑1 (MQ=255)
           aGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCg                                                                                                                                           >  1:3418267/1‑150 (MQ=255)
           aGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCg                                                                                                                                           <  1:942758/150‑1 (MQ=255)
             cgcCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGtt                                                                                                                                         <  1:3755069/150‑1 (MQ=255)
              gcCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTg                                                                                                                                        >  2:4492197/1‑150 (MQ=255)
                 aGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTcc                                                                                                                                     >  1:3182096/1‑150 (MQ=255)
                 aGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTcc                                                                                                                                     >  1:2430802/1‑150 (MQ=255)
                 aGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCGGCGCATCCAATTTAATGACTGCGGACAAGGCGACGAAACCGGTGTcc                                                                                                                                     >  1:1805245/1‑150 (MQ=255)
                     cAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGACGACGAAACCGTTGTCctgct                                                                                                                                 >  1:4030648/1‑150 (MQ=255)
                      agagAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCctgctg                                                                                                                                >  2:3977387/1‑150 (MQ=255)
                       gagaACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCtgctgc                                                                                                                               >  1:1735882/1‑150 (MQ=255)
                       gagaACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCtgctgc                                                                                                                               <  2:3182830/150‑1 (MQ=255)
                         gaACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCAt                                                                                                                             <  1:2115271/150‑1 (MQ=255)
                         gaACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCAt                                                                                                                             >  2:2083411/1‑150 (MQ=255)
                           aCTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGACGACGAAACCGTTGTCCTGCTGCATgg                                                                                                                           >  1:1481407/1‑150 (MQ=255)
                               aaTATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTcc                                                                                                                       >  2:900877/1‑150 (MQ=255)
                               aaTATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTcc                                                                                                                       >  1:2051255/1‑150 (MQ=255)
                               aaTATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGGTGTCCTGCTGCATGGTTcc                                                                                                                       >  2:4360937/1‑150 (MQ=255)
                                 tatGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCgg                                                                                                                     <  2:4526932/150‑1 (MQ=255)
                                     cAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCg                                                                                                                 <  1:4281093/150‑1 (MQ=255)
                                     cAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCg                                                                                                                 >  1:4250807/1‑150 (MQ=255)
                                       ggAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTATAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCggg                                                                                                               >  1:1611280/1‑150 (MQ=255)
                                        gagaAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGt                                                                                                              <  2:3531169/150‑1 (MQ=255)
                                        gagaAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGt                                                                                                              >  2:1145962/1‑150 (MQ=255)
                                             gatgatGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTAc                                                                                                         >  1:4718072/1‑150 (MQ=255)
                                              atgatgAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTAct                                                                                                        <  1:643190/150‑1 (MQ=255)
                                              atgatgAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTAct                                                                                                        <  1:122471/150‑1 (MQ=255)
                                              atgatgAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTAct                                                                                                        <  1:2727242/150‑1 (MQ=255)
                                               tgatgaGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTActg                                                                                                       >  1:2768973/1‑150 (MQ=255)
                                               tgatgaGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCGGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTActg                                                                                                       >  1:2079662/1‑150 (MQ=255)
                                                gatgaGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTActgg                                                                                                      >  2:3187078/1‑150 (MQ=255)
                                                gatgaGTTATCAGCCACAAACCGAAGACGCCACCAGCCGTTTTCTTAATGTAGAAGAAGCGGGTAAAACGCTGCGCATACATTTTAATTACAGCGGACAAGGCCAAGAAACCGTTGTCCTGCTGCATGGTTAAGGCCCGGGTGCTaaatg                                                                                                      >  1:1400915/1‑146 (MQ=255)
                                                 atgaGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTActggc                                                                                                     <  2:572802/150‑1 (MQ=255)
                                                    aGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTActggctgg                                                                                                  <  1:3383549/150‑1 (MQ=255)
                                                      ttATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGc                                                                                                >  1:1847516/1‑150 (MQ=255)
                                                      ttATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCGGAATGTAGAAGAAGCGGGTAAAACGCTGCGCCTCCATTTTAATGACTGCGGCCAAGGCGACGAAACCCTTGTCCTGCGGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGc                                                                                                >  1:2874993/1‑150 (MQ=255)
                                                          cAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAAc                                                                                            >  1:3320585/1‑150 (MQ=255)
                                                           aGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGAGCTACTGGCTGGGCGAACt                                                                                           <  1:848808/150‑1 (MQ=255)
                                                             ccACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTc                                                                                         <  2:4078774/150‑1 (MQ=255)
                                                              cacaAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAATGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCa                                                                                        >  2:163840/1‑150 (MQ=255)
                                                              cacaAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGTCTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCGACTGGCTGGGCGGACTTCa                                                                                        >  1:39568/1‑150 (MQ=255)
                                                              cacaAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCa                                                                                        >  1:3649849/1‑150 (MQ=255)
                                                              cacaAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCa                                                                                        <  1:3054686/150‑1 (MQ=255)
                                                              cacaAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGACGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCa                                                                                        <  1:1002687/150‑1 (MQ=255)
                                                              cacaAACCGAAGCCGCCACCAGCCGTGTTCTGAATGTAGAAGAAGCGGGGAAAACGCTGCGCATCCCTTTTAAAGACTGCGGACCACGCGACGAAACCCGTGTCCTGCTGCCTGGGTCCCGCCCCGGTGCCACAGGCGGGGCGCACGTCa                                                                                        >  1:2341589/1‑150 (MQ=255)
                                                                 aaaCCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGcc                                                                                     >  2:1309872/1‑150 (MQ=255)
                                                                 aaaCCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGcc                                                                                     <  2:1495868/150‑1 (MQ=255)
                                                                 aaaCCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAAACGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGcc                                                                                     >  1:4391439/1‑150 (MQ=255)
                                                                  aaCCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCg                                                                                    >  1:1525773/1‑150 (MQ=255)
                                                                   aCCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGc                                                                                   >  1:1668232/1‑150 (MQ=255)
                                                                             ccaccaGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATc                                                                         <  2:3468924/150‑1 (MQ=255)
                                                                              caccaGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATcc                                                                        <  2:4354884/150‑1 (MQ=255)
                                                                                     cGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTa                                                                 >  2:2540372/1‑150 (MQ=255)
                                                                                     cGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTa                                                                 >  2:2696480/1‑150 (MQ=255)
                                                                                         ttCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAgg                                                             <  1:4190073/150‑1 (MQ=255)
                                                                                          tCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGc                                                            >  1:668725/1‑150 (MQ=255)
                                                                                            tGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCgg                                                          <  2:695120/150‑1 (MQ=255)
                                                                                            tGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCgg                                                          <  2:4159195/150‑1 (MQ=255)
                                                                                                  tAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGACGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATc                                                    <  1:351870/150‑1 (MQ=255)
                                                                                                   agaagaAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCg                                                   <  2:4160871/150‑1 (MQ=255)
                                                                                                    gaagaaGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCgg                                                  >  2:3321691/1‑150 (MQ=255)
                                                                                                    gaagaaGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCgg                                                  >  2:4641645/1‑150 (MQ=255)
                                                                                                           cGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGACCTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATcc                                           >  2:1304670/1‑150 (MQ=255)
                                                                                                           cGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATcc                                           >  2:2527827/1‑150 (MQ=255)
                                                                                                            gggTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCct                                          >  1:2558718/1‑150 (MQ=255)
                                                                                                               tAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCctgct                                       <  2:424196/150‑1 (MQ=255)
                                                                                                               tAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCctgct                                       >  2:1955441/1‑150 (MQ=255)
                                                                                                                         cgcATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCg                             >  1:1175340/1‑150 (MQ=255)
                                                                                                                          gcATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCgg                            <  1:2304006/150‑1 (MQ=255)
                                                                                                                           cATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCggg                           >  2:3762148/1‑150 (MQ=255)
                                                                                                                            aTCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGt                          >  2:3902524/1‑150 (MQ=255)
                                                                                                                            aTCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGt                          >  1:462416/1‑150 (MQ=255)
                                                                                                                              ccATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTg                        <  1:5069950/150‑1 (MQ=255)
                                                                                                                                aTTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTggg                      <  2:4502372/150‑1 (MQ=255)
                                                                                                                                  tttAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGc                    >  1:1529922/1‑150 (MQ=255)
                                                                                                                                   ttAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCa                   >  1:4074692/1‑150 (MQ=255)
                                                                                                                                   ttAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCa                   <  2:2316499/150‑1 (MQ=255)
                                                                                                                                    tAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCaa                  >  1:4862064/1‑150 (MQ=255)
                                                                                                                                    tAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCaa                  >  2:2382800/1‑150 (MQ=255)
                                                                                                                                    tAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCaa                  <  1:1197822/150‑1 (MQ=255)
                                                                                                                                      aTGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAaga                >  2:4410123/1‑150 (MQ=255)
                                                                                                                                         aCTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCg             >  2:4987296/1‑150 (MQ=255)
                                                                                                                                           tGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGAt           >  2:2016896/1‑150 (MQ=255)
                                                                                                                                             cGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTc         >  2:519949/1‑150 (MQ=255)
                                                                                                                                               gACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCgg       >  2:170607/1‑150 (MQ=255)
                                                                                                                                                  aaGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCg    >  1:1931838/1‑150 (MQ=255)
                                                                                                                                                  aaGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCg    <  2:2079625/150‑1 (MQ=255)
                                                                                                                                                  aaGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCg    <  2:2815062/150‑1 (MQ=255)
                                                                                                                                                  aaGACGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCTATTCGGTCg    >  2:2993141/1‑150 (MQ=255)
                                                                                                                                                   aGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGt   >  2:3644113/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTT  >  CP009273/365423‑365720

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