Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,173 C→A 7.0% 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.0% 950.0 / 9.0 329intergenic (‑107/+488)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (81/207);  new base A (26/14);  total (107/222)
Fisher's exact test for biased strand distribution p-value = 8.30e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 333 total aligned reads displayed.

CCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGAT  >  CP009273/360024‑360322
                                                                                                                                                     |                                                                                                                                                     
gttcaccgacaaaccacagataaaaccataggcccttactttctaatgagactttgttattattttttgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATc                                                                                                                                                       >  2:5239167‑M1/81‑150 (MQ=255)
   caccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGTTCGCAACGTTCAAATCCGc                                                                                                                                                    <  2:1746705‑M1/73‑1 (MQ=255)
     ggCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGTTCAAATCCGCTc                                                                                                                                                  >  1:781449/1‑150 (MQ=255)
     ggCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                                                                  >  2:4816160/1‑150 (MQ=255)
       gacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                                                <  1:1553132‑M1/77‑1 (MQ=255)
        ttACCACCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGATTCCCGATCCCCAATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                                                                                               <  2:1338125/150‑1 (MQ=255)
           ccATCCCGCGCCACCATCCAGTGCTGGAGCTCCGATTATACGCTGTCAAACTTGATAATTCCGACCCCCATTTCCTGGATTTTTTGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcg                                                                                                                                            <  1:2033147/150‑1 (MQ=21)
           aacaacagataaaacgaaaggcccagtcttgcgactgagcctttcgtgttctttgatgccgggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCcggcg                                                                                                                                            <  2:1429436‑M1/81‑1 (MQ=255)
            cATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCGCACCccgcgg                                                                                                                                           >  1:2355952/1‑150 (MQ=255)
             aTCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                                          <  2:675851/150‑1 (MQ=255)
                 agataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTg                                                                                                                                      <  2:2612360‑M1/87‑1 (MQ=255)
                 agataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCCGGCGGGTTTg                                                                                                                                      <  2:1098711‑M1/87‑1 (MQ=255)
                 aGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTg                                                                                                                                      >  2:1558769/1‑150 (MQ=255)
                      cacaaTCCTGTCCTGTAGCTCCTATGATACGCTGAAAAACATCAGAATCCCGACCCCCAATCCCTGGCTGTTTTTGGCGGGCGCCCGGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                                 <  1:3141222/150‑1 (MQ=17)
                       aacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGGCAACGTTCGAATACGCGACCGGCGGctttgtgctca                                                                                                                                <  1:224264‑M1/93‑6 (MQ=255)
                         cgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCGCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATGTGTCCAACTc                                                                                                                              <  2:2136969‑M1/95‑1 (MQ=255)
                         caTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTCGCAACGTTCAAAAACGCGCCCCGCGGATGTTTCCGACTc                                                                                                                              >  1:3333157/1‑150 (MQ=25)
                         caTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTAGCAACCGTCAAATCCCCTCACCCCCGAATTGTGCTAcgc                                                                                                                              >  1:4420169/1‑148 (MQ=25)
                          gaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCctacgca                                                                                                                             <  2:2345127‑M1/96‑4 (MQ=255)
                           aaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTGAAATCCGCTCCCGGCGGGGTTGTGCTActccg                                                                                                                            <  2:3642261‑M1/97‑3 (MQ=255)
                            ccAGTGCTGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAgg                                                                                                                           <  2:4976624/150‑1 (MQ=255)
                            aaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAgg                                                                                                                           <  2:4271903‑M1/98‑1 (MQ=255)
                            aaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCCGGGAGCGGATTTGTCCTACTCAgg                                                                                                                           <  2:403721‑M1/98‑1 (MQ=255)
                                 gCAGGAGCTCCGATGATAAGCGGGCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                                      <  2:4742226/150‑1 (MQ=255)
                                    ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                                   <  1:1732570/150‑1 (MQ=255)
                                      ctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCACCCGGCGGAATTGTCCTACTCAGGAGAGcgttag                                                                                                                 <  2:1655912‑M1/108‑3 (MQ=255)
                                        ttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                                               <  2:4171100‑M1/110‑1 (MQ=255)
                                        cTCCGATGATAAGCTGTCACACATGAGAAGTCCGATCCCCAATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                                               <  2:2909324/150‑1 (MQ=255)
                                        cTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAAAGTTAAAATCCGCTGCCGGCGGGTTTGGGCTACTCAGGGGCGCGGGCAcc                                                                                                               >  1:4195606/1‑150 (MQ=18)
                                             ctgagcctttcgattaatttgatcccgggggaattAATGCCTGGCTGTTTTGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaa                                                                                                          >  1:5093768‑M1/36‑150 (MQ=255)
                                             atgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACCTTAAAAACCGCTCCCCGCGGATTTTTGCTAATAAGGAGAGCGTTAACCGacaa                                                                                                          >  1:4619678/1‑150 (MQ=18)
                                              tgatAAGCTGTAAAACATGAGAATTCCGATCCCCAATTCCTGGATGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaa                                                                                                         <  1:3752558/150‑1 (MQ=255)
                                              tgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCGACTCAGGAGAGCGTTCACCGacaaa                                                                                                         <  2:284198‑M1/116‑1 (MQ=255)
                                               gatAAGCTGTCAAACATGAGACTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                                        <  1:5232061/150‑1 (MQ=255)
                                                atAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTGGCCAAGTTGAAATACCCTCCCCGCGGGTTTGTGCTACACAGGTGGGCGTTCACCCccaaaaa                                                                                                       >  1:122191/1‑148 (MQ=18)
                                                 gcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCGCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTGAAATCCGCTCCCGGCGGATTTGTCCGACTCCGGCGAGCGTTcaccgacacaccg                                                                                                      <  2:2408169‑M1/119‑6 (MQ=255)
                                                  cctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAGTTGTCCTACTCAGGAGAGCGTTCACCGACAaacaac                                                                                                     <  2:4955026‑M1/120‑1 (MQ=255)
                                                  cctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAAACCCACCCCGCGGATTTGTCCTACACAGGAGAGCGTTCACCCACAaacaac                                                                                                     <  2:618152‑M1/120‑1 (MQ=255)
                                                  aaGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAAAGTTCAAATACGCGCCCGGCGGAGTTTTCCTACTAAGGAGGGGGGTCACCGccaaaaaaa                                                                                                     >  1:1607351/1‑149 (MQ=17)
                                                   aGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCCAACTTAAAATACCCGCCCCGCGGGGTTGTGCTACTCACGTGCGCGGTCACCCGCAAacaaca                                                                                                    >  1:4051099/1‑150 (MQ=17)
                                                    gCTGTCAAACATGAGAATTCCTATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                                   <  2:1891538/150‑1 (MQ=255)
                                                     cTGTCAAACATGAGAATTCCGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGa                                                                                                  <  2:5106407/150‑1 (MQ=255)
                                                      tGTCAAACATGAGATTTCCGATCCCCATTTCCTGGAAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                                                 <  1:2504017/150‑1 (MQ=255)
                                                      tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTTCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                                                 <  2:3996065/150‑1 (MQ=255)
                                                      tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                                                 <  2:858016/150‑1 (MQ=255)
                                                       cgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACACAGGAGAGCGTTCACCGACAAACAACAGATa                                                                                                <  2:749975‑M1/125‑1 (MQ=255)
                                                        gttttatttgatgcgtggggaattAATGCGGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                                               >  2:558198‑M1/25‑150 (MQ=255)
                                                         ttttattttaggcctgggtaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                                              >  2:3131710‑M1/24‑150 (MQ=255)
                                                           aaCATGAGAATTCCGATCCCCAATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                                            <  2:3538496/150‑1 (MQ=255)
                                                            aCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                                           <  1:2575338/150‑1 (MQ=255)
                                                             atttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                                                <  2:2514091‑M1/109‑1 (MQ=255)
                                                                 gttgcgggggtatttAATTCGTGGCTGTTTAGGGCGGGGGTCCTCCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                                                                                      >  1:3662591‑M1/16‑150 (MQ=255)
                                                                 atatttCCGATCCCCAATTCGTGGCTGTTTATGGCGGGGGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                                                                                      <  1:169635/146‑1 (MQ=255)
                                                                   aagtCCGAGCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                                                    <  2:2442886/147‑1 (MQ=255)
                                                                    aTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTAGCCACGTTCAAAAACGCTCCCGGCGGCTTTGTCCGACTCCGGAGAGCGGTCACCCACAAACAACAGAAAAAACGAAAGGccc                                                                                   >  1:4257784/1‑150 (MQ=14)
                                                                    aTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                                                                                   <  2:128378/150‑1 (MQ=255)
                                                                    aTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCTTCAAAAACGCGCCCCGCGGCTTTGTGCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAAGTAAGGccc                                                                                   >  1:299674/1‑150 (MQ=14)
                                                                     ttCCGATCCCCATTTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                                                  <  2:2214305/150‑1 (MQ=255)
                                                                     ttCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                                                  <  2:1957165/150‑1 (MQ=255)
                                                                     gtCCGATCCCCACTGCCTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAGAAAACGAAAGGCCCa                                                                                  <  2:1599061/149‑1 (MQ=11)
                                                                     cctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                                                  <  2:3449936‑M1/139‑1 (MQ=255)
                                                                       ccGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTAGCAAAGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCACGAGAGCGTTCACCGACCAACCACAGATAAAACGAAAGGCCCCGt                                                                                >  2:4909713/1‑150 (MQ=14)
                                                                        ggctcattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc                                                                               >  2:3715924‑M1/9‑150 (MQ=255)
                                                                         gATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                              >  2:1259335/1‑150 (MQ=38)
                                                                          ggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTGCTACTCAGGAGAGCGTTCACCCACAAACAACACATAAAAAGAAAGGCCCCgtcat                                                                             <  1:3432400‑M1/144‑3 (MQ=255)
                                                                          ggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                             <  2:3193745‑M1/144‑1 (MQ=255)
                                                                          aTCCCCTATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                             <  2:5073106/150‑1 (MQ=25)
                                                                            aattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAGTTGTCCGACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                           <  2:3023519‑M1/146‑1 (MQ=255)
                                                                             cccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                                          <  1:4702437/150‑1 (MQ=35)
                                                                              ttAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCCGCGGATTTGTACTACTCAGGAGAGCGTTCACCGACAAACAACACAAAAAAAGAAAGGCCCAGTCATTCGa                                                                         <  2:1495017‑M1/148‑1 (MQ=255)
                                                                              ccAATTCCTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                                         <  1:2934422/150‑1 (MQ=21)
                                                                               tAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                                        <  2:3497273‑M1/149‑1 (MQ=255)
                                                                                 aTTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                                      <  1:2027226/150‑1 (MQ=17)
                                                                                  ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                                     <  2:504813/150‑1 (MQ=31)
                                                                                  ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                                     <  2:3915092/150‑1 (MQ=31)
                                                                                   tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                                                  <  1:2100425/104‑1 (MQ=11)
                                                                                   tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                                                  >  2:2100425/1‑104 (MQ=11)
                                                                                   tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAg                                                                    <  1:1259335/150‑1 (MQ=11)
                                                                                             tttgggggggctTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGttt                                                          <  1:1031007/143‑1 (MQ=0)
                                                                                                            tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                                           <  1:3635476/150‑1 (MQ=11)
                                                                                                            tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTAAAATACGCTCCCGGCGGATTTTTCCTACTCAGGAGAGCGTTAACCGACAAACAACAGATAAAAAGAAAGGCCCAGTCTTTCGACTGAGCCGTTGGTTTTTTTTGTTGCCTggg                                           >  2:2719511/1‑150 (MQ=11)
                                                                                                             gCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAAGCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACACCAGAGAACACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                                          <  2:3367298/150‑1 (MQ=11)
                                                                                                              cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGggga                                         <  2:1524688/150‑1 (MQ=32)
                                                                                                               ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATGTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                                        <  2:443482/150‑1 (MQ=17)
                                                                                                               ccGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCAAATCCGCTCCCGGCGGATTTGTCCTACGCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTGTCGGTGTATTTGATGCCTGgggac                                        >  2:2620380/1‑149 (MQ=11)
                                                                                                                cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                       <  1:1995109/150‑1 (MQ=33)
                                                                                                                 gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTAACCGACAAACAAAAGATAAAACGAAATGCCCAGTCTTTAGACTGAGCCTTTCGTTTTATTTGATTCCTGgggactt                                      >  2:4773578/1‑150 (MQ=11)
                                                                                                                      ccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgggaa                                 <  1:957228/150‑1 (MQ=38)
                                                                                                                       cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACACCAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgggaat                                <  1:5004096/150‑1 (MQ=255)
                                                                                                                          cGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaa                             <  1:1487453/150‑1 (MQ=255)
                                                                                                                             gCCGTTGCTTCCCAACGTTCCACACCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCCTACACCGACAAACAACAGATAAAACGAACGGCCCAGACTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTc                          <  1:4820164/150‑1 (MQ=14)
                                                                                                                               cGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCtt                        >  2:4111254/1‑150 (MQ=255)
                                                                                                                               cGTTGCTTCGCAACGGTCTTATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCtt                        <  1:3175754/150‑1 (MQ=255)
                                                                                                                                 ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGa                      >  1:1176388/1‑150 (MQ=255)
                                                                                                                                  tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGaa                     >  1:1809320/1‑150 (MQ=255)
                                                                                                                                     ttCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGAc                  >  2:1986998/1‑150 (MQ=255)
                                                                                                                                     ttCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGACTTTCGACTGAGCCTTTCGTGTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGAc                  >  2:4615836/1‑150 (MQ=255)
                                                                                                                                          aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTGAAGACGAAAg             >  1:710644/1‑150 (MQ=255)
                                                                                                                                                   aTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTg    >  2:4240215/1‑150 (MQ=255)
                                                                                                                                                    tCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGa   >  1:2813943/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGAT  >  CP009273/360024‑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: