Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,422,671 A→T 5.1% N306Y (AAC→TAC)  moeA → molybdenum::molybdopterin ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,422,6710AT 5.1% 254.7 / 4.0 78N306Y (AAC→TAC) moeAmolybdenum::molybdopterin ligase
Reads supporting (aligned to +/- strand):  ref base A (37/37);  new base T (2/2);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

GAGGCGTCTCGGTGGGCGAGGCCGACTACCTCGGTGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGACATGCCGGCCGGGTTTGACTGGCCCAAGCCGGGTACGCGCCGCGAATACCTGC  >  NC_002947/2422525‑2422818
                                                                                                                                                  |                                                                                                                                                   
gAGGCGTCTCGGTGGGCGAGGCCGACTACCTCGGTGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAAcc                                                                                                                                                  >  1:719699/1‑150 (MQ=255)
  ggCGTCTCGGTGGGCGAAGACGACTACCTAGGTGCAGCCCTGTGCGAGGCCGGCGAACTGGAGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCg                                                                                                                                                <  2:1007897/150‑1 (MQ=255)
       ctcGGTGGGCGAGGCCGACTACCTCGGTGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTc                                                                                                                                           <  2:280756/150‑1 (MQ=255)
            tGGGCGAGGCCGACTACCTCGGTGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCt                                                                                                                                      <  2:719699/150‑1 (MQ=255)
                 gAGGCCGACTACCTCGGTGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTc                                                                                                                                 >  1:486778/1‑150 (MQ=255)
                   ggCCGACTACCTCGGTGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCAc                                                                                                                               >  1:416796/1‑150 (MQ=255)
                   ggCCGACTACCTCGGTGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCGGACCTTCGGCCACTTCCAGGGCGTGCCGGTGGTTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCAc                                                                                                                               >  2:293160/1‑150 (MQ=255)
                       gACTACCTCGGTGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTAGGCCAATTAAAGGGCGTGGCTTGGATCGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTc                                                                                                                           >  1:1028816/1‑150 (MQ=255)
                        aCTACCTCGGTGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCg                                                                                                                          <  1:1335343/150‑1 (MQ=255)
                        aCTACCTCGGTGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCg                                                                                                                          <  1:186971/150‑1 (MQ=255)
                          tACCTCGGTGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGt                                                                                                                        <  2:245135/150‑1 (MQ=255)
                           aCCTCGGTGCCTCCCTCCGCGAGGCCGCCCAACTCGCGTTGTGGAAGCTTGCGATCAAACCGGGCAAGCCGCTGACCTTCTTCCACTTCCACGGCGTGCCGGTGATTGGCTTGCCGGCTAACCCGGCCTCGACCTTCGTCATCTTCGTTc                                                                                                                       <  2:88141/150‑1 (MQ=255)
                            ccTCGGTGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTct                                                                                                                      >  1:815613/1‑150 (MQ=255)
                            ccTCGGTGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTct                                                                                                                      >  2:1318332/1‑150 (MQ=255)
                            ccTCGGTGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTct                                                                                                                      >  2:1318014/1‑150 (MQ=255)
                             cTCGGTGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTaa                                                                                                                                                    >  1:390652/1‑119 (MQ=255)
                             cTCGGTGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTaa                                                                                                                                                    <  2:390652/119‑1 (MQ=255)
                             cTCGGTGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCa                                                                                                                                <  2:322224/139‑1 (MQ=255)
                             cTCGGTGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCa                                                                                                                                >  1:322224/1‑139 (MQ=255)
                             cTCGGTGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTctg                                                                                                                     >  2:34813/1‑150 (MQ=255)
                                  tGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGAc                                                                                                                <  1:1358951/150‑1 (MQ=255)
                                     cgccCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAg                                                                                                             <  1:1301311/150‑1 (MQ=255)
                                      gccCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCTGTCTGCTGACCAgg                                                                                                            <  2:486778/150‑1 (MQ=255)
                                      gccCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAgg                                                                                                            <  1:1317342/150‑1 (MQ=255)
                                       cccTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGc                                                                                                           >  1:432169/1‑150 (MQ=255)
                                       cccTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGc                                                                                                           <  1:1318014/150‑1 (MQ=255)
                                       cccTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGc                                                                                                           <  1:1318332/150‑1 (MQ=255)
                                             gcgAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAAc                                                                                                                                                   >  2:843553/1‑104 (MQ=255)
                                             gcgAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAAc                                                                                                                                                   <  1:843553/104‑1 (MQ=255)
                                                      gCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGCCCTTGGGCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGcgtc                                                                                            <  2:296047/150‑1 (MQ=255)
                                                        gAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGcgtcgt                                                                                          >  2:1199642/1‑150 (MQ=255)
                                                              gCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCg                                                                                                                          >  1:1322139/1‑112 (MQ=255)
                                                              gCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCg                                                                                                                          <  2:1322139/112‑1 (MQ=255)
                                                                gTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTg                                                                                                                                    >  1:462694/1‑100 (MQ=255)
                                                                gTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTg                                                                                                                                    <  2:462694/100‑1 (MQ=255)
                                                                 ttGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCtgctgc                                                                                               >  1:346737/1‑136 (MQ=255)
                                                                 ttGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCtgctgc                                                                                               <  2:346737/136‑1 (MQ=255)
                                                                   gtgGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGTTGGTGAACCCGGCCTCGACCTTGGTCACATTCGTTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGAt                                                                               >  1:673745/1‑150 (MQ=255)
                                                                    tgGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCTACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTTGACCTTGTTCACGTTCGGTCTGCTAACCAGGCCCGACCTGCTGCGTCGTCAGGGTGTGAtt                                                                              <  2:426197/150‑1 (MQ=255)
                                                                    tgGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGAtt                                                                              <  1:305266/150‑1 (MQ=255)
                                                                     ggAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTctgct                                                                                                                   >  1:121206/1‑112 (MQ=255)
                                                                     ggAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTctgct                                                                                                                   <  2:121206/112‑1 (MQ=255)
                                                                       aaGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACctg                                                                                                   <  2:929424/126‑1 (MQ=255)
                                                                       aaGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACctg                                                                                                   >  1:929424/1‑126 (MQ=255)
                                                                           tGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCa                                                                       >  1:110121/1‑150 (MQ=255)
                                                                             gCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACg                                                                     >  2:1353980/1‑150 (MQ=255)
                                                                               gATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGcc                                                                   >  2:1333392/1‑150 (MQ=255)
                                                                                    aaCCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCa                                                                                        >  1:302314/1‑124 (MQ=255)
                                                                                    aaCCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCa                                                                                        <  2:302314/124‑1 (MQ=255)
                                                                                     aCCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGgt                                                                                    >  1:1229095/1‑127 (MQ=255)
                                                                                     aCCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGgt                                                                                    <  2:1229095/127‑1 (MQ=255)
                                                                                            aaGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGAc                                                      >  2:339709/1‑150 (MQ=255)
                                                                                              gCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGAcc                                                                                                                                       >  2:861968/1‑67 (MQ=255)
                                                                                              gCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGAcc                                                                                                                                       <  1:861968/67‑1 (MQ=255)
                                                                                              gCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAg                                                                                                             <  1:302371/93‑1 (MQ=255)
                                                                                              gCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAg                                                                                                             >  2:302371/1‑93 (MQ=255)
                                                                                              gCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGACAt                                                    <  2:1033917/150‑1 (MQ=255)
                                                                                                        ttCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGAc                                                                           >  1:610475/1‑117 (MQ=255)
                                                                                                        ttCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGAc                                                                           <  2:610475/117‑1 (MQ=255)
                                                                                                          cGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGACATGCCGGCCGGGtt                                        <  1:339709/150‑1 (MQ=255)
                                                                                                            gCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGACATGCCGGCCGGGTTTg                                      >  2:121542/1‑150 (MQ=255)
                                                                                                             ccACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGACATGCCGGCCGGGTTTGa                                     <  1:34813/150‑1 (MQ=255)
                                                                                                             ccACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCGGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGACATGCCGGCCGGGTTTGa                                     >  2:1052870/1‑150 (MQ=255)
                                                                                                                 ttCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGACATGCCGGCCGGGTTTGACTgg                                 <  1:1353980/150‑1 (MQ=255)
                                                                                                                    cAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTAGGTCTGCTGACCAGCTCTTTCCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGACATGCCGGCCGGGTTTGACTGGccc                              >  1:884127/1‑150 (MQ=255)
                                                                                                                              cGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGACATGCCGGCCGGGTTTGACTGGCCCAAGCCGGGTa                    >  1:585780/1‑150 (MQ=255)
                                                                                                                                  gATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGACATGCCGGCCGGGTTTGACTGGCCCAAGCCgggaccccg                >  2:678731/1‑144 (MQ=255)
                                                                                                                                     tGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTAAATGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGACATGCCGGCCGGGTTTGACTGGCCCAAGCCGGGTACGcgccg             >  1:1098569/1‑150 (MQ=255)
                                                                                                                                       gCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGACATGCCGGCCGGGTTTGACTGGCCCAAGCCGGGTACGCGCcgcg           >  1:690464/1‑150 (MQ=255)
                                                                                                                                         aggCCGGGTTACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGcc                                                                                                          <  2:39735/51‑1 (MQ=255)
                                                                                                                                         aggCCGGGTTACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGcc                                                                                                          >  1:39735/3‑53 (MQ=255)
                                                                                                                                          ggCCGGGTTACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCTGTTCGACATgccggc                                              <  2:1417163/111‑1 (MQ=255)
                                                                                                                                          ggCCGGGTTACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGACATgccggc                                              >  1:1417163/2‑112 (MQ=255)
                                                                                                                                            ccGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGACATGCCGGCCGGGTTTGACTGGCCCAAGCCGGGTACGCGCCGCGAATAc      <  2:673745/150‑1 (MQ=255)
                                                                                                                                             cGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGCCATGCCGGCCGGGTTTGACTGGCCCaa                            >  1:446732/1‑127 (MQ=255)
                                                                                                                                             cGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGACATGCCGGCCGGGTTTGACTGGCCCaa                            <  2:446732/127‑1 (MQ=255)
                                                                                                                                              gggTAACCCGGCCTCGACCTTAGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGACATGCCGGCCGGGTTTGACTGGCCCAAGCCGGGTACGCGCCGCGaa         <  2:452377/145‑1 (MQ=255)
                                                                                                                                              gggTAACCCGGCCTCGACCTTAGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGACATGCCGGCCGGGTTTGACTGGCCCAAGCCGGGTACGCGCCGCGaa         >  1:452377/1‑145 (MQ=255)
                                                                                                                                                gTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGACATGCCGGCCGGGTTTGACTGGCCCAAGCCGGGTACGCGCCGCGAATACCTgc  <  2:432169/150‑1 (MQ=255)
                                                                                                                                                gTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGACATGCCGGCCGGGTTTGACTGGCCCAAGCCGGGTACGCGCCGCGAATACCTgc  <  2:1028816/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GAGGCGTCTCGGTGGGCGAGGCCGACTACCTCGGTGCCGCCCTGCGCGAGGCCGGCGAACTGGCGTTGTGGAAGCTGGCGATCAAACCGGGCAAGCCGCTGACCTTCGGCCACTTCCAGGGCGTGCCGGTGATTGGCTTGCCGGGTAACCCGGCCTCGACCTTGGTCACCTTCGGTCTGCTGACCAGGCCCTACCTGCTGCGTCGTCAGGGTGTGATTGACGTCACGCCGTTGCGGTTCGACATGCCGGCCGGGTTTGACTGGCCCAAGCCGGGTACGCGCCGCGAATACCTGC  >  NC_002947/2422525‑2422818

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