Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,172 T→A 5.2% intergenic (‑106/+489) lacZ ← / ← lacI pseudogene, truncated/DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273360,1720TA 5.2% 864.0 / 10.9 291intergenic (‑106/+489)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (83/187);  new base A (9/10);  total (92/199)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 296 total aligned reads displayed.

CCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCT  >  CP009273/360024‑360316
                                                                                                                                                    |                                                                                                                                                
ccAGCGGCTTACCATCCAGCGCCACCATCCCGTGCCGGAGCTCCGATGATACGCTGACAAACATGAGAATTCCGATCCCCAATTCGTGGCTGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAt                                                                                                                                                  <  1:1440381/149‑1 (MQ=37)
 ttcaccgaccaacaacagataaaacgaaaggcccagtctttcgactgagccttgcgttttatttgatgccgggggacgtCATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCGCCGGGCCGGTGCTTCGCAACGGTCCAATc                                                                                                                                                 <  2:3303321‑M1/70‑1 (MQ=255)
  aGCGGCTTACCATCCAGCGCCACCATCCCGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCACCGTTCAAATcc                                                                                                                                                >  2:2403581/1‑149 (MQ=255)
     ggCTTACCATCCAGCGCCACCCTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                                             <  1:2262212/149‑1 (MQ=255)
     ggCTAACCACCCAGCCCCACCACCCAGTGCAGGAGCTCCGATTATACGCTGTAAAACATGAGAATTCCGAACCCCATTCCCTGGCTGTTTTGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                                             <  1:3244886/149‑1 (MQ=21)
      gCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTAGCAAAGTTAAAATCAGcac                                                                                                                                            >  1:2144851/1‑147 (MQ=37)
      cgacaaacaccagataaaacgaccggcccagtctttcgactgagcctttctttttatttgatgcctggggaataAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                                                            >  1:4273440‑M1/75‑149 (MQ=255)
       cTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTcc                                                                                                                                           >  2:1888868/1‑149 (MQ=255)
        taaccaaccaGCGCCACCACCCAGTGCTGTAGCTCCGATTATATGGTGTCAAACATGAGAATTCCGAGCCCCATTTCCTGGCTGTTTTTGGGGGGGGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                                          <  1:4389899/147‑1 (MQ=21)
         caaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgccgggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                                                                                         <  2:3005885‑M1/78‑1 (MQ=255)
             caacagataaaacgaaaggcccagtctttcgactgagcctggcggtgtatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCCACGGTCAAAACCGCACCcggcgg                                                                                                                                     <  2:3931066‑M1/82‑1 (MQ=255)
             aTCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                                                                                                                     <  2:2186554/149‑1 (MQ=255)
              tCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                                    <  2:86732/149‑1 (MQ=255)
              aacagataaaccgacacgcccaggcgctcgtctgagcctctcgatttatttgttgcctggggaattAATTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                                    >  2:2985535‑M1/67‑149 (MQ=255)
               ccAGCGCCACCACCCTGTGCTGGAGCTCCTATTATAAGCTGTCAAACTTGAGATTACCGACCCCCAAGCCGTGGCTGTTTTGGGGGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                                   <  1:488661/149‑1 (MQ=18)
               acagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                                   >  2:2440444‑M1/66‑149 (MQ=255)
               acagaaaaaacgaacggcccagtctttctactgagcctttctttttatttgaggcgtgggtaataAATTCCTGGATTTTTTTGGCGGGCGTCCTCCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                                   >  1:590623‑M1/66‑149 (MQ=255)
                cagataaaacgaaagccccagtctttctactgaccctttctttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                                  >  1:3403390‑M1/65‑149 (MQ=255)
                cAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                                  <  2:1803650/149‑1 (MQ=255)
                      caccaTCCAGTGCAGGAGCTCCGTTGATAAGGTGTCAAACATGTGAATTCCGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGATTTGTCCt                                                                                                                            <  2:2243349/149‑1 (MQ=21)
                       accaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTTCTTCGCAACCTTCAAATCCCCGCCCGGCGGATTTGTccga                                                                                                                           >  1:3774846/1‑147 (MQ=255)
                       accaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                           <  2:3586717/149‑1 (MQ=255)
                       accaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                           <  2:3402718/149‑1 (MQ=255)
                        ccaTCCAGTGCAGGAGCTCCAATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCTTCAAAAACCCTCCCCGCGGATTTGTACTaa                                                                                                                          >  1:3554217/1‑148 (MQ=25)
                           aaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCGCCCGGCGGAGTTTGCCTACTCa                                                                                                                       <  2:13652‑M1/96‑1 (MQ=255)
                            aaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAg                                                                                                                      <  2:1526748‑M1/97‑1 (MQ=255)
                             aggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAgg                                                                                                                     >  2:744367‑M1/52‑149 (MQ=255)
                             aggcccagtctttcgactgagccgttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAgg                                                                                                                     <  2:3038213‑M1/98‑1 (MQ=255)
                              ggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattTATTCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATGTGGCCTTCTCaggc                                                                                                                    <  2:3363458‑M1/99‑2 (MQ=255)
                              aGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTGCTACTCAggg                                                                                                                    >  1:3036980/1‑148 (MQ=255)
                               gTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCACTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                                   <  2:2034450/149‑1 (MQ=255)
                                 ccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGGAGTTTATGGCGGGCGTCCTGCCCGCCACCCACCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCAACTCAGGagag                                                                                                                 <  2:3307116‑M1/102‑1 (MQ=255)
                                 ccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACGTTCAAATCCGCTCCCCGCGGAGGTGTCCTCCTCAGGagag                                                                                                                 <  2:1673134‑M1/102‑1 (MQ=255)
                                 ccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTTCGCAACGGTCAAATCCGCTCCCGGCGGAGTTGTCCTACGAAGGagag                                                                                                                 <  2:2358573‑M1/102‑1 (MQ=255)
                                   agtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAAACCGCACCCCGCGGATTTGTACTACTCAGGAGAGCg                                                                                                               <  1:4230194‑M1/104‑1 (MQ=255)
                                   agtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAAAGTTAAAAAACGCACCCCGCGGCGTTTTACTACTCCTGAGAGCg                                                                                                               <  1:4440601‑M1/104‑1 (MQ=255)
                                    gtctttcgactgagcctttcgttttatttgatgcctggggaattAATACCTGGAACTTAATGGCGTGCGTCATGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGGTCAAATCCGCTACCGGCGCATTTGTCCTACTCAGGCGAGCGt                                                                                                              <  2:4008466‑M1/105‑1 (MQ=255)
                                     tctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                             <  2:284229‑M1/106‑1 (MQ=255)
                                     tctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGtt                                                                                                             <  2:4134762‑M1/106‑1 (MQ=255)
                                      aGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                                            <  2:2521299/149‑1 (MQ=255)
                                       gCTCCTATTATAAGCTGTCAAACATGAGAATTCCGATCCCCATTTCCTGGCTTTTTTGGGGGGGGGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                                           <  1:1286688/149‑1 (MQ=18)
                                        ttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAGCCGCTCCCGGCGGATTTGTCCTACTCAGGAGGGCGttctc                                                                                                          <  2:2574891‑M1/109‑3 (MQ=255)
                                         tCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                                         <  2:4519084/149‑1 (MQ=255)
                                         tCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                                         <  2:2720598/149‑1 (MQ=255)
                                         tCCGATGATAAGCTGTCAAACATGAGAATTCCGAGCCCCACTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                                         <  1:2330879/149‑1 (MQ=255)
                                            gatgatATGCTGTCAAACATTTGAATTCCTATCCCCAATTCCTGGTATTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGac                                                                                                      <  1:663504/149‑1 (MQ=17)
                                             ctgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGaca                                                                                                     <  2:2386402‑M1/114‑1 (MQ=255)
                                              tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTAGCAACGTTAAAAACCCCTCCCGGCGGATTTTTACTAATCAGGAGAGCGTTCACCGacca                                                                                                    >  1:615143/1‑147 (MQ=18)
                                                agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTCAAAAACCCACCCCGCGGATTTGTCCAACTCAGGCGAGCGGTCACCCacaaac                                                                                                  <  1:2481812‑M1/117‑1 (MQ=255)
                                                 tAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACCGTCAAAACCGCACCCGGCGGATTTGTCCAACTCAGGAGCGCGGTCACCCccaaaca                                                                                                 >  1:2503730/1‑149 (MQ=18)
                                                  cctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaacaa                                                                                                <  2:3381086‑M1/119‑1 (MQ=255)
                                                   aGCTGTCAAACATGAGAATTCCGATCCCCAATTCCGGGCAGTTTAGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAaacaac                                                                                               <  1:3357289/149‑1 (MQ=255)
                                                    gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                                              <  2:312567/149‑1 (MQ=255)
                                                      tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTAAAAACAGCTCCCGGCGGATTTGTCCTACACAGGAGAGCGTTAACCGAAAAacaacaca                                                                                            >  1:4137603/1‑147 (MQ=17)
                                                          tttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                                        <  2:4329319‑M1/127‑1 (MQ=255)
                                                          tttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAACGTTCAAATCCGCACCCGGCGGATTTGTCCTACACAGGAGAGCGTTCACCCACAAACAACAGAaaaa                                                                                        <  2:4384619‑M1/127‑1 (MQ=255)
                                                          aaaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                                        >  2:2643781/1‑149 (MQ=255)
                                                          aaaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                                        <  1:4561698/149‑1 (MQ=255)
                                                          aaaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                                        >  2:1956082/1‑149 (MQ=255)
                                                           aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTGCGACTCAGGAGAGCGTTCACCGACAAACAACAGATaaat                                                                                       >  2:1796009/1‑148 (MQ=255)
                                                            tatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                                      <  2:606257‑M1/129‑1 (MQ=255)
                                                              tttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCACCCGGCGGATTTGTCCTACTCAGGCGTGCGTTCACCGACAAAAAAAAGATAAAACGa                                                                                    <  2:1205796‑M1/131‑1 (MQ=255)
                                                               tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                                                                                   <  2:376302/149‑1 (MQ=255)
                                                                  gaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGATTTTTCCTACTCCGGAGAGCGTTCACCGACAAACCACAGATAAAACCAAAgg                                                                                >  1:2674830/1‑149 (MQ=14)
                                                                  atgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGTGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                                                                                <  2:2285740‑M1/135‑1 (MQ=255)
                                                                    tgctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                                              >  2:2674830‑M1/13‑149 (MQ=255)
                                                                    aTTCCGATCCCCATTTCCTGGCTGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACTTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                                              <  1:2447199/149‑1 (MQ=255)
                                                                        ggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                                  <  2:3875270‑M1/117‑1 (MQ=255)
                                                                        ggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGt                                                                          <  2:4473563‑M1/141‑1 (MQ=255)
                                                                         ggtaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgaga                                                                                                                  >  1:1127501‑M1/8‑108 (MQ=255)
                                                                         gATCCCCAATTCCTGGCAGTTTTTGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc                                                                         <  1:2903584/149‑1 (MQ=255)
                                                                         gATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc                                                                         <  2:1737802/149‑1 (MQ=38)
                                                                          ggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTGCACCGACAAACAACAGCTAAAACGAAAGGCCCagtgt                                                                        <  2:465002‑M1/143‑3 (MQ=255)
                                                                            ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCttt                                                                      <  1:4612905/149‑1 (MQ=35)
                                                                            ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCttt                                                                      <  2:4569188/149‑1 (MQ=35)
                                                                             cccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                     <  1:1843847/149‑1 (MQ=35)
                                                                               tAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                                   >  1:1033106‑M1/2‑149 (MQ=255)
                                                                                acttCCTGGCAGTTTTTGGCGGGGGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                                  <  1:3912533/147‑1 (MQ=9)
                                                                                 tttCCTGGATGTTTCTGGGGGGGGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                                 <  1:2084385/148‑1 (MQ=11)
                                                                                 aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                                    <  2:350203/146‑1 (MQ=32)
                                                                                  ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                                <  1:2782587/149‑1 (MQ=31)
                                                                                   tCCTGGCTGTTTCGGGCGGGCGTCCGCCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAACCAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                               <  1:1673134/149‑1 (MQ=11)
                                                                                                          ccggCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAGCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCt                                        <  1:3671135/146‑1 (MQ=1)
                                                                                                              cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                                    <  2:3203026/149‑1 (MQ=31)
                                                                                                                 gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                 <  1:2010691/149‑1 (MQ=33)
                                                                                                                 gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                 >  2:1129372/1‑149 (MQ=33)
                                                                                                                 gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                 <  1:1419457/149‑1 (MQ=33)
                                                                                                                   cACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCGGgggacgtg                               >  2:3856579/1‑145 (MQ=17)
                                                                                                                           gggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGGCTGGGGAATTGGGaattaa                       >  2:1055080/1‑149 (MQ=255)
                                                                                                                                gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGCTAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCtt                  <  1:3307116/149‑1 (MQ=255)
                                                                                                                                gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCtt                  >  2:4207838/1‑149 (MQ=255)
                                                                                                                                gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAAGTGGGAATTAATTCtt                  >  2:4170970/1‑149 (MQ=255)
                                                                                                                                  tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGa                <  1:1169269/149‑1 (MQ=255)
                                                                                                                                     ttCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGa             <  1:4665265/149‑1 (MQ=255)
                                                                                                                                     gtCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGa             <  1:2324116/148‑1 (MQ=255)
                                                                                                                                          aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGaaa        >  2:3579005/1‑149 (MQ=255)
                                                                                                                                          aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTGCGTTTTATTTGCTGCCGGGGGAATTGGGAATTAATTCTTGAAGACGaaa        >  2:3057102/1‑149 (MQ=255)
                                                                                                                                                ctaaTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCt  <  1:3701265/147‑1 (MQ=255)
                                                                                                                                                cAACTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCt  <  1:3856579/149‑1 (MQ=255)
                                                                                                                                                cAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCt  <  2:2481812/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCT  >  CP009273/360024‑360316

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