Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,181 G→C 6.8% intergenic (‑115/+480) 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,1810GC 6.8% 936.0 / 8.8 309intergenic (‑115/+480)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (87/194);  new base C (16/9);  total (103/206)
Fisher's exact test for biased strand distribution p-value = 1.51e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 332 total aligned reads displayed.

TTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTAT  >  CP009273/360032‑360330
                                                                                                                                                     |                                                                                                                                                     
ttACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                                                                                                       <  2:2108213/150‑1 (MQ=255)
ttACCACCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGATTCCCGATCCCCAATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                                                                                                       <  2:1338125/150‑1 (MQ=255)
   ccATCCCGCGCCACCATCCAGTGCTGGAGCTCCGATTATACGCTGTCAAACTTGATAATTCCGACCCCCATTTCCTGGATTTTTTGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcg                                                                                                                                                    <  1:2033147/150‑1 (MQ=21)
    acaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGACGTTGGTTAGCAAAGTTGAAATACGCTCCcggcgg                                                                                                                                                   <  1:406680‑M1/82‑1 (MQ=255)
     caacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCTAATACCCtcccggccgc                                                                                                                                                  <  2:4830681‑M1/83‑4 (MQ=255)
        cAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCACCGTTCAAATCCGCTCCCGGCGGAttt                                                                                                                                               <  2:2610752/150‑1 (MQ=255)
         agataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCCGGCGGGTTTg                                                                                                                                              <  2:1098711‑M1/87‑1 (MQ=255)
           cgcCACCATCCAGTGCTGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                                            <  2:1832871/150‑1 (MQ=255)
              caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                                         <  2:1965716/150‑1 (MQ=255)
               acaaTCCAGTGCAGGTGCTCCGATGATAAGCGGACAAACCTGAGAATTCCGATCCCCACTGCCTGGCTGTTTCTGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                                                        <  1:3262191/147‑1 (MQ=18)
                 caTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGGGTTGTCCTACTc                                                                                                                                      >  2:3720313/1‑150 (MQ=255)
                 caTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTAGCAACCGTCAAATCCCCTCACCCCCGAATTGTGCTAcgc                                                                                                                                      >  1:4420169/1‑148 (MQ=25)
                  gaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGCGTTGTCCTTCTCa                                                                                                                                     <  2:1468201‑M1/96‑1 (MQ=255)
                  aTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                                                     <  2:4326633/150‑1 (MQ=255)
                   aacggccctgtctttcgactgagcctttgtttttatttgttgccgggggaattAATCCCTGGCTTTTTAGGGGGGGCTTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAg                                                                                                                                    >  2:3072047‑M1/54‑150 (MQ=255)
                     aggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCCCGCGGAGTTGTCCAACTCAGga                                                                                                                                  <  2:49359‑M1/99‑1 (MQ=255)
                         ccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                                              <  2:5000475‑M1/103‑1 (MQ=255)
                           aGGAGCTCCGATGATACGCTGTCACACAGGAGAATTCCGATCCCCACTTCCGGGCAGTTTCGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTTCGCACCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGt                                                                                                                            <  2:2021302/150‑1 (MQ=18)
                            ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                                           <  1:1732570/150‑1 (MQ=255)
                              ctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCACCCGGCGGAATTGTCCTACTCAGGAGAGcgttag                                                                                                                         <  2:1655912‑M1/108‑3 (MQ=255)
                               tttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGt                                                                                                                                             >  1:391486‑M1/42‑129 (MQ=255)
                               tttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                                                        >  2:3490665‑M1/42‑150 (MQ=255)
                                ttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                                                       <  2:4171100‑M1/110‑1 (MQ=255)
                                cTCCGATGATAAGCTGTCACACATGAGAAGTCCGATCCCCAATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                                                       <  2:2909324/150‑1 (MQ=255)
                                   gactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGac                                                                                                                    >  2:1553132‑M1/38‑150 (MQ=255)
                                   gactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCAACACAGgagagcgttaacagaa                                                                                                                    <  2:4777783‑M1/113‑8 (MQ=255)
                                    tttgatAAGCTGTAAAACATGATAATACCGATCCCCAATTCCTGGAATTTTTTGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGaca                                                                                                                   <  1:3705224/148‑1 (MQ=17)
                                     ctgagcctttcgattaatttgatcccgggggaattAATGCCTGGCTGTTTTGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaa                                                                                                                  >  1:5093768‑M1/36‑150 (MQ=255)
                                     atgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACCTTAAAAACCGCTCCCCGCGGATTTTTGCTAATAAGGAGAGCGTTAACCGacaa                                                                                                                  >  1:4619678/1‑150 (MQ=18)
                                      tgatAAGCTGTCAAACATGATAATTCCGATCCACAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaa                                                                                                                 <  2:1512717/150‑1 (MQ=255)
                                      tgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGGTTTTTCCTACTCAGGAGAGCGGTAACCGacaaa                                                                                                                 <  2:283956‑M1/116‑1 (MQ=255)
                                       gatAAGCTGTCAAACATGAGACTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                                                <  1:5232061/150‑1 (MQ=255)
                                        atAAGTTGTCAAACTTGTGTATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaaca                                                                                                               <  2:2390107/150‑1 (MQ=255)
                                        atAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTGGCCAAGTTGAAATACCCTCCCCGCGGGTTTGTGCTACACAGGTGGGCGTTCACCCccaaaaa                                                                                                               >  1:122191/1‑148 (MQ=18)
                                        agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGGGCGTTCACCGacaaaca                                                                                                               <  2:351292‑M1/118‑1 (MQ=255)
                                         tAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaacaa                                                                                                              <  1:347998/150‑1 (MQ=255)
                                          cctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAaacaac                                                                                                             <  2:4957713‑M1/120‑1 (MQ=255)
                                           aGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                                                            <  2:1572099/150‑1 (MQ=255)
                                             cTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGa                                                                                                          <  2:2478463/150‑1 (MQ=255)
                                              tgtattaatttgatgccgggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                                                         >  2:4051099‑M1/27‑150 (MQ=255)
                                              tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTTCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                                                         <  2:3996065/150‑1 (MQ=255)
                                              tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                                                         <  2:858016/150‑1 (MQ=255)
                                                tCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGGTTCGCCACCGGCAAATACGCTACCCGGGGCGTTGTGCGTCTCGGGGGGGCGGTGACCGCCAACCAACAGATaa                                                                                                       >  1:2710280/1‑150 (MQ=12)
                                                gttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                                                       <  2:5048236‑M1/126‑1 (MQ=255)
                                                 cAAACATGAGAATTCCTATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                                                      <  2:1892567/150‑1 (MQ=255)
                                                 cAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                                                      <  2:227479/150‑1 (MQ=255)
                                                   aaCATGAGAATTCCGATCCCCAATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                                                    <  2:3538496/150‑1 (MQ=255)
                                                    aCATGAGAATTCCGATCCCCAAGTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCTGATTTGTCCTACTCAGGAGAGCGTTCACAGACAAACAACAGATAAAAct                                                                                                   <  2:679625/150‑2 (MQ=255)
                                                         atatttCCGATCCCCAATTCGTGGCTGTTTATGGCGGGGGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                                                                                              <  1:169635/146‑1 (MQ=255)
                                                           aagtCCGAGCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                                                            <  2:2442886/147‑1 (MQ=255)
                                                           aaTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                                                            <  2:647937/150‑1 (MQ=255)
                                                           aaTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAAAACGCTACCGGCGGAGTTTTACTAAACACGAGTGCGTTCAACGACAAACAACAGATAAATAGAAAGGcc                                                                                            >  1:2979311/1‑150 (MQ=14)
                                                            aTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTAGCCACGTTCAAAAACGCTCCCGGCGGCTTTGTCCGACTCCGGAGAGCGGTCACCCACAAACAACAGAAAAAACGAAAGGccc                                                                                           >  1:4257784/1‑150 (MQ=14)
                                                            aTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                                                                                           <  2:128378/150‑1 (MQ=255)
                                                            aTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCTTCAAATACCCTACCCGCGGATTTTTGCTACTCAGGAGAGCGTTCACCCACCAACAAAAGAAAAAAACAAATGccc                                                                                           >  1:2924290/1‑150 (MQ=14)
                                                                 gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                                      <  2:5240823‑M1/143‑1 (MQ=255)
                                                                 gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCCCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACCACAGATAAAAAGAAAGGCCCCGTCt                                                                                      <  2:4447683‑M1/143‑1 (MQ=255)
                                                                 gctcacaaaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                                      <  2:1741991/143‑1 (MQ=12)
                                                                 gATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                                      >  2:1259335/1‑150 (MQ=38)
                                                                  ggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                                     <  2:3193745‑M1/144‑1 (MQ=255)
                                                                  aTCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                                     <  2:465677/150‑1 (MQ=38)
                                                                    ccccAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                                   <  1:809345/150‑1 (MQ=25)
                                                                    ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                                   >  2:4764185/1‑150 (MQ=35)
                                                                     cccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                                                  <  1:4702437/150‑1 (MQ=35)
                                                                        aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCGTCAAATACCCTCCCGGCGGCTTTGTCCGACTCAGGAGAGCGTTAACCCACAAACAACACATAAAAACAAAGGCCCACTCATTAGAca                                                                               >  1:1664773/1‑149 (MQ=11)
                                                                         aTTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                                              <  1:2027226/150‑1 (MQ=17)
                                                                          ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                                             <  2:3915092/150‑1 (MQ=31)
                                                                          ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                                             <  2:507778/150‑1 (MQ=31)
                                                                          ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTCAAAAACGCTCCCGGCGGCTTTGTACTACTCAAGAGAGCGTTCACCGACAAACAACAGATAATAAGAAAGGCCCAGTGtttataaagg                                                                             >  2:2207769/1‑143 (MQ=11)
                                                                           tCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAg                                                                            <  1:316717/150‑1 (MQ=11)
                                                                           tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                                                          <  1:2100425/104‑1 (MQ=11)
                                                                           tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAg                                                                            <  1:1259335/150‑1 (MQ=11)
                                                                               ggaaGTTTAGGGGGGGCGTCCTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCtt                                                                        <  1:2249009/147‑1 (MQ=0)
                                                                                     tttgggggggctTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGttt                                                                  <  1:1031007/143‑1 (MQ=0)
                                                                                                    tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAAGAAAGGCCCAGTCTTTCGACTGAGCCATTCGTTTTTTTTGATGCCTggg                                                   >  2:2719508/1‑150 (MQ=2)
                                                                                                    tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTACCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                                                   <  1:3373703/150‑1 (MQ=11)
                                                                                                    tGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                                                   <  1:3635517/150‑1 (MQ=11)
                                                                                                     gCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAAGCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACACCAGAGAACACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                                                  <  2:3367298/150‑1 (MQ=11)
                                                                                                      cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAACCACCAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGggga                                                 <  2:1050448/150‑1 (MQ=12)
                                                                                                      cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGggga                                                 <  1:1570713/150‑1 (MQ=32)
                                                                                                       ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGACTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                                                <  2:2979311/150‑1 (MQ=17)
                                                                                                       ccGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCAAATCCGCTCCCGGCGGATTTGTCCTACGCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTGTCGGTGTATTTGATGCCTGgggac                                                >  2:2620380/1‑149 (MQ=11)
                                                                                                        cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                               <  1:1995109/150‑1 (MQ=33)
                                                                                                        cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAGTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTGTTATTTGATGCCTGgggact                                               >  2:364911/1‑148 (MQ=12)
                                                                                                               cTCCGTGCCGTTGCTTCTCAACGTTCAAATCCGCTCCCTGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgggaat                                        <  1:2499868/150‑1 (MQ=255)
                                                                                                               cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACACCAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgggaat                                        <  1:5004096/150‑1 (MQ=255)
                                                                                                                  cGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaa                                     <  1:1487453/150‑1 (MQ=255)
                                                                                                                     gCCGTTGCTTCCCAACGTTCCACACCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCCTACACCGACAAACAACAGATAAAACGAACGGCCCAGACTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTc                                  <  1:4820164/150‑1 (MQ=14)
                                                                                                                       cGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                               <  2:5117320/135‑1 (MQ=33)
                                                                                                                       cGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCtt                                >  2:573094/1‑150 (MQ=255)
                                                                                                                       cGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCtt                                >  2:4111254/1‑150 (MQ=255)
                                                                                                                        gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCGGGGGACTTGGGAATTAATTCTTg                               >  1:4222299/1‑150 (MQ=255)
                                                                                                                         ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGa                              >  1:1176388/1‑150 (MQ=255)
                                                                                                                          tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGaa                             >  1:1809320/1‑150 (MQ=255)
                                                                                                                               cGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgg                                            >  2:4255321/1‑130 (MQ=37)
                                                                                                                                  aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTGAAGACGAAAg                     >  1:710644/1‑150 (MQ=255)
                                                                                                                                  aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAg                     <  1:2181693/150‑1 (MQ=255)
                                                                                                                                           aTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTg            >  1:2531695/1‑150 (MQ=255)
                                                                                                                                            tCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGa           <  1:3026391/150‑1 (MQ=255)
                                                                                                                                                     ggcggATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTAt  <  2:3136994/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTAT  >  CP009273/360032‑360330

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: