Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,181 G→C 11.9% 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,1810GC11.9% 313.0 / 4.6 105intergenic (‑115/+480)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (43/48);  new base C (12/2);  total (55/50)
Fisher's exact test for biased strand distribution p-value = 8.88e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

TACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTAT  >  CP009273/360033‑360330
                                                                                                                                                     |                                                                                                                                                     
tattccgcaaacaactgagaaatttttaattcgctctgaaactgatggcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggcg                                                                                                                                                       <  1:697412‑M1/2‑1 (MQ=255)
tattccgcaaacaactgagaaatttttaattcgctctgaaactgatggcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggcg                                                                                                                                                       <  2:1059298‑M1/2‑1 (MQ=255)
  ttccgcaaacaactgagaaatttttaattcgctctgaaactgatggcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggcggc                                                                                                                                                     <  2:671397‑M1/4‑1 (MQ=255)
     aacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                                                 >  1:918860‑M1/67‑150 (MQ=255)
      gcaaacaactgagaaatttttaattcgctctgaaactgatggcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGAt                                                                                                                                                 <  2:1019171‑M1/8‑1 (MQ=255)
       cagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                                                                                               >  1:485566‑M1/65‑150 (MQ=255)
          ataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                                            >  1:557201‑M1/62‑150 (MQ=255)
           gcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCCCCCCCCCGGCCGGTGCTTCGCAAAGTTCAAAAACCCCCCCCGCGGGttttttc                                                                                                                                           >  1:110743/1‑146 (MQ=25)
           gcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCCCCCGGCCGGTGCTTCCCAACGTTCAAATCCCCCCCCCGCGGGTTTGTcc                                                                                                                                           >  1:8206/1‑150 (MQ=25)
             caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                                         <  1:973151/150‑1 (MQ=255)
             caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                                         <  1:321110/150‑1 (MQ=255)
             caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCCGCCCGTGCTTCCCCACGCTCCAAACCCCCCCCCGCCGGTTTTTTctt                                                                                                                                          >  2:1021208/1‑148 (MQ=255)
               ccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGCCGCTTCTTCGCAAAGTTCCAAACCCCCCCCGGCGGGTTTGTCCTAcc                                                                                                                                       >  1:679494/1‑149 (MQ=25)
                 gaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                                                     >  1:1021208‑M1/55‑150 (MQ=255)
                 gaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCACCCCCCCGGCCGGTGCTTCGCCACGTTCAAAAACCCTCCCCGCGGGGTTGTGCTACTc                                                                                                                                      <  2:725094‑M1/95‑1 (MQ=255)
                 agaaatttttaattcgctctgaaactgatggcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACt                                                                                                                                       <  1:295817‑M1/18‑1 (MQ=255)
                 agaaatttttaattcgctctgaaactgatggcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACt                                                                                                                                       <  1:927183‑M1/18‑1 (MQ=255)
                 agaaatttttaattcgctctgaaactgatggcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACt                                                                                                                                       >  2:101035‑M1/132‑149 (MQ=255)
                 aTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCCCCCGGCCGTTGCTTCGCAACGTTCAAATCCGCCCCCCGCGGGTTTGTCCTACt                                                                                                                                       >  1:772160/1‑148 (MQ=25)
                  tCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAg                                                                                                                                    <  2:579129/150‑1 (MQ=255)
                  aaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTTCCCGCCACCCCCCCGGCCGTTGCTTTCCAACGGTCAAATCCGCTCCCCGCGGGTTTTTTCTActccc                                                                                                                                    <  2:1034543‑M1/97‑3 (MQ=255)
                  aaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAg                                                                                                                                    >  1:868600‑M1/54‑150 (MQ=255)
                    aatttttaattcgctctgaaactgatggcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAg                                                                                                                                    <  1:977343‑M1/21‑1 (MQ=255)
                        ccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                                              >  2:763561‑M1/48‑150 (MQ=255)
                         cAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                                              <  2:132856/149‑1 (MQ=255)
                         cAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                                              <  1:713371/149‑1 (MQ=255)
                          taattcgctctgaaactgatggcgtaaaaggaaagaaagaagaacctttataccatgtcaagaagtaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCg                                                                                                                             <  1:1028386‑M1/28‑1 (MQ=255)
                          taattcgctatgaaaatgatggcgtagaagagaagaaataataccatatataccatgtcccgtattaccgcctttgagtgagatgataccgctagacgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                                              >  1:911900‑M1/123‑149 (MQ=255)
                             aGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAAAGTTCAAATCCGCTCCCCGCGGATTTTTCCTACTCAGGAGAGCGTTCa                                                                                                                         >  1:28180/1‑150 (MQ=21)
                              gCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCCCCCGGCCGTTGCTTCGCAAAGTTCAAATCCGCTCCCCGCGGGTTTGTCCTTCTCAGGGGAGCGGTCAc                                                                                                                        >  1:508459/1‑150 (MQ=18)
                               cgctctgaaactgatggcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                                                        <  1:1001021‑M1/33‑1 (MQ=255)
                               cTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCCCCCCCCCGGCCCTTGCTTTGCAAAGGTTAAATCCCCCCCCCGCGGGTTTTTCCCACTCCAGGGAGCGGTcccc                                                                                                                       >  2:905258/1‑150 (MQ=16)
                                   gatgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCACCCCCCGCGCCCTTGGTTCCCAACATTCAAATCCGCCCCCCGCGGGGTTGTCCTACCCAGGAGAGCGTTCACCCac                                                                                                                    >  1:802288/1‑149 (MQ=18)
                                   actgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGaca                                                                                                                   >  1:514288‑M1/37‑150 (MQ=255)
                                     tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaa                                                                                                                 <  2:391400/150‑1 (MQ=255)
                                      aaactgatggcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaa                                                                                                                 <  1:965562‑M1/40‑1 (MQ=255)
                                       aactgatggcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaa                                                                                                                 <  1:675990‑M1/40‑1 (MQ=255)
                                          tgatggcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaacaa                                                                                                              <  1:534460‑M1/43‑1 (MQ=255)
                                          ctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                                                            >  1:358203‑M1/30‑150 (MQ=255)
                                          aGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAaacaac                                                                                                             <  1:324052/149‑1 (MQ=255)
                                             tggcgtagaagcgaagaaataataccctttataccatgttccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaacaa                                                                                                              <  2:276626‑M1/43‑1 (MQ=255)
                                             tcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGa                                                                                                          >  2:772160‑M1/27‑149 (MQ=255)
                                             tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGTCCTGCCCCCCCCCCCCCCCGCCGCTGCTTTTCAACGTTCAAATACCCCCCCCGCGGGGTTTTCCTCCTCCCGAGAGCGGTTTCCGcccaacaaacaa                                                                                                          >  2:631983/1‑145 (MQ=255)
                                              gTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATTGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                                                         >  2:183484/1‑149 (MQ=255)
                                               tCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                                                        <  1:667217/149‑1 (MQ=255)
                                               tCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                                                        <  2:273143/149‑1 (MQ=255)
                                               gcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                                                        >  1:1019171‑M1/102‑150 (MQ=255)
                                                 gtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                                                       <  2:562546‑M1/50‑1 (MQ=255)
                                                    gaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                                                   >  2:421897‑M1/97‑150 (MQ=255)
                                                     aagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                                                   <  1:505394‑M1/54‑1 (MQ=255)
                                                     aagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                                                   <  1:1053925‑M1/54‑1 (MQ=255)
                                                      tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaaa                                                                                                <  2:369342/150‑1 (MQ=255)
                                                        gatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAg                                                                                               >  1:725094‑M1/16‑149 (MQ=255)
                                                         atgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                                                                                              >  2:902484‑M1/15‑149 (MQ=255)
                                                           aTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                                                                                           <  1:747962/150‑1 (MQ=255)
                                                             aaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                                                            <  1:305801‑M1/61‑1 (MQ=255)
                                                                gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc                                                                                       >  2:843657‑M1/8‑149 (MQ=255)
                                                                gATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                                      <  1:123600/150‑1 (MQ=38)
                                                                 ggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                                     >  1:693744‑M1/7‑150 (MQ=255)
                                                                 aataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                                      <  2:666670‑M1/67‑1 (MQ=255)
                                                                  gaattaaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGCTGCTTTGCAACGTTCAAATCCCCCCCCCCCGGGTTTTTTCTACCCAGGGGGGCGGTCACCGACAAACAACAGACAAAAAGAAAAGCGCAGTCttt                                                                                    >  2:918860/6‑150 (MQ=2)
                                                                  ataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                                     >  2:1053925‑M1/83‑150 (MQ=255)
                                                                   taccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCttt                                                                                    <  1:1018742‑M1/69‑1 (MQ=255)
                                                                   ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGTTGCTTCCCAACGTTCAAATCCCCTCCCCGCGGGTTTTTTCTACTCAGGAGAGCGTTTACCGACAAAAAAAAAATAAAACGAAAAGCCCAGTCtttt                                                                                   >  1:819042/1‑149 (MQ=11)
                                                                   aattaaTTCCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCTCCCGGCCGTTGCTTCTCAACGTTCAAAACCACTCCCCGCGGGGTTTTTCTACTCAGGAGAGCGGTTACCGACAAACAAAAGAAAAAAAAAAAAGACCAGTCtttt                                                                                   >  2:715020/5‑149 (MQ=2)
                                                                    cccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                                   <  1:1025007/149‑1 (MQ=35)
                                                                     ccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                                                  >  1:442860‑M1/80‑150 (MQ=255)
                                                                       ctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                                                <  1:367027‑M1/73‑1 (MQ=255)
                                                                        aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                                              <  1:838730/150‑1 (MQ=32)
                                                                        aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                                              <  2:28180/150‑1 (MQ=32)
                                                                          tataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                                               <  2:167738‑M1/74‑1 (MQ=255)
                                                                            taccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAg                                                                            >  2:546305‑M1/73‑149 (MQ=255)
                                                                            taccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGc                                                                           <  2:87247‑M1/78‑1 (MQ=255)
                                                                              ccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCt                                                                         <  2:117343‑M1/80‑1 (MQ=255)
                                                                                atgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCttt                                                                       <  1:67505‑M1/82‑1 (MQ=255)
                                                                                atgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCttt                                                                       >  2:867589‑M1/69‑150 (MQ=255)
                                                                                  gtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCg                                                                     <  1:300381‑M1/84‑1 (MQ=255)
                                                                                          ttaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTAtt                                                              >  1:403282‑M1/59‑149 (MQ=255)
                                                                                              cgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTAtttgctg                                                         >  1:833401‑M1/55‑147 (MQ=255)
                                                                                              cgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATg                                                         >  1:833404‑M1/55‑150 (MQ=255)
                                                                                               gcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGc                                                        <  2:428022‑M1/97‑1 (MQ=255)
                                                                                                cctttgagtgagctgatgccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGcc                                                       >  1:673770‑M1/53‑150 (MQ=255)
                                                                                                 ctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCt                                                      <  1:719407‑M1/99‑1 (MQ=255)
                                                                                                     cccGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                                                   <  1:148880/148‑1 (MQ=11)
                                                                                                     cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                                                   <  2:679494/148‑1 (MQ=11)
                                                                                                      ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                                                >  1:855477/1‑150 (MQ=32)
                                                                                                        gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                              >  1:660989/1‑150 (MQ=34)
                                                                                                                  gggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaat                                    >  1:281741/1‑150 (MQ=255)
                                                                                                                     ccGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCt                                 >  1:592401/1‑150 (MQ=255)
                                                                                                                     ccGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCt                                 <  1:80602/150‑1 (MQ=255)
                                                                                                                           cTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGa                           <  1:961/150‑1 (MQ=255)
                                                                                                                            ttCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGAc                          >  1:320078/1‑150 (MQ=255)
                                                                                                                            ttCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGAc                          >  2:275366/1‑150 (MQ=255)
                                                                                                                              cGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGAc                          <  2:887896/148‑1 (MQ=255)
                                                                                                                               gCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGa                        >  2:340614/1‑149 (MQ=255)
                                                                                                                                        cAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTc               >  1:425999/1‑150 (MQ=255)
                                                                                                                                          aaTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGt             >  1:467985/1‑150 (MQ=255)
                                                                                                                                          aaTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGt             >  2:517396/1‑150 (MQ=255)
                                                                                                                                                   ccGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCt    <  2:65247/150‑1 (MQ=255)
                                                                                                                                                     ggcggATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTAt  <  2:855477/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTAT  >  CP009273/360033‑360330

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