Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_012971 788,122 A→G intergenic (+62/‑3) lom → / → ECD_RS23105 outer membrane beta‑barrel protein Lom/prophage tail fiber N‑terminal domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_012971788,1220AG91.7% 115.0 / ‑3.2 36intergenic (+62/‑3)lom/ECD_RS23105outer membrane beta‑barrel protein Lom/prophage tail fiber N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (22/11);  total (23/13)
Fisher's exact test for biased strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTTGTTGATATTGCTTATGAAGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCG  >  NC_012971/787975‑788420
                                                                                                                                                   |                                                                                                                                                                                                                                                                                                          
tggccattgaccttgcttatgaaggtTCCGGCAGTGGTGACTGGCGCACTGACGGTTTCATCGTGGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttac                                                                                                                                                    >  2:22090‑M2/13‑156 (MQ=255)
  gccatttaccttgcgtatgaagGTTCCGGCAGTGGTGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaanctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacag                                                                                                                                                  >  2:575062‑M2/11‑154 (MQ=255)
     attgaccttgcttatgAAGGTTCCGGCAGTGGTGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcat                                                                                                                                               >  1:79686‑M2/8‑151 (MQ=255)
                ttATGAAGGTTCCGGCAGTGGTGACTGGCGCACGGACGGTTTCATCGTGGGTGTCGGTTATAAATTCTGATTAGCCANGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcactacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgag                                                                                                                                     >  1:154344‑M2/1‑140 (MQ=255)
                         ttcCGGCAGTGGTGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccaaggtactggtgaacacggtagcatatgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggcaagt                                                                                                                           <  2:755412‑M2/299‑170 (MQ=255)
                             ggCAGTGGTGACTTGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCNCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtac                                                                                                                         >  2:627309‑M2/1‑127 (MQ=255)
                             ggCAGTGGTGACTTGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtac                                                                                                                         >  2:805674‑M2/1‑127 (MQ=255)
                             ggCAGTGGTGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtac                                                                                                                         >  2:805424‑M2/1‑127 (MQ=255)
                             ggCAGTGGTGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtac                                                                                                                         >  1:72106‑M2/1‑127 (MQ=255)
                             ggCAGTGGTGACTGGCGCACTGACGGGTTCAGCGTGGGTGTCGGTGATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCNCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtac                                                                                                                         >  2:627299‑M2/1‑127 (MQ=255)
                                aGTGGTGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgt                                                                                                                    >  2:702314‑M2/1‑124 (MQ=255)
                                 gtggtgACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagcagggcgttacagcatggatgtcgagtatggcaagtacagcgtc                                                                                                                   <  2:95022‑M2/300‑178 (MQ=255)
                                     tgaCTGGCGCACTGACTGTTTCATCGTGGGTGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccantacagaactgcactattcaactgaaggcaaagngtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccc                                                                                                               >  2:616661‑M2/2‑119 (MQ=255)
                                              cacTGACGGTTTCATCGTGGGTGTCGGTTATAAATTCTGATTAGCAAGGTAACACAGTGTTATGACGGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaataccgaactgcactcttcaactgaaggcaaagagtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaacatgctggtt                                                                                                       <  2:175014‑M2/298‑190 (MQ=255)
                                                  gACGGTTTCATCGTGGGTGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgt                                                                                                                    >  1:542656‑M2/1‑106 (MQ=255)
                                                        ttCATCGTGGGTGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaantccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttcc                                                                                            <  1:170990‑M2/300‑201 (MQ=255)
                                                                      cGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaacctgctggttgaaggttttccgccttcacatgca                                                                               <  2:646520‑M2/299‑214 (MQ=255)
                                                                                ttCTGATGATCCGGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgcttaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacaacgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattac                                                                    >  2:375969‑M2/1‑76 (MQ=255)
                                                                                   tGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccg                                                                  >  1:242031‑M2/1‑73 (MQ=255)
                                                                                    gATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGATTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtaaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtc                                                                >  1:665674‑M2/1‑72 (MQ=255)
                                                                                         gCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatga                                                           <  1:695551‑M2/300‑234 (MQ=255)
                                                                                            aGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaa                                                          >  1:751072‑M2/1‑64 (MQ=255)
                                                                                            aGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgcttgttgaaggttttccgccttcacatgccggaaccattaccgtctatgaa                                                          <  1:319081‑M2/298‑235 (MQ=255)
                                                                                            aGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaagg                                                        >  1:549778‑M2/1‑64 (MQ=255)
                                                                                            aGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaa                                                          >  1:649036‑M2/1‑64 (MQ=255)
                                                                                            aGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaa                                                          >  2:44941‑M2/1‑64 (MQ=255)
                                                                                              gTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctaaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgacggcaaagcgtaacagcaccagggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaacgtatatgaaggtt                                                      <  2:206416‑M2/300‑239 (MQ=255)
                                                                                                     aGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaaggttccagacc                                               >  1:89806‑M2/1‑55 (MQ=255)
                                                                                                     aGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaaggttccagac                                                >  1:101413‑M2/1‑55 (MQ=255)
                                                                                                      gtgtTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaacgtctatgaaggttccagacct                                              <  1:99502‑M2/300‑247 (MQ=255)
                                                                                                      gtgtTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagcagggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgacttcacatgccggaaccattaccgtctatgaaggttccagacct                                              <  2:10682‑M2/300‑247 (MQ=255)
                                                                                                         ttATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGAGTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtggggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaaggttccagacctggta                                           >  1:201264‑M2/1‑51 (MQ=255)
                                                                                                                     gaagatatcgtaccggggccgtcagctccgaatatggcagcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaaggttccagacctggtacgctgaatg                                 >  1:535885‑M2/30‑39 (MQ=255)
                                                                                                                        gatatcgtaccggggccgtcagctccgaatatggcagcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaaggttccagacctggtacgctgaatgatttt                            <  1:617557‑M2/274‑265 (MQ=255)
                                                                                                                         atatcgtaccggggccgtcagctccgaatatggcagcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaaggttccagacctggtacgctgaatgattttc                           <  1:11122‑M2/275‑266 (MQ=255)
                                                                                                                          tatcgtaccggggccgtcagctccgaatatggcagcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaac                                                                          >  1:72330‑M2/25‑34 (MQ=255)
                                                                                                                          tatcgtaccggggccgtcagctccgaatatggcagcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcatctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaac                                                                          <  2:72330‑M2/228‑219 (MQ=255)
                                                                                                                            tcgtaccggggccgtcagctccgaatatggcagcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcgctattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaaggttccagacctggtacgctgaatgattttctcg                        >  1:417302‑M2/23‑32 (MQ=255)
                                                                                                                            tcgtaccggggccgtcagctccgaatatggcagcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaaggttccagacctggtacgctgaatgattttctcg                        >  1:417190‑M2/23‑32 (MQ=255)
                                                                                                                              gtaccggggccgtcagctccgaatatggcagcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaaggttccagacctggtacgccgaatgattttctcggt                      <  1:227137‑M2/280‑271 (MQ=255)
                                                                                                                                acctgtgccgtcagctccgaatatggcagcattaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcttaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatttcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgcctccacatgccggaaccattaccgtctatgaaggttccagacctggtacgctgaatgattttctcggtgc                    >  2:167036‑M2/19‑28 (MQ=255)
                                                                                                                                   ggggccgtcagctccgaatatggcagcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaaggttccagacctggtacgctgaatgattttctcggtgccat                 <  2:73696‑M2/285‑276 (MQ=255)
                                                                                                                                     ggccgtcagctccgaatatggcagcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaaggttccagacctggtacgctgaatgattttctcggtgccatga               >  1:483654‑M2/14‑23 (MQ=255)
                                                                                                                                             gctccgaatatGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatcggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaaggttccagacctggtacgctgaatgattttctcggtgccatgaaggaagat        <  2:783454‑M2/294‑285 (MQ=255)
                                                                                                                                             gctccgaatatGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaaggttccagacctggtacgctgaatgattttctcggtgccatgacggaagat       >  1:793356‑M2/6‑15 (MQ=255)
                                                                                                                                                gTGGATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctctgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaaggttccagacctggtacgctgaatgattttctaggtgccatgacggaagatgat    <  2:796282‑M2/300‑289 (MQ=255)
                                                                                                                                                  gAATATGGCAgcagtaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctatgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaaggttccagacctggtacgctgaatgattttctcggtgccatgacggaagatgatgt  <  2:189752‑M2/300‑291 (MQ=255)
                                                                                                                                                  gAATATGGCAgcagtaaaaatctaaggtgtgctgaaagatggtgcgggaaaaccaatacagaactgcactattcaactgaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctatgaaaatccggatgaagccgggcgttacagcatggatgtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaaggttccagacctggtacgctgaatgattttctcggtgccatgaaggaagatgatgt  <  2:186811‑M2/300‑291 (MQ=255)
                                                                                                                                                  gAATAGGGCAgcattaaaaatctcaggtgtgctgaaagatggtgcgggaaaaccaatacgggattgtactatttaatggaaggcaaagcgtaacagcaccacggtactggtgaacacggtagcctccgaaaatccggatgaagccgggcgttacagcatggaggtcgagtatggccagtacagcgtcaccctgctggttgaaggttttccgccttcacatgccggaaccattaccgtctatgaaggttccagacctggtacgctgaatgattttctcggtgccatgacggaagatgatgt  >  2:699105‑M2/1‑10 (MQ=255)
                                                                                                                                                   |                                                                                                                                                                                                                                                                                                          
TCGTTGTTGATATTGCTTATGAAGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCG  >  NC_012971/787975‑788420

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