Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,181 G→C 8.6% 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 8.6% 721.0 / 11.7 234intergenic (‑115/+480)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (68/143);  new base C (13/10);  total (81/153)
Fisher's exact test for biased strand distribution p-value = 3.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 242 total aligned reads displayed.

TTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTAT  >  CP009273/360032‑360330
                                                                                                                                                     |                                                                                                                                                     
acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGTTCCAATCCGCACCcg                                                                                                                                                       <  2:2555649‑M1/78‑1 (MQ=255)
 tACCCTCCAGCGCCACCATCCTGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcgg                                                                                                                                                      <  1:3629873/150‑1 (MQ=255)
   aacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTTCTTCGCAACGTTCAAATCCGCTCCcggcg                                                                                                                                                    <  2:1868874‑M1/81‑1 (MQ=255)
     caacagataaaacgaaaggccccgtctttcgactgagcctttcgtgttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGGGCGTCGCACCGTTCAACTCCGCTCccggcggc                                                                                                                                                  <  2:3032383‑M1/83‑2 (MQ=255)
     aTCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                                                  <  2:4134265/150‑1 (MQ=255)
      tCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGATAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                                                 <  1:2816749/150‑1 (MQ=255)
      tCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGATTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                                                 <  1:4304120/150‑1 (MQ=255)
      tCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                                                 <  2:4211995/150‑1 (MQ=255)
           ataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                                            <  2:4261715‑M1/89‑1 (MQ=255)
             ccaccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCt                                                                                                                                          >  1:3585621/1‑150 (MQ=255)
              aaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                                         <  1:3891289‑M1/92‑1 (MQ=255)
               aacgaaaggcccagtctttcgactgagcctttcgttttatttgatgccgggggaattCATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTTCCGTTCCGGTCCCGGCGGAGTTGTCCTAc                                                                                                                                        <  2:4178909‑M1/93‑1 (MQ=255)
                acgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTAGCAACGTTCAAATAAGCTCCCGGCGGATttgtccaaca                                                                                                                                       <  2:1473315‑M1/94‑5 (MQ=255)
                 cgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTGAGCCACGTTCAAAACCGGTCCCGGCGGGTGTGTCCtacac                                                                                                                                      <  1:2717108‑M1/95‑3 (MQ=255)
                 caTCCAGTGCAGGCGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTCGCAACGTTAAAAAACGCGCCCCGCGGCTTTGTCCTAcac                                                                                                                                      >  1:2393267/1‑148 (MQ=21)
                  aTCCTGTGCTGGAGCTCCGATGATATGCTGACAAACATGAGAATTCCGACCCCCAATTCCTGGCAGTTTTGGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                                                     <  1:2318901/150‑1 (MQ=18)
                    ccAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCGCTACCGGCCGATTTGTCCTACTCacg                                                                                                                                   >  1:2453286/1‑148 (MQ=255)
                    aaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGTGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAgg                                                                                                                                   <  2:3978897‑M1/98‑1 (MQ=255)
                    aaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGTGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCGCTCCCGGCGGATTTGTACTACTCAgg                                                                                                                                   <  2:3337770‑M1/98‑1 (MQ=255)
                        tGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCCACCTTCAAATACCCAACCCGCGGATTTGTCCGACACACgagat                                                                                                                               >  1:4178292/1‑149 (MQ=21)
                        cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGCGTTGTCCTACTCAGGagag                                                                                                                               <  2:451555‑M1/102‑1 (MQ=255)
                        cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGGTCAATGCCGCTCCCGGCGGATTGGTCCTACTCAGGagag                                                                                                                               <  2:2624265‑M1/102‑1 (MQ=255)
                         ccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACCATAAAAAACGCTCCCCGCCGCTTTGTGCAACACAGGCGAGc                                                                                                                              <  1:4155751‑M1/103‑1 (MQ=255)
                          cAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCg                                                                                                                             >  2:3012410/1‑150 (MQ=255)
                           agtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                                            <  2:4260217‑M1/105‑1 (MQ=255)
                            ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                                           <  2:1777384/150‑1 (MQ=255)
                              aGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                                                         <  2:3813323/150‑1 (MQ=255)
                               gCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCTTTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                                                        <  1:2672244/150‑1 (MQ=255)
                                 tCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGTGCGTTCACCg                                                                                                                      >  2:2418664/1‑150 (MQ=255)
                                      tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTGGCAAACTTGAAATACGCAACCCGCGGATTTGTGCAAATCAAGGGAGCGGTGACCCacaaa                                                                                                                 >  2:370002/1‑150 (MQ=18)
                                      tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACGTTCAAATCCGCGCCCGGCGGCTTTGTCCGACGCAGGAGCGCGTTCACCGacaaa                                                                                                                 >  2:804647/1‑150 (MQ=18)
                                       tatAAGCTGTAAAACAGGACAATTCCTATCCCAAATACGTGGCATTTTATGGCGGGCGTCCGCCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                                                <  1:456309/149‑1 (MQ=14)
                                        agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaaca                                                                                                               <  2:2504106‑M1/118‑1 (MQ=255)
                                        agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAACTTAAAAAACGCTCCCGGCGGGTTTTTCCAACTCAGGAGAGCGTTCACccacaaaaa                                                                                                               <  1:942977‑M1/118‑3 (MQ=255)
                                         tAAGCTCACAAACATGAGAATACCGACCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaacaa                                                                                                              <  1:1951313/150‑1 (MQ=255)
                                          cctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGCAAaacaac                                                                                                             <  2:1570209‑M1/120‑1 (MQ=255)
                                          aaGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTCGCAACCTTCAAATACGCGCCCCGCGGATTTGTGCAACAAAGGAGAGCGGTAAACCacaaacaaa                                                                                                             >  1:3521524/1‑149 (MQ=17)
                                            gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                                           <  2:92621/150‑1 (MQ=255)
                                            gCTGTCAAACATGAGAAGTCCGATCCCCAATTCCTGGATGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                                           <  1:1674896/150‑1 (MQ=255)
                                             ttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCGCTCCCGGCGGATTTGTCCTACACAGGAGTGCGTTCACCGACAAACCACAGa                                                                                                          <  1:3995563‑M1/123‑1 (MQ=255)
                                                tCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                                                       <  2:3915038/150‑1 (MQ=255)
                                                 ttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCGACTCAGGAGAGCGTTCACCGACCAACCCCAGATaaa                                                                                                      <  2:3705473‑M1/127‑1 (MQ=255)
                                                 cAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTAAAATCCGCACCCCGCGGATTTGTCCTACTCAGGAGATCGTTAAACGACACACAACAGATaaa                                                                                                      >  2:3641361/1‑150 (MQ=17)
                                                  aaaaATGAGAATTCCTATCCCAAATTCCTGGCAGTTTTTGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                                                     <  2:2180778/150‑1 (MQ=14)
                                                   aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                                                    <  2:3485365/150‑1 (MQ=255)
                                                   aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTAGCAAACTTAAAAAACGCGACCCGCGGAGTTTTGCAACTAAGGAGAGCGTTAAACGACCAACAACACAAAAAAc                                                                                                    >  1:2253553/1‑150 (MQ=12)
                                                    tatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                                                   >  1:4260219‑M1/21‑150 (MQ=255)
                                                        tgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCTTCAAAAACGCTCCCGGCGGATTTGTCCTACACAGGAGAGCGGTCAACCCCCAACCACAGATAAaaagaaag                                                                                               <  2:764703‑M1/134‑1 (MQ=255)
                                                         gatgcctggggatttAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGAAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                                                                                              >  2:2953981‑M1/16‑150 (MQ=255)
                                                          atgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCCGCGGCTTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAAGAAAGGc                                                                                             <  2:3372534‑M1/136‑1 (MQ=255)
                                                           tgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTTCGCAACGGTGAAATACCCTCCCCGCGGCTTTGTGCTAATCAGGAGAGCGTTCACCCACAAACAAAAGATAAAAACAAAGGcc                                                                                            <  2:35423‑M1/137‑1 (MQ=255)
                                                             ttCCGATCCCCAATTCCTGGCAGTTTATGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                                                          <  2:564014/150‑1 (MQ=255)
                                                              ctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAg                                                                                         <  2:310959‑M1/140‑1 (MQ=255)
                                                               tggggaattAATCCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGCTTTGTCCTAATCAGGAGAGGGTTAACCGACAAACAACAGATAAAAAGAAAGGCCCAGt                                                                                        <  1:2463217‑M1/141‑1 (MQ=255)
                                                                ggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                                                        >  1:1051458‑M1/9‑133 (MQ=255)
                                                                ggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                                                        <  2:1051097‑M1/125‑1 (MQ=255)
                                                                 ggtaattAATTCGTGGCTGTTTTTGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGTTTTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                                      >  1:3337770‑M1/8‑150 (MQ=255)
                                                                 gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                                      >  1:1766945‑M1/8‑150 (MQ=255)
                                                                 gATCCCCAATTCCTGGCTGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                                      <  2:1105129/150‑1 (MQ=255)
                                                                  aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                                     <  2:2155885/150‑1 (MQ=38)
                                                                  aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                                     <  2:4136492/150‑1 (MQ=38)
                                                                  aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACGGTCAAAACCCCTCCCGGCGGCGTTGTCCCACGCCGGCGAGCGGTCCCCGACAAACCACAGATAAAACCAAAGGCCCCGTCtt                                                                                     >  2:3679593/1‑150 (MQ=11)
                                                                     cccAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                                                  <  1:3970826/150‑1 (MQ=25)
                                                                     attAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGCCAAACACCAGCTAAAACGCAAGGCCCAGTCTTTCg                                                                                  <  2:4199050‑M1/147‑1 (MQ=255)
                                                                      ttAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                                                 >  2:2497959‑M1/3‑150 (MQ=255)
                                                                      ccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                                                 <  2:1743458/150‑1 (MQ=34)
                                                                       cAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                                                <  1:2882237/150‑1 (MQ=33)
                                                                        aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGGCAAAGGTCAAAACCGCGCCCGGCGGGTTTGTCCTACACAGGAGGGCGTTCACCGCCCAACCACAGATAAAACGCAAGGGCCAGTCTTTCGTCt                                                                               >  1:1999055/1‑150 (MQ=11)
                                                                                      tagggcgggcgCCCGCCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACCCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGACAGGCCCAGTCTTTCGACTGAGCCTTTCGtttt                                                                 <  1:2958077/147‑1 (MQ=1)
                                                                                                     gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAGCCGCTCCCGGCGGATTGGTCCTCCTCAGGAGAGCGTTCACCGACAAACAACAGATAAACCGAAAGGCCCAGTCTTTCGACTGAGCCGTTCGTTTTATTTGATGCCTgggg                                                  <  2:2931970/150‑1 (MQ=11)
                                                                                                       ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                                                <  1:2718167/150‑1 (MQ=32)
                                                                                                        cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                               <  1:3012410/150‑1 (MQ=33)
                                                                                                         gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGACTGGTCCGACTCAGGAGAGCGTTCACCGACAAACAACAGATACAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                              <  2:2094991/150‑1 (MQ=11)
                                                                                                         gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                              <  1:3453626/150‑1 (MQ=34)
                                                                                                                         ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATGGGGAATTAATTCTTGa                              >  1:2884406/1‑150 (MQ=255)
                                                                                                                          tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattg                                             >  1:3339450/1‑134 (MQ=35)
                                                                                                                             ttCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTGAAGAc                          >  2:636387/1‑150 (MQ=255)
                                                                                                                              tCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACg                         <  2:1635946/150‑1 (MQ=255)
                                                                                                                               cGCAACGTTCAACTCCGCGCCCGGCGGATTTGTCCTACGCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGa                        <  1:136133/150‑1 (MQ=255)
                                                                                                                                gCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTGGTTTTATTTGATGCCTGGGGAAGTGGGAATTAATTCTTGTAGACGaa                       >  1:4094784/1‑150 (MQ=255)
                                                                                                                                 gaaCTTTAAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAACCAAAAGATAAAACGACAGGCCCAGTCTTTAGACTGAGCCTTTCGTTTTATTTGATGCCTGTGGAATTGGGAATTAATTCTTGAAGACGaaa                      <  1:1360970/149‑1 (MQ=17)
                                                                                                                                   aCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATACAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAgg                    <  2:1913019/150‑1 (MQ=255)
                                                                                                                                    cGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTGTCGTGTTATTTGATGCCTGGGGAAGTGGGAATTAATTCTTGAAGACGAAAggg                   >  1:599657/1‑150 (MQ=255)
                                                                                                                                        cAAACCCGCCCCCGGCGGATTGGTCCTACTCCGGAGAGCGTACACCGACAAACCACAGATAAAACGACAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGTAATTAATTCTTGAAGACGAAAGGGCCTc               <  1:764703/150‑1 (MQ=17)
                                                                                                                                          aaTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGa                        >  2:136131/1‑139 (MQ=255)
                                                                                                                                            tCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTGTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGa           >  1:1769771/1‑150 (MQ=255)
                                                                                                                                             ccGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATCAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATGGGGAATTAATTCTTGAAGACGAACGGGCCTCGTGAt          >  2:2708582/1‑150 (MQ=255)
                                                                                                                                             ccGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGAt          <  1:4285960/150‑1 (MQ=255)
                                                                                                                                              cGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATa         >  2:2791204/1‑150 (MQ=255)
                                                                                                                                              cGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATa         >  1:3141250/1‑150 (MQ=255)
                                                                                                                                               gCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATAc        >  1:1042341/1‑150 (MQ=255)
                                                                                                                                                cgcccGGCGGATTTGTCCTACTCCGGAGAGCGTTAACCGAAAAACAACAGATACAAGGAATGGCCCTGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTATTTCTTTAAGACGAAAGGGCCTCGTGATACg       <  1:2561761/148‑1 (MQ=18)
                                                                                                                                                cTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGACTTAATTCTTGAAGACGAAAGGGCCTCGTGATACg       >  1:878425/1‑150 (MQ=255)
                                                                                                                                                cTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTTATTCTTGAAGACGAAAGGGCCTCGTGATCCg       >  2:3789762/1‑150 (MQ=255)
                                                                                                                                                cTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACg       <  2:3593255/150‑1 (MQ=255)
                                                                                                                                                 tCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTATTTCTTGAAGACGATAGGGCCTCGTGATACGc      >  2:3146632/1‑150 (MQ=255)
                                                                                                                                                 tCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGc      >  2:86559/1‑150 (MQ=255)
                                                                                                                                                 tCCCGGCGGATTTGTCCAACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCGTTCGACAGAGCCTTTCGTGTTATTTGAAGCCTGGGGCATTGGGACGTAATTCTTGAAGACGAACTGGCCTCGTGATACGc      >  1:1271556/1‑150 (MQ=17)
                                                                                                                                                  cccGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGcc     >  2:2810593/1‑150 (MQ=255)
                                                                                                                                                     ggcggATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACCACAGATACAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTAt  <  2:2340252/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: