Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,173 C→A 7.6% intergenic (‑107/+488) 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,1730CA 7.6% 808.4 / 5.9 303intergenic (‑107/+488)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (71/193);  new base A (13/26);  total (84/219)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 313 total aligned reads displayed.

CAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGAT  >  CP009273/360025‑360322
                                                                                                                                                    |                                                                                                                                                     
cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                                      <  2:1477563/150‑1 (MQ=255)
cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                                      >  2:2302672/1‑150 (MQ=255)
cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                                      <  2:504651/150‑1 (MQ=255)
 tgCGGCTTACCATCCCGCGCAACCATCCTGTGCTGGCGCTCCGATTATACGGTGTCAAACTTGTTATTTCCGATCCCAATTTCCTGGATTTTTTGGGCGGGGGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                                                                     <  1:2191665/149‑1 (MQ=16)
 aGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCg                                                                                                                                                     >  2:4009583/1‑150 (MQ=255)
   cggataaccaccCCGCCCCACCATCCCGTGCTGGAGCTCCTTTTATACGCTTACAAACATGAGAATTCCGATCCCAATTTCCTGGCTGTTTTGGGGGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                                                   <  1:1017592/144‑1 (MQ=16)
   cGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGTTAAAAAACGCt                                                                                                                                                   >  1:1366841/1‑150 (MQ=255)
      cTTACCATCCCGCGCCACCATCCAGTGCAGGAGCTCCTATTATATGGTGTAAACCAGGACAATTCCGATCCCCAATTCCTGGCTGTTTTTGGCGGGCGCCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                                                <  1:4136734/150‑1 (MQ=21)
      cTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                                                <  2:1884180/150‑1 (MQ=255)
      cTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCCACGGTCAAAAACGCGccc                                                                                                                                                >  1:2923325/1‑150 (MQ=37)
       ttACCATCCCGCGCCACCATCCCGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                                                                                               <  1:1150123/150‑1 (MQ=255)
       ttACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTTCAAATCCGCTCCcg                                                                                                                                               >  1:3523605/1‑150 (MQ=255)
       ttACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                                                                                               <  2:1784126/150‑1 (MQ=255)
       acaaacaccagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                                                                                               >  2:1401996‑M1/73‑150 (MQ=255)
       acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTAAAAACCGCTCCcg                                                                                                                                               <  1:381131‑M1/78‑1 (MQ=255)
         aaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc                                                                                                                                             <  1:3934100‑M1/80‑1 (MQ=255)
          ccATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcg                                                                                                                                            >  2:3689504/1‑150 (MQ=255)
           cATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGTTCGCAACGTTCAAATACGCAGCccgcgg                                                                                                                                           >  1:889113/1‑150 (MQ=37)
           acaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTCAAATACGCTCCcggcgg                                                                                                                                           <  1:397113‑M1/82‑1 (MQ=255)
             tCCCGCGCCACCATCCCGGGCTGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                                         <  2:2028457/150‑1 (MQ=25)
                  cgcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                                    <  2:2170952/150‑1 (MQ=255)
                   gcaaccaTCCAGTGCTGGAGCTCCTTTGATAAGCTTAAAAACATGATAATTCCGATCCACAATTCCTGGCTTTTTAGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                                                                                   <  1:322433/147‑1 (MQ=18)
                     caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                                 >  2:3945301/1‑150 (MQ=255)
                         gaaaggcccagtctttcgactgtgcctttcgttttttttgaggcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                                             >  2:3144611‑M1/55‑150 (MQ=255)
                           ccAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAgg                                                                                                                           <  2:4635622/150‑1 (MQ=255)
                             aGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                                         <  2:2934346/150‑1 (MQ=255)
                             aGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                                         <  2:546274/150‑1 (MQ=255)
                              gTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATGCCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgaga                                                                                                                        <  2:4408496/150‑1 (MQ=255)
                               cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                                                                                       <  2:4697391‑M1/102‑1 (MQ=255)
                                  tgtccgtcgacggagcctttcgttttatttgctgcccggggcaatAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                                    >  2:4380229‑M1/46‑150 (MQ=255)
                                  agtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAAAGTTAAAATACCCACCCCGCGGATTTGTACTACTAAGGAGagccg                                                                                                                    <  2:1213657‑M1/105‑3 (MQ=255)
                                   ggcggtcgcctgagccgttcgatttatttgcttgcgggggacttCATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                                   >  2:3000465‑M1/45‑150 (MQ=255)
                                    gAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                                                  <  2:3922301/150‑1 (MQ=255)
                                    gAGCTCCGATGAGAAGCTGTCAAACATGAGAATTCCGACCCCCACTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                                                  <  2:3739667/150‑1 (MQ=255)
                                         ccGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGa                                                                                                             <  2:2729873/150‑1 (MQ=255)
                                          gactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGac                                                                                                            <  2:403784‑M1/113‑1 (MQ=255)
                                          gactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCGCTCCCCGCGGAATTGTCCAACTCAGGAGAGCGGTCACCGac                                                                                                            <  2:2441445‑M1/113‑1 (MQ=255)
                                           gatgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGaca                                                                                                           <  1:1113639/150‑1 (MQ=255)
                                            atgatAATCTGTCAAAAATGAGAATTCCTATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaa                                                                                                          <  2:4470655/150‑1 (MQ=255)
                                             tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaa                                                                                                         <  2:3309842/150‑1 (MQ=255)
                                             tgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAAACTTCAAATACGCACCCCGCGGATTTGTCCTACACACGAGAGcgttcaccaaaaaa                                                                                                         <  2:475382‑M1/116‑7 (MQ=255)
                                                gcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaacaa                                                                                                      <  2:3554466‑M1/119‑1 (MQ=255)
                                                  ctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAACGTTCAAATACGCTCCCCGCGGATTTTTCCAACAAAGGGGAGCGTTCACCCCCAAacaaca                                                                                                    <  1:4666575‑M1/121‑1 (MQ=255)
                                                  ctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTAGCAACGTTCAAATCCCCTCCCGGCGGATTTGTCCAACACAGGAGAGCGTTCACCGACAAaccaca                                                                                                    <  2:69668‑M1/121‑1 (MQ=255)
                                                   tttcgttttatttgatgcctggggaattAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                                   >  1:4372808‑M1/29‑150 (MQ=255)
                                                   gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                                   <  2:2510245/150‑1 (MQ=255)
                                                     tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                                                 <  1:3586796/150‑1 (MQ=255)
                                                      ggttttttttgatgcctggggatttATTTCCTGGCTGTTTAGGGGGGGCGTCGTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                                                >  2:1636401‑M1/26‑150 (MQ=255)
                                                      gTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                                                <  1:1991219/150‑1 (MQ=255)
                                                         tttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                                             <  2:3544033‑M1/128‑1 (MQ=255)
                                                         aaaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                                             <  2:341695/150‑1 (MQ=255)
                                                           aaaTGAGAATTCCGATCCCCATTTCCTGGCTGTTTTGGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                                           <  2:3625047/148‑1 (MQ=14)
                                                            cATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACGTTCAAAACCGCTCCCCGCGGAATTGTCCGACTCAGGAGAGCGGTCCCCCACCCACAACCGATAAAAccc                                                                                          >  2:421248/1‑148 (MQ=17)
                                                            atttgatgcctgggtaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                                                >  1:4254919‑M1/20‑128 (MQ=255)
                                                                agaATTCCGATCCCCAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                                                                                      <  1:4331309/150‑1 (MQ=255)
                                                                agaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                                                                                      <  1:1928470/150‑1 (MQ=255)
                                                                agaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                                                                                      <  2:2254259/150‑1 (MQ=255)
                                                                  tgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                                                    >  1:3175801‑M1/14‑150 (MQ=255)
                                                                  tgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAAACCCCTCCCCGCGGAGTTGTCCTACACAGGAGAGCGGTAACCCACAAACAACAAATAAAACGAAAGGcc                                                                                    <  1:749751‑M1/137‑1 (MQ=255)
                                                                   tttCCGATCCCCATTTCCGTGCTGTTTTGGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                                                                                   <  1:1748445/149‑1 (MQ=11)
                                                                   agtCCGATCCCCAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                                                                                   <  1:1973461/148‑1 (MQ=255)
                                                                   aTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAAACGCTACCGGCGGCTTTTTCCTACGAAGGAGAGCGTTGACCCACAAAAAAAAGATAAAAAGAAAGGgcc                                                                                   >  1:4770574/1‑150 (MQ=14)
                                                                    cctgggtaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                                             >  1:2335457‑M1/12‑139 (MQ=255)
                                                                    cctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAACGTTCAAATCCGCTCCCGGCGGATTTGTACTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                                                  <  2:3365854‑M1/139‑1 (MQ=255)
                                                                      tggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGt                                                                                <  2:2648218‑M1/141‑1 (MQ=255)
                                                                       ggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGCTTTGTCCTACGCAGGTGAGCGTTCACCGACACACAACAGATAAAACGAAAGGCCCCGTc                                                                               <  2:4481179‑M1/142‑1 (MQ=255)
                                                                       cGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc                                                                               <  2:1809671/150‑1 (MQ=39)
                                                                        gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCCGCGGCTTTGTCCTCCGCAGGGGCGCGGTCCCCGGCAAACGACAGATTAAACGTAAGGGCCagtaa                                                                              <  2:1297792‑M1/143‑3 (MQ=255)
                                                                        gATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                              <  2:1069449/150‑1 (MQ=38)
                                                                         aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCGCACCCGGCGGATTTGTCCTACGCCGGAGAGCCGTCACCCCCCAACAACAGAAAAAACCAAACGCCCACTcgg                                                                             >  1:3144611/1‑148 (MQ=11)
                                                                           ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                           <  1:662539/150‑1 (MQ=35)
                                                                           ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                           <  2:1156642/150‑1 (MQ=35)
                                                                           aattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                           <  2:1052443‑M1/146‑1 (MQ=255)
                                                                            cccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTTCTTCGCAACGTTCAAAACCCCTCCCGGCGGCTTTGTCCTACACAGGAGAGCGGTCACCGCCACAAAACACATAAAAAGAAATGCCCACTCTTTCg                                                                          >  1:2685095/1‑150 (MQ=11)
                                                                             ccAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                                         <  2:895076/150‑1 (MQ=17)
                                                                              cAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                                        <  2:655548/150‑1 (MQ=33)
                                                                               aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                                       <  2:3194550/150‑1 (MQ=32)
                                                                               aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTAACCGACAAACAACAGATAAAACGAAATGCCCACTCTTTCGACt                                                                       >  2:2078956/1‑150 (MQ=255)
                                                                                aTTCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                                      <  1:3049190/150‑1 (MQ=17)
                                                                                                      cgccccgcccgcccCCCTCCGGGCCGTTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAACCAACAGATAAACCGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGc                                                <  1:1621746/144‑1 (MQ=1)
                                                                                                           tGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTTCAAATCCGCTCCCGGCGGCTTTGTGCTACTCAGGAGGGCGTTAAACGACAAACAACAGATAAAAAGAAAGGCCCAGTCTTTCGACGGGGCCGTTGGGTTTATTTGTTGCCTggg                                           >  2:2864862/1‑150 (MQ=11)
                                                                                                           tGCCCGCCACCCTCCGGGCCGGTGCTTCGCAACGTTCAAAGCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGGTCACCGACAAACACCAGATAACACGAAAGGCCCCGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                                           <  2:1129548/150‑1 (MQ=11)
                                                                                                           ggCCCGCCACCCTCCGGGCCGGTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                                           <  1:2982170/149‑1 (MQ=11)
                                                                                                           gccccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                                          <  1:3262759/148‑1 (MQ=14)
                                                                                                            gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                                          <  2:1549759/150‑1 (MQ=31)
                                                                                                               cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                       <  1:459487/150‑1 (MQ=33)
                                                                                                                   aCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGCCAAACAACAGATAAAACGAAAGGCCCCGTCTTTCGCCTGAGCCTTTCGTTTTATTTGATGCCTGgggacttggg                                   >  2:3853336/1‑150 (MQ=255)
                                                                                                                      cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAAGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCGGgggaattgggaat                                >  2:1528505/1‑150 (MQ=255)
                                                                                                                        ccGGGCCGTTGCTTCGAAACGTTAAAATCAGCTCCCGGCGGATTTGTCCTACTCAGTAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaatta                              <  1:403728/150‑1 (MQ=255)
                                                                                                                             ccGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTGTCGTTTGATTTGATGCCTGGGGAAGTGGGAATTAATTCt                         >  1:2028464/1‑150 (MQ=255)
                                                                                                                                ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTGa                      >  2:1931769/1‑150 (MQ=255)
                                                                                                                                         aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTCTTTGATGCCGGGGGACTTGGGACTGAATTCTTGAAGACGAAAg             >  1:535322/1‑150 (MQ=14)
                                                                                                                                          aCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAgg            >  2:4321362/1‑150 (MQ=255)
                                                                                                                                              tCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCt        <  2:2365883/150‑1 (MQ=255)
                                                                                                                                              tCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCt        >  1:3271633/1‑150 (MQ=255)
                                                                                                                                                 caTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGt     <  1:1913330/149‑1 (MQ=255)
                                                                                                                                                  aTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTAGTTGATGCCTGGGGAAGTGGGACTTAATTCTTGAAGACGAAAGGGCCTCGTg    >  2:972383/1‑150 (MQ=255)
                                                                                                                                                    ccGCTCCCGGGGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAAAAACATATAAAACTAAAGGCCCTGTCTTTCGACTGAGCCTTTCGTTTTTTTTGATGCCTGGGGAATTGGGATTTAATTCTTGAAGACGAAAGGGCCTCGTGAt  <  1:2441445/150‑1 (MQ=18)
                                                                                                                                                    ccGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGGCTGGGGCATTGGGAATTAATTCTTGAAGACGAAAGGGCCGCGTGAt  >  1:2809959/1‑150 (MQ=255)
                                                                                                                                                    ccGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGAt  >  2:4553765/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGAT  >  CP009273/360025‑360322

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: