Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,129 T→G 10.3% intergenic (‑63/+532) 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,1290TG10.3% 253.5 / 2.2 97intergenic (‑63/+532)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (46/35);  new base G (11/2);  total (59/38)
Fisher's exact test for biased strand distribution p-value = 7.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGA  >  CP009273/359980‑360255
                                                                                                                                                     |                                                                                                                                 
cTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGt                                                                                                                                   >  1:996790/1‑150 (MQ=255)
 cgttcaaatccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTc                                                                                                                                  >  1:81404‑M1/124‑150 (MQ=255)
  gttcaaatccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTcc                                                                                                                                 >  1:439180‑M1/123‑150 (MQ=255)
   tttACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTcc                                                                                                                                 <  1:266889/149‑1 (MQ=255)
   tttACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCt                                                                                                                                <  1:85718/150‑1 (MQ=255)
    tcaaatccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGCCCTg                                                                                                                               <  1:687100‑M1/30‑1 (MQ=255)
      aaatccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGc                                                                                                                              >  1:927152‑M1/119‑149 (MQ=255)
      aCCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGc                                                                                                                              <  1:664371/149‑1 (MQ=255)
       aatccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGccc                                                                                                                            <  2:1010426‑M1/33‑1 (MQ=255)
        cTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCg                                                                                                                           <  1:442899/150‑1 (MQ=255)
        atccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGGGGCCTgcccc                                                                                                                           <  2:452142‑M1/34‑2 (MQ=255)
          ccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGcc                                                                                                                             >  2:169687‑M1/115‑146 (MQ=255)
          ccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGGCCTGcc                                                                                                                             <  1:169687‑M1/32‑1 (MQ=255)
           cgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCa                                                                                                                        <  1:379871‑M1/37‑1 (MQ=255)
            tgGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCCGcccccccc                                                                                                                       >  1:186426/1‑148 (MQ=255)
            gctcccggcggatttgtcctactcaggagagcgttcaccgaaaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCAc                                                                                                                       >  2:419140‑M1/113‑150 (MQ=255)
                ccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCt                                                                                                                    >  1:368277‑M1/109‑149 (MQ=255)
                 cGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTcc                                                                                                                  <  1:1020585/150‑1 (MQ=255)
                  ggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCg                                                                                                                 >  2:98792‑M1/107‑150 (MQ=255)
                   aCATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCgg                                                                                                                <  1:724115/150‑1 (MQ=255)
                        cAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGt                                                                                                           <  1:101323/150‑1 (MQ=255)
                        cAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCCCCCGGcccg                                                                                                           >  2:912662/1‑148 (MQ=255)
                             gCACTTCACCGCTTGCCAGCGGCTTTCCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                      <  1:906831/150‑1 (MQ=255)
                             gCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCCGCCCTTGcct                                                                                                      >  2:423787/1‑148 (MQ=255)
                             gCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCCGCCGCTGCtt                                                                                                      >  2:665365/1‑150 (MQ=255)
                                cTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGTCCTGGCCCCCCCCCCCCCGGCCCGTGCTTccc                                                                                                   >  2:505395/1‑148 (MQ=255)
                                  tCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGGCCGCCACCCCCCCGGCCGCTGCTTCGCaa                                                                                                >  1:676633/1‑150 (MQ=255)
                                     ccGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGCTGCTTCGCAACGt                                                                                           >  1:126647/1‑150 (MQ=38)
                                      cGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                          <  1:581883/150‑1 (MQ=255)
                                         ttGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCCCTGCTTCGCCACGGTCa                                                                                        >  2:498679/1‑149 (MQ=38)
                                             ttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                   >  2:379871‑M1/80‑150 (MQ=255)
                                             ttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                   >  1:948069‑M1/80‑150 (MQ=255)
                                             ttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                   >  1:134017‑M1/80‑150 (MQ=255)
                                              tcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGGGTCTTGGCCCCCCCCCCCCCGGCCGTTGCTTCCCCAACgttcaatacc                                                                                   <  2:950180‑M1/72‑5 (MQ=255)
                                                 ggCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                               <  1:980988/150‑1 (MQ=255)
                                                   cTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCCCCCCCCCGGCCGTTGCTTCTCCACGTTCAAAACCGCcccc                                                                             >  2:5488/1‑150 (MQ=255)
                                                    acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                            >  1:264962‑M1/73‑150 (MQ=255)
                                                     caaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcgg                                                                           >  1:657522‑M1/72‑150 (MQ=255)
                                                      aCCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCCCAACGCTTCAAACCCCCCCCcgc                                                                           >  1:3175/1‑149 (MQ=37)
                                                       ccATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc                                                                          <  1:112397/149‑1 (MQ=255)
                                                       ccATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCCCCCCCCCGGCCGGTGCTTCGCAACGCTCAAAACCGCTCCccgcg                                                                         >  1:968804/1‑150 (MQ=255)
                                                        cATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCCGGCGCTGCTTCCCCACGGTTCAATCCCCTccccccgg                                                                        >  2:264962/1‑145 (MQ=255)
                                                        acaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                                                        >  2:400182‑M1/69‑150 (MQ=255)
                                                        acaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                                                        >  2:76578‑M1/69‑150 (MQ=255)
                                                           ccAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGCCCCCTCCCCCCCCCCCCCGGCCGTTGTTTCGCAACGGTCAAAACCCCCCCCCCCGGact                                                                     >  1:837883/1‑148 (MQ=255)
                                                            cagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                    >  1:950180‑M1/65‑150 (MQ=255)
                                                            cagagaaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGAGGACTGCCCGCCACCCTCAGGGAAGTTGATTAGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                    >  1:3966‑M1/65‑150 (MQ=255)
                                                            cAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                     <  1:776989/149‑1 (MQ=255)
                                                             aGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTg                                                                   <  1:835517/150‑1 (MQ=255)
                                                                 aaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGCTGCTTCCCAACATTTAAATCCCCTCCCCGCGGGTTTGTCCt                                                               <  1:614721‑M1/91‑1 (MQ=255)
                                                                  aaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCt                                                               >  1:458259‑M1/59‑149 (MQ=255)
                                                                  aaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                              >  1:274018‑M1/59‑150 (MQ=255)
                                                                   aacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCCCCCGGGCGGTGCTTCGCAAAGTTCAAAACCGCCCCcggcgggttttttcttc                                                             <  1:644563‑M1/93‑8 (MQ=255)
                                                                   aacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                             >  2:482312‑M1/58‑150 (MQ=255)
                                                                    acgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                                            >  2:687100‑M1/57‑150 (MQ=255)
                                                                     caTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTc                                                           <  2:996790/150‑1 (MQ=255)
                                                                      gaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                          >  1:452142‑M1/55‑150 (MQ=255)
                                                                       tCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCAACTCAg                                                         <  2:978393/150‑1 (MQ=255)
                                                                        aaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAgg                                                        >  2:499014‑M1/53‑150 (MQ=255)
                                                                           gTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGCGGGCCTCCCCCCCCCCCTTCCGGCCCTTGGTTTTCAACCACTAAATACGCACCCGGCGGGTTAGTGCTCCTCTggg                                                       >  1:23570/1‑147 (MQ=255)
                                                                            tGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                    <  1:936530/150‑1 (MQ=255)
                                                                            tGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                    <  2:158313/150‑1 (MQ=255)
                                                                            cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                    >  1:246373‑M1/49‑150 (MQ=255)
                                                                               agtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAACCTTCAAATCCGCTCCCCGCGGGTTTTTTCTACTCAGGAGAGcgg                                                 <  1:318099‑M1/105‑2 (MQ=255)
                                                                               aGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCCGCCCGCCCCCCCCCCCGGCGTTGGTTTGCGACGGTTCAAAACCCTCCCCGCGGGTTTTTTCTTCCCCGGAGAGCgg                                                 >  2:658904/1‑149 (MQ=255)
                                                                                 tctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                               >  2:487660‑M1/44‑150 (MQ=255)
                                                                                 tctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGCCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                               >  2:126009‑M1/44‑150 (MQ=255)
                                                                                  ggCTCCGATAATAAGCTGTCAAACATGAGATTTCCGAAGCCCAATACCGGGCAGTTCATGCCGGGCCTCCTCCCCGCCGCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACACAGGAGAGCGTTCa                                              <  1:324663/149‑1 (MQ=17)
                                                                                   gCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                             <  1:863818/150‑1 (MQ=255)
                                                                                   gCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                             <  1:505395/150‑1 (MQ=255)
                                                                                    cTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                            <  1:385816/150‑1 (MQ=255)
                                                                                      cgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCg                                           >  1:62996‑M1/39‑149 (MQ=255)
                                                                                       cGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGa                                          <  2:968760/149‑1 (MQ=255)
                                                                                        actgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGaca                                        >  1:344418‑M1/37‑150 (MQ=255)
                                                                                             tAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaacaa                                   <  1:377948/150‑1 (MQ=255)
                                                                                                    tCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTTGCAAAGTTCAAATCCCCCCCCCGCGGGGTTGTCCTACTCAGGAGAGCGGTCACCCccaaaaaaaaaaaa                             >  1:98792/1‑143 (MQ=17)
                                                                                                    gttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                            >  2:610020‑M1/25‑150 (MQ=255)
                                                                                                        tatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGt                                                                                                                                   >  1:976145‑M1/21‑46 (MQ=255)
                                                                                                        tatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGt                                                                                                                                   <  2:976145‑M1/26‑1 (MQ=255)
                                                                                                        tatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTg                                                                                                                               >  2:345181‑M1/21‑50 (MQ=255)
                                                                                                        tatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTg                                                                                                                               <  1:345181‑M1/30‑1 (MQ=255)
                                                                                                        aCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                        <  1:924767/150‑1 (MQ=255)
                                                                                                             gatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAg                    >  2:968804‑M1/16‑149 (MQ=255)
                                                                                                              atgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                   >  1:1010426‑M1/15‑149 (MQ=255)
                                                                                                               aaTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                 <  1:154389/150‑1 (MQ=255)
                                                                                                                 cctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                >  2:325786‑M1/12‑149 (MQ=255)
                                                                                                                    cGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGGTGCTTCGCAACGTTTAAATCCCCTCCCCGCGGGTTTGTCCTTCTCAGGGGAGCGGTCACCCACAAAAAACAAAAAAAACGAAAAGCCCAggc            >  1:400182/1‑148 (MQ=14)
                                                                                                                     gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt           >  2:186426‑M1/8‑150 (MQ=255)
                                                                                                                     gaaaattAGTTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt           >  2:114582‑M1/8‑150 (MQ=255)
                                                                                                                     gATCCCCAATTCCTGGCAGTTTATGGCGGGGGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt           <  2:837883/150‑1 (MQ=38)
                                                                                                                     gATCCCCAATTCCTGGCAGTTTATGGCGGGGGGCCTGCCCGCCCCCCCCCCGGCCGTTGCTTCGCCACGGTCAAAAACGCCCCCCGCGGGTTTGTTCCACTCCGGAGAGCGGTCACCCACAAAAAACAGAAAAAAAGAAAAAGCCAGTCt           >  2:368620/1‑150 (MQ=255)
                                                                                                                      aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCGGTGCTTCGCAAAGGTCAAATCCCCCCCCCGCGGGTTTGTCCTACTCAGGAGAGCGGTCACCGCCAAACAACAAATAAAAAGAAAAGccccatcct          >  1:76578/1‑144 (MQ=11)
                                                                                                                       gaattAATTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCttt         >  2:596963‑M1/6‑150 (MQ=255)
                                                                                                                        ccccAATTCCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGCCCTTGCTTTGCAACGCTCCAATCCCCCCCCCGCGGGTTTGTTCTTCCCCGGAGGGGGGTCACCGACAAACAAAAGATAAAAAAAAAAGCCCACTCTTTc        >  2:246373/1‑150 (MQ=2)
                                                                                                                        aattaaTTCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCCGTGCTTCTCCAAGTTCAAAATCCCCCCCCGCGGGTTTTTTCTACTCCGGGGAGCGGTCACCGCCAAACAACAAAAAAAAAAAAAAGACCAGACTTTc        >  2:164473/5‑150 (MQ=2)
                                                                                                                        aattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc        >  1:1011971‑M1/5‑150 (MQ=255)
                                                                                                                         cccAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCGGTGCTTCGCAACGTTCAAAACCCCTCCCCGCGGGGTTTTTCTACTCAGGGGAGCGTTCACCCCCAAAAAACAAATAAAAAAAAAAGCCCAGGCTTTc        >  1:308148/1‑149 (MQ=2)
                                                                                                                           taaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGGTTCTTCGCAACGTTCAAAACCCCCCCCCGCGGGTTTTTCCTTCTCCAGAGAGCGTTCCCCCACAAACAACAGAAAAAAAGAAAAGCCCACGCTTTCCAc     >  1:968443/2‑150 (MQ=2)
                                                                                                                            aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc     <  1:423787/149‑1 (MQ=32)
                                                                                                                              ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa  <  1:5488/150‑1 (MQ=31)
                                                                                                                                                     |                                                                                                                                 
CTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGA  >  CP009273/359980‑360255

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