Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,840,330 A→G 6.3% A288A (GCA→GCG)  ynjE → molybdopterin synthase sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,840,3300AG 6.3% 196.6 / 3.1 64A288A (GCA→GCG) ynjEmolybdopterin synthase sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (32/27);  new base G (2/2);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CGCAGATTATGCTTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGT  >  NC_000913/1840185‑1840465
                                                                                                                                                 |                                                                                                                                       
cgcAGATTATGCTTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACaa                                                                                                                                      >  4:189739/1‑149 (MQ=255)
       tatgctttatgctGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGt                                                                                                                               >  5:102989/1‑149 (MQ=255)
           ctttatgctGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGAt                                                                                                                           <  1:312137/149‑1 (MQ=255)
              tatgctGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATg                                                                                                                          <  2:343472/147‑1 (MQ=255)
                         gAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGaa                                                                                                              >  4:206828/1‑148 (MQ=255)
                          aaaGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCGCAACCGCAGTTGATGCTCGATATGGaa                                                                                                              >  1:146015/1‑147 (MQ=255)
                          aaaGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCGCAACCGCAGTTGATGCTCGATATGGaa                                                                                                              >  1:140139/1‑147 (MQ=255)
                          aaaGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGata                                                                                                                   >  4:360894/1‑142 (MQ=255)
                          aaaGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGata                                                                                                                   <  3:360894/142‑1 (MQ=255)
                          aaaGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGaa                                                                                                              >  7:88410/1‑147 (MQ=255)
                           aaGATGTGCGCCTGCTGGATGGCGGCGGGCAAACCTGGTCCGACGCGGGACTGCGTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAg                                                                                                           <  7:372745/149‑1 (MQ=255)
                                       tgctgGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGActg                                                                                               >  7:282814/1‑149 (MQ=255)
                                          tgGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGActgctg                                                                                            >  1:60549/1‑149 (MQ=255)
                                             aTGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCAt                                                                                         <  5:207580/149‑1 (MQ=255)
                                             aTGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCAt                                                                                         <  8:105659/149‑1 (MQ=255)
                                                  ggcTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCa                                                                                    <  2:60549/149‑1 (MQ=255)
                                                   gcTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGGGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAg                                                                                   <  3:230131/149‑1 (MQ=255)
                                                     tGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGa                                                                                 >  4:142934/1‑149 (MQ=255)
                                                             cTGGTCCGACACGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGt                                                                                                                               >  5:258202/1‑95 (MQ=255)
                                                             cTGGTCCGACACGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGt                                                                                                                               <  6:258202/95‑1 (MQ=255)
                                                                 tCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGt                                                                     <  6:102989/149‑1 (MQ=255)
                                                                 tCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGt                                                                     >  2:367665/1‑149 (MQ=255)
                                                                            aCTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGAt                                                                                >  8:197543/1‑127 (MQ=255)
                                                                            aCTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGAt                                                                                <  7:197543/127‑1 (MQ=255)
                                                                                ccTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATtcgttcgt                                                       <  8:282814/148‑1 (MQ=255)
                                                                                 cTGTTGAGCGCGGAACGCCGCCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTg                                                      >  7:264986/1‑148 (MQ=255)
                                                                                 cTGTTGAGCGCGGAACGCCGCCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTg                                                      <  8:264986/148‑1 (MQ=255)
                                                                                   gTTGATCGCGGAACGCCACCGAAAGTGAAAGCGGAACCCCATTTCGTCGTGAAGATCCCGGCACAACAGCAGTTAATCCTCGATAAGGAACAGGCCCGTGGACTGATGCCACGCCAGGATGCAACGCTGGAGACCAGTCGTTCGGGGc                                                    >  2:290753/1‑148 (MQ=255)
                                                                                    ttGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGc                                                    >  2:97060/1‑147 (MQ=255)
                                                                                        gcgcGTAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAAt                                              >  5:155365/1‑149 (MQ=255)
                                                                                         cgcgGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAAtt                                             >  8:246450/1‑149 (MQ=255)
                                                                                          gcgGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAAtt                                             >  4:179291/1‑148 (MQ=255)
                                                                                            ggAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTAt                                          >  3:103507/1‑149 (MQ=255)
                                                                                             gAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCGCAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATc                                         <  2:146015/149‑1 (MQ=255)
                                                                                             gAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCGCAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATc                                         <  2:140139/149‑1 (MQ=255)
                                                                                                 gCCACCGAAAGTGAAATCGGAACCGGATTTCGGCGTGTAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCGGTa                                     <  2:20625/149‑1 (MQ=255)
                                                                                                  ccACCGAAAGTGTAAGCGGAACCGGATTTCGGCTTGAAGATCCGGGCACATCCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGTTGCATCGCTGGTGAGCATTGGTTCGTGGCCAGAATTTATCGGTAc                                    <  1:290753/149‑1 (MQ=255)
                                                                                                   cACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCCCAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCGGTAc                                    >  1:39590/1‑148 (MQ=255)
                                                                                                   cACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGATCAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTgg                                                     <  5:210971/131‑1 (MQ=255)
                                                                                                   cACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTgg                                                     >  6:210971/1‑131 (MQ=255)
                                                                                                    aCCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACCGGCCCGTGGCCTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGGCAGAATTTATCGGTACgc                                  >  7:321881/1‑148 (MQ=255)
                                                                                                      cGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAgg                                                                                  <  5:145975/99‑1 (MQ=255)
                                                                                                      cGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAgg                                                                                  >  6:145975/1‑99 (MQ=255)
                                                                                                      cGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCGGTACGAc                                 >  6:284266/1‑148 (MQ=255)
                                                                                                       gAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCg                                        >  7:353087/1‑140 (MQ=255)
                                                                                                       gAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCg                                        <  8:353087/140‑1 (MQ=255)
                                                                                                          aGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCGGTACGACCAGCg                            <  3:206828/149‑1 (MQ=255)
                                                                                                            tgaacacGGAACCGGATTTAGTCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTCGACTGCTGCAACGCCAGGCTGCATCGCTGGTGAGCAATCGTACGTTGCCAGAATTTATCGGTACGACCAGCgg                           >  4:39810/7‑148 (MQ=255)
                                                                                                              aaaGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTCATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCg                                        <  7:267208/133‑1 (MQ=255)
                                                                                                              aaaGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTCATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCg                                        >  8:267208/1‑133 (MQ=255)
                                                                                                                    gAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCGGTACGACCAGCg                            >  7:214387/1‑139 (MQ=255)
                                                                                                                    gAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCGGTACGACCAGCg                            <  8:214387/139‑1 (MQ=255)
                                                                                                                    gAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCGGTACGACCAGCGGTTACAGcga                  >  1:335576/1‑147 (MQ=255)
                                                                                                                     aaCCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCGGTACGACCAGCGGTTACAGCta                  >  3:222191/1‑148 (MQ=255)
                                                                                                                        cGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCGGTACGACCAGCGGTTACAGCTATAtt              <  4:103507/149‑1 (MQ=255)
                                                                                                                        cGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCGGTACGACCAGCGGTTACAGCTATAtt              <  1:367665/149‑1 (MQ=255)
                                                                                                                         ggATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCGGTACGACCAGCGGTTACAGCTATAtt              >  1:276892/1‑148 (MQ=255)
                                                                                                                          gATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCGGTACGACCAGCGGTTACAGCTATATTaa            <  4:222191/149‑1 (MQ=255)
                                                                                                                                 gCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCGGTACGACCAGCGGTTACAGCTATATTAAACCaaaa     <  3:189739/149‑1 (MQ=255)
                                                                                                                                   gTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTAt                                          <  5:17785/110‑1 (MQ=255)
                                                                                                                                   gTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTAt                                          >  6:17785/1‑110 (MQ=255)
                                                                                                                                    tGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGt  >  2:285813/1‑149 (MQ=255)
                                                                                                                                    tGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGt  <  2:39590/149‑1 (MQ=255)
                                                                                                                                     gAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGt  >  3:149676/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                       
CGCAGATTATGCTTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGTTGATGCTCGATATGGAACAGGCCCGTGGACTGCTGCATCGCCAGGATGCATCGCTGGTGAGCATTCGTTCGTGGCCAGAATTTATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGT  >  NC_000913/1840185‑1840465

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