Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 548,374 G→A 47.1% A6A (GCG→GCA)  ylbE → NO‑induced DUF1116 protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913548,3740GA47.1% 16.6 / 111.0 85A6A (GCG→GCA) ylbENO‑induced DUF1116 protein
Reads supporting (aligned to +/- strand):  ref base G (21/24);  new base A (22/18);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

TGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAG  >  NC_000913/548277‑548470
                                                                                                 |                                                                                                
tgCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCa                                                                                                >  2:900355/1‑100 (MQ=255)
 gCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACaa                                                                                               >  2:1610413/1‑100 (MQ=255)
  cATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAg                                                                                              >  2:312183/1‑100 (MQ=255)
  cATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAg                                                                                              <  1:449415/100‑1 (MQ=255)
   aTTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGAGATTCCCATGTTTACATCAGTGGCACAAGc                                                                                             >  2:392221/1‑100 (MQ=255)
    ttACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGcc                                                                                            <  1:1361920/100‑1 (MQ=255)
    ttACCAATGGTCGCCAGTCGACGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGcc                                                                                            <  1:2038574/100‑1 (MQ=255)
     tACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGcc                                                                                            >  1:1148088/1‑99 (MQ=255)
      aCCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCaa                                                                                          >  1:409912/1‑100 (MQ=255)
      aCCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCaa                                                                                          >  2:2196122/1‑100 (MQ=255)
       ccAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAAt                                                                                         <  2:1126122/100‑1 (MQ=255)
         aaTGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAAtgc                                                                                       >  1:853058/1‑100 (MQ=255)
         aaTGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAAtgc                                                                                       <  2:1001933/100‑1 (MQ=255)
             gTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCg                                                                                   >  1:82883/1‑100 (MQ=255)
              tCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCgg                                                                                  >  1:1420194/1‑100 (MQ=255)
              tCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTAATAGAACGTTTGCAATAAGGGCTTCCCATGTTTACATAAGTTGCACAAGCCAAtgctgc                                                                                    >  2:385212/1‑98 (MQ=255)
              tCGCCAGTCGCCGGTAGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCgg                                                                                  <  1:1130926/100‑1 (MQ=255)
               cGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGt                                                                                 <  2:1914305/100‑1 (MQ=255)
                gCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGGTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGt                                                                                 <  2:667684/99‑1 (MQ=255)
                  cAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTAt                                                                              >  1:2214186/1‑100 (MQ=255)
                   aGTCGCCGGTGGCAATAAAAAACTGGCTCCTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTCACATCAGTGGCGCAAGCCAATGCTGCGGTTATc                                                                             <  2:1435935/100‑1 (MQ=255)
                    ctCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCg                                                                            >  2:700184/2‑100 (MQ=255)
                     tCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGa                                                                           <  1:1887835/100‑1 (MQ=255)
                      cGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGaa                                                                          >  1:1316848/1‑100 (MQ=255)
                       gCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGaa                                                                          >  1:924148/1‑99 (MQ=255)
                       gCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAAc                                                                         >  1:148982/1‑100 (MQ=255)
                       gCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAAc                                                                         >  1:520271/1‑100 (MQ=255)
                             ggCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATtc                                                                   >  2:940651/1‑100 (MQ=255)
                             ggCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATtc                                                                   <  2:520268/100‑1 (MQ=255)
                               cAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATtcgt                                                                 >  1:719137/1‑100 (MQ=255)
                               cAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATtcgt                                                                 >  1:794769/1‑100 (MQ=255)
                               cAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATtc                                                                   >  1:770255/1‑98 (MQ=255)
                               cAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATtcgt                                                                 >  2:1027652/1‑100 (MQ=255)
                                  tAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTcgc                                                              <  1:364408/100‑1 (MQ=255)
                                        aCTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGt                                                        <  1:1116676/100‑1 (MQ=255)
                                        aCTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGt                                                        <  1:1027663/100‑1 (MQ=255)
                                         cTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGt                                                        <  2:1150685/99‑1 (MQ=255)
                                          tGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTcc                                                      <  1:1446313/100‑1 (MQ=255)
                                            gCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTcc                                                      <  2:1169097/98‑1 (MQ=255)
                                              tCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCacac                                                  <  1:24575/100‑1 (MQ=255)
                                                 tttATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACActg                                                >  2:2108408/1‑99 (MQ=255)
                                                    attaGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACActggctg                                            <  2:1366802/100‑1 (MQ=255)
                                                    attaGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACActggctg                                            <  1:2338141/100‑1 (MQ=255)
                                                     ttaGAACGTTTGCAATAAGGTGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTTCGGTTATCGAACAAATTCGTCGCGCTCGTCCACActggctgg                                           >  2:1551660/1‑100 (MQ=255)
                                                     ttaGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACActggctgg                                           <  1:1942051/100‑1 (MQ=255)
                                                     ttaGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACActggctgg                                           >  1:1718599/1‑100 (MQ=255)
                                                       aGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCTCTCGTCCACActggctg                                            >  1:746741/1‑97 (MQ=255)
                                                        gAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGAtg                                        >  2:215371/1‑100 (MQ=255)
                                                        gAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGAtg                                        <  2:719129/100‑1 (MQ=255)
                                                        gAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGAtg                                        >  1:1035362/1‑100 (MQ=255)
                                                          aCGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGAtgtg                                      <  2:768782/100‑1 (MQ=255)
                                                          aCGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGAtgtg                                      >  2:2389656/1‑100 (MQ=255)
                                                           cGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGAtgtg                                      <  1:74107/99‑1 (MQ=255)
                                                            gTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCa                                    >  1:861137/1‑100 (MQ=255)
                                                            gTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCa                                    >  1:41484/1‑100 (MQ=255)
                                                            gTTTGCAATAAGGGGTTCCCATGTTTACAACAGTGGCACAAACCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGACTAGAtg                                        >  2:2169342/1‑96 (MQ=255)
                                                              ttGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCaa                                   >  2:603257/1‑99 (MQ=255)
                                                               tGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAAcc                                 <  1:1478749/100‑1 (MQ=255)
                                                                gCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCg                                >  1:1924940/1‑100 (MQ=255)
                                                                 cAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCgg                               >  1:118576/1‑100 (MQ=255)
                                                                  aaTAAGGGGTTCCGATGTTTAAATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGc                              <  1:1551678/100‑1 (MQ=255)
                                                                  aaTAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGc                              <  1:1992549/100‑1 (MQ=255)
                                                                  aaTAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGc                              >  2:1235015/1‑100 (MQ=255)
                                                                    tAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGctt                            >  1:1101463/1‑100 (MQ=255)
                                                                       ggggTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGcttctt                         >  1:1307033/1‑100 (MQ=255)
                                                                        gggTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCttcttc                        <  2:82881/100‑1 (MQ=255)
                                                                         ggTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCa                       <  2:668890/100‑1 (MQ=255)
                                                                         ggTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCa                       <  2:882627/100‑1 (MQ=255)
                                                                         ggTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCttcttc                        >  1:2118078/1‑99 (MQ=255)
                                                                          gTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCAc                      <  1:1359963/100‑1 (MQ=255)
                                                                          gTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCAc                      <  1:1610433/100‑1 (MQ=255)
                                                                          gTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCAc                      >  1:1586889/1‑100 (MQ=255)
                                                                           ttCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACt                     <  1:900366/100‑1 (MQ=255)
                                                                           ttCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACt                     <  1:1422464/100‑1 (MQ=255)
                                                                               cATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATc                 <  1:1326934/100‑1 (MQ=255)
                                                                                aTGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATc                 <  2:1404872/99‑1 (MQ=255)
                                                                                aTGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCa                >  1:1631158/1‑100 (MQ=255)
                                                                                   tttACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCg             <  2:558153/100‑1 (MQ=255)
                                                                                   tttACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGc              >  1:834685/1‑99 (MQ=255)
                                                                                      aCATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAAc          <  2:1751046/100‑1 (MQ=255)
                                                                                      aCATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAAc          <  1:1581731/100‑1 (MQ=255)
                                                                                      aCATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAAc          <  2:1101453/100‑1 (MQ=255)
                                                                                       cATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACt         <  1:915425/100‑1 (MQ=255)
                                                                                          cAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTaaa      >  2:1049032/1‑100 (MQ=255)
                                                                                              ggCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAg  <  2:1439080/100‑1 (MQ=255)
                                                                                                 |                                                                                                
TGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAG  >  NC_000913/548277‑548470

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