Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,030,376 C→A 8.8% E364* (GAA→TAA)  recJ ← ssDNA exonuclease, 5' ‑‑> 3'‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,030,3760CA 8.8% 1193.4 / 54.1 316E364* (GAA→TAA) recJssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base C (137/151);  new base A (11/17);  total (148/168)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 325 total aligned reads displayed.

CTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATCGAGTTCATTTGCCAGCACGCGCGCTT  >  CP009273/3030227‑3030525
                                                                                                                                                     |                                                                                                                                                     
   cAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCggg                                                                                                                                                    <  1:146226/150‑1 (MQ=255)
   cAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCggg                                                                                                                                                    >  2:3673674/1‑150 (MQ=255)
   cAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCggg                                                                                                                                                    <  1:2014292/150‑1 (MQ=255)
    aGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCgggg                                                                                                                                                   >  1:3471036/1‑150 (MQ=255)
      cgcATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTg                                                                                                                                                 >  2:2455378/1‑150 (MQ=255)
        cATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGAt                                                                                                                                               <  1:1232372/150‑1 (MQ=255)
         aTCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATa                                                                                                                                              >  2:1156199/1‑150 (MQ=255)
         aTCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATa                                                                                                                                              >  1:884838/1‑150 (MQ=255)
         aTCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATa                                                                                                                                              >  1:3154942/1‑150 (MQ=255)
           cACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACa                                                                                                                                            >  2:2476445/1‑150 (MQ=255)
           cACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACa                                                                                                                                            <  1:2216917/150‑1 (MQ=255)
            aCGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACAt                                                                                                                                           >  2:402474/1‑150 (MQ=255)
             cGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACAtt                                                                                                                                          <  2:3003909/150‑1 (MQ=255)
             cGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACAtt                                                                                                                                          <  1:4634397/150‑1 (MQ=255)
              gCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTg                                                                                                                                         <  2:3854519/150‑1 (MQ=255)
                 tatGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCa                                                                                                                                      <  2:4397812/150‑1 (MQ=255)
                 tatGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCa                                                                                                                                      <  2:458550/150‑1 (MQ=255)
                        ccccTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCccgcc                                                                                                                               >  1:4124913/1‑150 (MQ=255)
                        ccccTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCccgcc                                                                                                                               <  1:2638188/150‑1 (MQ=255)
                            tGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGt                                                                                                                           >  2:1330942/1‑150 (MQ=255)
                              aaTGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAg                                                                                                                         <  2:288874/150‑1 (MQ=255)
                              aaTGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAg                                                                                                                         <  2:288414/150‑1 (MQ=255)
                                tGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCg                                                                                                                       >  2:897307/1‑150 (MQ=255)
                                 ggAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCgt                                                                                                                      >  2:3582355/1‑150 (MQ=255)
                                 ggAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCgt                                                                                                                      <  1:525123/150‑1 (MQ=255)
                                  gAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCgtg                                                                                                                     <  2:4618720/150‑1 (MQ=255)
                                  gAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCgtg                                                                                                                     <  1:810619/150‑1 (MQ=255)
                                    gCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTc                                                                                                                   <  2:4237540/150‑1 (MQ=255)
                                     cGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTAGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCa                                                                                                                  >  2:925251/1‑150 (MQ=255)
                                      gACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCAc                                                                                                                 <  1:2476445/150‑1 (MQ=255)
                                      gACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCAc                                                                                                                 <  2:2268271/150‑1 (MQ=255)
                                           gAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTg                                                                                                            >  2:3174480/1‑150 (MQ=255)
                                             aCCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTgcg                                                                                                          >  1:2718890/1‑150 (MQ=255)
                                                tttCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTc                                                                                                       >  2:3593292/1‑150 (MQ=255)
                                                     gCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGtt                                                                                                  >  1:2180252/1‑150 (MQ=255)
                                                        tACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTtct                                                                                               <  1:651195/150‑1 (MQ=255)
                                                         aCCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTtctc                                                                                              <  1:314942/150‑1 (MQ=255)
                                                            gTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGca                                                                                           >  2:1810872/1‑150 (MQ=255)
                                                            gTCAACTGCTGGCGCAAACTCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAAAGCCCTGATGCCATTCGGTGTGATACATTGCCAGCCCGCCGGGTAGGGTGTCACGGCTGCGCTCCAGTTTCTCGca                                                                                           <  1:1802790/150‑1 (MQ=255)
                                                             tCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCGGTTTCTCGcac                                                                                          >  1:451837/1‑150 (MQ=255)
                                                              cACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGcaca                                                                                         >  1:236336/1‑150 (MQ=255)
                                                                ccTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAgg                                                                                       <  1:304188/150‑1 (MQ=255)
                                                                      ggCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAgg                                                                                 <  2:4082396/150‑1 (MQ=255)
                                                                      ggCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAgg                                                                                 <  2:3291817/150‑1 (MQ=255)
                                                                         gcAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCt                                                                              <  2:2391546/150‑1 (MQ=255)
                                                                          cAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCtt                                                                             <  2:3131932/150‑1 (MQ=255)
                                                                          cAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCtt                                                                             <  2:37660/150‑1 (MQ=255)
                                                                           aaaCGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTc                                                                            <  1:2440856/150‑1 (MQ=255)
                                                                                cgATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTAGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAAttt                                                                       <  2:4189683/150‑1 (MQ=255)
                                                                                cgATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTAGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAAttt                                                                       >  1:3668152/1‑150 (MQ=255)
                                                                                 gATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTg                                                                      >  1:592691/1‑150 (MQ=255)
                                                                                      ccGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTc                                                                 <  2:4119442/150‑1 (MQ=255)
                                                                                        ggACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCt                                                               >  2:2204272/1‑150 (MQ=255)
                                                                                        ggACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCt                                                               >  2:4439499/1‑150 (MQ=255)
                                                                                         gACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCtt                                                              <  2:3599651/150‑1 (MQ=255)
                                                                                         gACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCtt                                                              <  2:1055673/150‑1 (MQ=255)
                                                                                         gACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCtt                                                              <  2:504294/150‑1 (MQ=255)
                                                                                         gACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCtt                                                              >  2:652434/1‑150 (MQ=255)
                                                                                          aCGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTg                                                             >  1:4048902/1‑150 (MQ=255)
                                                                                          aCGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTg                                                             <  1:4619901/150‑1 (MQ=255)
                                                                                             gTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTc                                                          <  2:775222/150‑1 (MQ=255)
                                                                                               gAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGa                                                        >  2:4517215/1‑150 (MQ=255)
                                                                                                 aaaCGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGAtc                                                      <  1:4247784/150‑1 (MQ=255)
                                                                                                  aaCGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGAtct                                                     >  1:3444/1‑150 (MQ=255)
                                                                                                  aaCGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGAtct                                                     <  1:3784574/150‑1 (MQ=255)
                                                                                                   aCGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGAtctc                                                    <  2:3165143/150‑1 (MQ=255)
                                                                                                    cGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATctct                                                   <  1:1194611/150‑1 (MQ=255)
                                                                                                    cGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATctct                                                   <  2:1477052/150‑1 (MQ=255)
                                                                                                      ctctTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCttt                                                 <  1:504959/150‑1 (MQ=255)
                                                                                                      ctctTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTAGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCttt                                                 <  2:1693248/150‑1 (MQ=255)
                                                                                                        ctTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTAACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTc                                               <  2:2192015/150‑1 (MQ=255)
                                                                                                         tttGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTcg                                              >  1:2107235/1‑150 (MQ=255)
                                                                                                             aTGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTc                                          >  1:4475229/1‑150 (MQ=255)
                                                                                                             aTGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTc                                          >  2:2677026/1‑150 (MQ=255)
                                                                                                              tGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCt                                         >  2:3270839/1‑150 (MQ=255)
                                                                                                               gcgcGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTg                                        <  2:3722924/150‑1 (MQ=255)
                                                                                                               gcgcGAAGCCAGAATACCGACAACGCCCTGATGCCATTAGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTg                                        <  2:3389958/150‑1 (MQ=255)
                                                                                                                   gAAGCCAGAATACCGACAACGCCCTGATGCCATTAGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGttt                                    >  2:3464829/1‑150 (MQ=255)
                                                                                                                    aaGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTa                                   <  1:2455384/150‑1 (MQ=255)
                                                                                                                    aaGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTa                                   >  2:4106510/1‑150 (MQ=255)
                                                                                                                    aaGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTGTCGCGTCTGGTTTa                                   >  2:2094528/1‑150 (MQ=255)
                                                                                                                     aGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAg                                  <  1:3305908/150‑1 (MQ=255)
                                                                                                                       ccAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAgcg                                <  1:1966383/150‑1 (MQ=255)
                                                                                                                        cAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAgcgc                               >  2:1115834/1‑150 (MQ=255)
                                                                                                                         aGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCa                              >  2:1262999/1‑150 (MQ=255)
                                                                                                                          gAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCAt                             >  2:1458949/1‑150 (MQ=255)
                                                                                                                          gAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCAt                             <  2:1393404/150‑1 (MQ=255)
                                                                                                                           aaTACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATc                            <  2:2630837/150‑1 (MQ=255)
                                                                                                                            aTACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATCg                           <  1:1143813/150‑1 (MQ=255)
                                                                                                                                    aaCGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATCGAGTTCAtt                   >  2:190620/1‑150 (MQ=255)
                                                                                                                                     aCGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATCGAGTTCAttt                  <  1:1636671/150‑1 (MQ=255)
                                                                                                                                         ccTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATCGAGTTCATTTGCCa              <  1:2860362/150‑1 (MQ=255)
                                                                                                                                          cTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATCGAGTTCATTTGCCAg             <  1:3593292/150‑1 (MQ=255)
                                                                                                                                                 cATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATCGAGTTCATTTGCCAGCAcgcgc      >  1:3985626/1‑150 (MQ=255)
                                                                                                                                                 cATTAGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATCGAGTTCATTTGCCAGCAcgcgc      >  1:3704468/1‑150 (MQ=255)
                                                                                                                                                   ttCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATCGAGTTCATTTGCCAGCACgcgcgc    >  2:2521737/1‑150 (MQ=255)
                                                                                                                                                    tCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATCGAGTTCATTTGCCAGCACGCGCGCt   <  2:197433/150‑1 (MQ=255)
                                                                                                                                                     cGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATCGAGTTCATTTGCCAGCACGCGCGCtt  <  2:4307384/150‑1 (MQ=255)
                                                                                                                                                     cGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATCGAGTTCATTTGCCAGCACGCGCGCtt  <  1:4490789/150‑1 (MQ=255)
                                                                                                                                                     cGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATCGAGTTCATTTGCCAGCACGCGCGCtt  <  1:1573588/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATCGAGTTCATTTGCCAGCACGCGCGCTT  >  CP009273/3030227‑3030525

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: