Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 365,572 G→A 7.0% 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 7.0% 868.0 / 38.2 255pseudogene (116/534 nt)mhpCpseudogene, IS30 insertion
Reads supporting (aligned to +/- strand):  ref base G (118/119);  new base A (11/7);  total (129/126)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 257 total aligned reads displayed.

TGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTA  >  CP009273/365424‑365721
                                                                                                                                                    |                                                                                                                                                     
tGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGc                                                                                                                                                      >  1:903937/1‑150 (MQ=255)
 ggCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGcg                                                                                                                                                     <  2:4258797/150‑1 (MQ=255)
     cGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGcgacga                                                                                                                                                 <  2:4154326/150‑1 (MQ=255)
       ttAAGCGCCAGAACAGAGAACTGAATATGCCGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGTATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGaaa                                                                                                                                               >  2:4388483/1‑150 (MQ=255)
       ttAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGaaa                                                                                                                                               >  1:2598396/1‑150 (MQ=255)
       ttAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGaaa                                                                                                                                               >  1:724643/1‑150 (MQ=255)
        tAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAAc                                                                                                                                              <  1:1902927/150‑1 (MQ=255)
        tAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAAc                                                                                                                                              <  1:3667095/150‑1 (MQ=255)
        tAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAAc                                                                                                                                              <  1:4105893/150‑1 (MQ=255)
         aaGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAAcc                                                                                                                                             >  1:1571812/1‑150 (MQ=255)
          aGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCg                                                                                                                                            >  1:4081497/1‑150 (MQ=255)
             gcCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTg                                                                                                                                         <  2:3542246/150‑1 (MQ=255)
             gcCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGATATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGACACCGTTg                                                                                                                                         <  2:252431/150‑1 (MQ=255)
              ccAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGt                                                                                                                                        <  1:1813673/150‑1 (MQ=255)
                  aaCAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCctg                                                                                                                                    >  1:1206793/1‑150 (MQ=255)
                  aaCAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCctg                                                                                                                                    >  1:2985569/1‑150 (MQ=255)
                     agagAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCctgctg                                                                                                                                 <  2:600312/150‑1 (MQ=255)
                        gaACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCAt                                                                                                                              >  1:34685/1‑150 (MQ=255)
                         aaCTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATg                                                                                                                             >  1:1069831/1‑150 (MQ=255)
                            tGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGtt                                                                                                                          >  2:1623118/1‑150 (MQ=255)
                            tGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGACGACGAAACCGTTGTCCTGCTGCATGGtt                                                                                                                          <  2:3739906/150‑1 (MQ=255)
                              aaTATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGACGACGAAACCGTTGTCCTGCTGCATGGTTcc                                                                                                                        <  2:3045827/150‑1 (MQ=255)
                              aaTATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGACGACGAAACCGTTGTCCTGCTGCATGGTTcc                                                                                                                        <  2:2948365/150‑1 (MQ=255)
                               atatGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCg                                                                                                                       <  2:2601239/150‑1 (MQ=255)
                                tatGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCgg                                                                                                                      <  2:4505374/150‑1 (MQ=255)
                                    cAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTTCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCg                                                                                                                  >  2:4406327/1‑150 (MQ=255)
                                     aGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCgg                                                                                                                 >  2:1778405/1‑150 (MQ=255)
                                      ggAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCggg                                                                                                                >  1:3866701/1‑150 (MQ=255)
                                       gagaAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGt                                                                                                               <  1:768628/150‑1 (MQ=255)
                                        agaagaTGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTg                                                                                                              <  1:2507455/150‑1 (MQ=255)
                                        agaagaTGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTg                                                                                                              <  2:1028051/150‑1 (MQ=255)
                                        agaagaTGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTg                                                                                                              >  1:2353728/1‑150 (MQ=255)
                                             atgatgAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTAct                                                                                                         <  2:1206793/150‑1 (MQ=255)
                                             atgatgAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTAct                                                                                                         <  2:2985569/150‑1 (MQ=255)
                                                  gaGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTActggctg                                                                                                    >  1:884618/1‑150 (MQ=255)
                                                   aGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTActggctgg                                                                                                   >  2:1428626/1‑150 (MQ=255)
                                                     ttATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGc                                                                                                 >  1:4459800/1‑150 (MQ=255)
                                                     ttATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGc                                                                                                 >  1:459253/1‑150 (MQ=255)
                                                     ttATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGc                                                                                                 >  2:4348266/1‑150 (MQ=255)
                                                        tCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGaa                                                                                              >  1:1473542/1‑150 (MQ=255)
                                                        tCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGaa                                                                                              <  2:1678663/150‑1 (MQ=255)
                                                         cAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAAc                                                                                             >  1:3895001/1‑150 (MQ=255)
                                                             cacaAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGACGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCa                                                                                         <  2:901258/150‑1 (MQ=255)
                                                                   ccGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCa                                                                                   <  1:2949604/150‑1 (MQ=255)
                                                                   ccGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCa                                                                                   <  1:2950352/150‑1 (MQ=255)
                                                                      aaGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAata                                                                                >  1:2814418/1‑150 (MQ=255)
                                                                       aGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGATGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAatat                                                                               <  2:2039973/150‑1 (MQ=255)
                                                                         ccgccACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCg                                                                             >  2:3721215/1‑150 (MQ=255)
                                                                         ccgccACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCg                                                                             <  1:2537440/150‑1 (MQ=255)
                                                                           gccACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGAt                                                                           <  2:4116481/150‑1 (MQ=255)
                                                                           gccACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGAt                                                                           <  1:205932/150‑1 (MQ=255)
                                                                            ccaccaGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAAAATCGATc                                                                          >  1:1373861/1‑150 (MQ=255)
                                                                             caccaGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATcc                                                                         <  2:2280328/150‑1 (MQ=255)
                                                                             caccaGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATcc                                                                         >  2:3345900/1‑150 (MQ=255)
                                                                              accaGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCg                                                                        >  1:3333102/1‑150 (MQ=255)
                                                                                 aGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTg                                                                     >  2:1408732/1‑150 (MQ=255)
                                                                                 aGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTg                                                                     <  2:2287142/150‑1 (MQ=255)
                                                                                  gCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTgg                                                                    <  1:1900868/150‑1 (MQ=255)
                                                                                     gTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTag                                                                 >  1:2809780/1‑150 (MQ=255)
                                                                                      ttttCTGAATGTAGAAGAAGCGGGTAAAACGCTGCTCATCCATTATAATGACTGCGGACCAGGCGACGAAACCGTTGTCCTGCTGAATGGTTCCGGCCCGGGTGCTCCTGGCTGGGCGAACTTCAGCCGCAAAATCGATCCGCTGTTaga                                                                >  2:2614601/1‑150 (MQ=255)
                                                                                      ttttCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAACCCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTaga                                                                <  2:4005997/150‑1 (MQ=255)
                                                                                      ttttCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTaga                                                                >  1:4361326/1‑150 (MQ=255)
                                                                                        ttCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAgg                                                              >  1:1174464/1‑150 (MQ=255)
                                                                                         tCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGc                                                             <  1:2559267/150‑1 (MQ=255)
                                                                                           tGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCgg                                                           <  1:1698370/150‑1 (MQ=255)
                                                                                            gAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCggg                                                          <  1:4208048/150‑1 (MQ=255)
                                                                                            gAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCggg                                                          <  1:4208020/150‑1 (MQ=255)
                                                                                             aaTGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGc                                                         >  1:1990804/1‑150 (MQ=255)
                                                                                                  agaagaAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCg                                                    >  1:611815/1‑150 (MQ=255)
                                                                                                  agaagaAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCg                                                    >  2:241778/1‑150 (MQ=255)
                                                                                                       aagCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGa                                               >  2:1913199/1‑150 (MQ=255)
                                                                                                       aagCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGa                                               >  1:3701780/1‑150 (MQ=255)
                                                                                                       aagCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGa                                               >  1:2291923/1‑150 (MQ=255)
                                                                                                       aagCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGa                                               >  1:2635187/1‑150 (MQ=255)
                                                                                                            ggTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCctg                                          <  2:2145867/150‑1 (MQ=255)
                                                                                                              tAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCctgct                                        <  1:1623118/150‑1 (MQ=255)
                                                                                                              tAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCctgct                                        <  1:383502/150‑1 (MQ=255)
                                                                                                                       gcgcATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTcc                               <  1:3061931/150‑1 (MQ=255)
                                                                                                                        cgcATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCg                              >  2:3170183/1‑150 (MQ=255)
                                                                                                                          cATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCggg                            >  1:4140794/1‑150 (MQ=255)
                                                                                                                           aTCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGt                           >  1:2269805/1‑150 (MQ=255)
                                                                                                                               aTTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTggg                       <  1:4210903/150‑1 (MQ=255)
                                                                                                                               aTTTTAATGACTGCGGACAAGGCGACGAAACCGGTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTggg                       <  2:554325/150‑1 (MQ=255)
                                                                                                                                 tttAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGc                     >  1:3534168/1‑150 (MQ=255)
                                                                                                                                  ttAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCa                    <  1:1895878/150‑1 (MQ=255)
                                                                                                                                     aTGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAaga                 <  1:734819/150‑1 (MQ=255)
                                                                                                                                      tGACTGCGGACAAGGCGACGACACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAagag                <  2:4281202/150‑1 (MQ=255)
                                                                                                                                      tGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAagag                >  2:2882086/1‑150 (MQ=255)
                                                                                                                                       gACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGc               <  1:2279366/150‑1 (MQ=255)
                                                                                                                                       gACTGCGGACAAGACGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGc               >  2:465405/1‑150 (MQ=255)
                                                                                                                                          tGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGAt            <  1:429477/150‑1 (MQ=255)
                                                                                                                                          tGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGAt            <  2:3595834/150‑1 (MQ=255)
                                                                                                                                           gCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGAtt           >  2:651224/1‑150 (MQ=255)
                                                                                                                                           gCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGAtt           >  2:651277/1‑150 (MQ=255)
                                                                                                                                            cGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTTCTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTc          <  2:2366897/150‑1 (MQ=255)
                                                                                                                                             ggACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCg         >  2:2585736/1‑150 (MQ=255)
                                                                                                                                              gACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCgg        <  1:1980070/150‑1 (MQ=255)
                                                                                                                                                cAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTc      >  2:1504570/1‑150 (MQ=255)
                                                                                                                                                  aGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGt    >  1:1004920/1‑150 (MQ=255)
                                                                                                                                                   ggCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGtt   <  1:3214508/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTA  >  CP009273/365424‑365721

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: