Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,260,888 C→T 9.0% A97V (GCC→GTC)  yjbM → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,260,8880CT 9.0% 191.7 / 6.9 66A97V (GCC→GTC) yjbMuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (31/29);  new base T (3/3);  total (34/32)
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.53e-01

CCCATTGAACAATGTACTTTTGATGTAGAAACTAAGAATGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAA  >  NC_000913/4260740‑4261033
                                                                                                                                                    |                                                                                                                                                 
cccATTGAACAATGTACTTTTGATGTAGAAACTAAGAATGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGATTATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTTGTAAtgc                                                                                                                                                   >  7:232131/1‑149 (MQ=255)
      gAACAATGTACTTTTGATGTAGAAACTAAGAATGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATTGGTTCTGGTTTTTTAACGTTTGTTAAtgcctgc                                                                                                                                               >  1:80862/1‑147 (MQ=255)
       aaCAATGTACTTTTGATGTAGAAACTAAGAATGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAGGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAg                                                                                                                                            >  7:213810/1‑149 (MQ=255)
       aaCAATGTACTTTTGATGTAGAAACTAAGAATGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAg                                                                                                                                            >  7:226758/1‑149 (MQ=255)
              tACTTTTGATGTAGAAACTAAGAATGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACatg                                                                                                                                     <  8:240426/149‑1 (MQ=255)
               aCTTTTGATGTAGAAACTAAGAATGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGATCCGTTAGCGAATATACATGTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACat                                                                                                                                      >  4:28211/1‑147 (MQ=255)
                cTTTTGATGTAGAAACTAAGAATGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACatgat                                                                                                                                   <  1:224258/149‑1 (MQ=255)
                cTTTTGATGTAGAAACTAAGAATGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACatga                                                                                                                                    <  5:144749/148‑1 (MQ=255)
                 ttttGATGTAGAAACTAAGAATGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACatgat                                                                                                                                   >  6:68557/1‑148 (MQ=255)
                 ttttGATGTAGAAACTAAGAATGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACatgat                                                                                                                                   >  5:12260/1‑148 (MQ=255)
                      aTGTAGAAACTAAGAATGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGacac                                                                                                                              >  7:194329/1‑147 (MQ=255)
                      aTGTAGAAACTAAGAATGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGacaa                                                                                                                              >  2:114281/1‑148 (MQ=255)
                         tAGAAACTAAGAATGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGc                                                                                                                          >  6:59795/1‑149 (MQ=255)
                          aGAAACTAAGAATGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTAAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAggt                                                                                                                         >  4:127801/1‑147 (MQ=255)
                           gAAACTAAGAATGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAGGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCt                                                                                                                         <  8:213810/148‑1 (MQ=255)
                            aaaCTAAGAATGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCtt                                                                                                                        >  2:116625/1‑148 (MQ=255)
                                tAAGAATGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCAAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAat                                                                                                                   >  3:118747/1‑149 (MQ=255)
                                 aaGAATGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAat                                                                                                                   >  1:172033/1‑148 (MQ=255)
                                    aaTGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAattatt                                                                                                               >  7:80045/1‑149 (MQ=255)
                                      tGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAAGCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTTGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTg                                                                                                             <  7:174552/149‑1 (MQ=255)
                                              ccTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCTAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATAATTTGCTTCGCt                                                                                                      <  4:118747/148‑1 (MQ=255)
                                              ccTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCt                                                                                                      <  1:116625/148‑1 (MQ=255)
                                               cTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGt                                                                                                    >  7:115641/1‑149 (MQ=255)
                                               cTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGt                                                                                                    >  7:160459/1‑149 (MQ=255)
                                                  cACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTaa                                                                                                 >  6:49893/1‑149 (MQ=255)
                                                     ggTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAAttt                                                                                              <  5:59795/149‑1 (MQ=255)
                                                          aGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACtt                                                                                         >  4:178580/1‑149 (MQ=255)
                                                                 ttATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCatat                                                                                   >  3:75717/1‑148 (MQ=255)
                                                                  tATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAAtgcct                                                                                                                                                 <  5:222834/85‑1 (MQ=255)
                                                                  tATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAAtgcct                                                                                                                                                 >  6:222834/1‑85 (MQ=255)
                                                                  tATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTc                                                                                 >  1:86113/1‑149 (MQ=255)
                                                                    tGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTa                                                                               <  1:13539/149‑1 (MQ=255)
                                                                     gAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAg                                                                              >  8:75140/1‑149 (MQ=255)
                                                                      aGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGt                                                                             <  8:115641/149‑1 (MQ=255)
                                                                              aaGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAattatt                                                                                                               >  7:141543/1‑107 (MQ=255)
                                                                              aaGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACg                                                                      >  1:117003/1‑148 (MQ=255)
                                                                               aGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATTGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTa                                                                        <  2:80862/145‑1 (MQ=255)
                                                                                 gATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCt                                                                                >  5:230698/1‑135 (MQ=255)
                                                                                 gATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCt                                                                                <  6:230698/135‑1 (MQ=255)
                                                                                 gATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACg                                                                      >  4:47539/1‑145 (MQ=255)
                                                                                 gATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGGTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGATCATGATGACAACAGCTTAATATTATTTGCTTCGCTGTTAATTTCACTGTCATATTCTAGTGCGTACg                                                                      >  8:93761/1‑145 (MQ=255)
                                                                                   tAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGAtt                                                                <  2:44237/149‑1 (MQ=255)
                                                                                      tGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTAt                                                             <  4:75717/149‑1 (MQ=255)
                                                                                        aaCCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGTCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCa                                                           >  6:8968/1‑149 (MQ=255)
                                                                                         aacTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAAAGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCa                                                           <  1:114281/146‑1 (MQ=255)
                                                                                         aCCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGTCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAg                                                          <  5:8968/149‑1 (MQ=255)
                                                                                         aCCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCCAGTGCCTACGCAGATTTATCAg                                                          <  4:115466/149‑1 (MQ=255)
                                                                                              tAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGaaa                                                       <  8:226758/147‑1 (MQ=255)
                                                                                                 cGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAgtg                                                  >  6:56195/1‑149 (MQ=255)
                                                                                                 cGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAgtg                                                  <  6:23150/149‑1 (MQ=255)
                                                                                                     tataCATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAGTGtata                                              <  5:56195/149‑1 (MQ=255)
                                                                                                        aCATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAGTGTATAtt                                            >  2:12139/1‑148 (MQ=255)
                                                                                                          aTTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAGGGTATATTAAt                                         >  7:84403/1‑149 (MQ=255)
                                                                                                           ttttCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAGTGTATATTAATg                                        <  7:44943/149‑1 (MQ=255)
                                                                                                                aaGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAGTGTATATTAATGaaaa                                    >  7:71983/1‑148 (MQ=255)
                                                                                                                aaGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAGTGTATATTAATGaaaa                                    <  8:151247/148‑1 (MQ=255)
                                                                                                                     aTAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATga                              <  3:178580/149‑1 (MQ=255)
                                                                                                                      tAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATgag                             <  5:68557/149‑1 (MQ=255)
                                                                                                                      tAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATgag                             <  8:80045/149‑1 (MQ=255)
                                                                                                                        ggTTCTGTGTTTTTAACGTTTGTAAATGCCTGCCATGAACATGATGACAACAGCTTAAAATTATTTGCATCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGAGTTATCAGAAACAGTGTATATTAATGAAAATAATgaga                            <  3:127801/148‑1 (MQ=255)
                                                                                                                               gTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAAGAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTAc                       <  3:28211/146‑1 (MQ=255)
                                                                                                                               gTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTAc                       <  3:134722/146‑1 (MQ=255)
                                                                                                                                         gTTTGTTAATGCCTTCAAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAGTGTATATTAATTAAAATAATGAGATCTAcc                      >  8:60181/1‑137 (MQ=255)
                                                                                                                                         gTTTGTTAATGCCTGCAAGTAACATGATGACAACATCTTAAAATTATTTGCTTCGATGTTAATTTCACTTTAATATTCTAGTGACTACGCAGATTTATCAGAAACATTGTATATTAATGAAAATAATGAGAGCTAcc                      <  7:60181/137‑1 (MQ=255)
                                                                                                                                                 aTGTCTGCCAGTAACATGATGACAACAGCTTAAAATTATTTTCTTCGCTTTTAATTTCACTTTCATATTATAGTGCCTACGCAGATTTATCAGAAACAGTGTATATTAATg                                        <  6:162743/111‑1 (MQ=255)
                                                                                                                                                 aTGTCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAGTGTATATTAATg                                        >  5:164347/1‑111 (MQ=255)
                                                                                                                                                 aTGTCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAGTGTATATTAATg                                        <  6:164347/111‑1 (MQ=255)
                                                                                                                                                 aTGTCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATAAGAAACAGTGTATATTAATg                                        >  5:162743/1‑111 (MQ=255)
                                                                                                                                                 aTGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTTATATTCTAGTGCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGaaaaa  >  4:98319/1‑149 (MQ=255)
                                                                                                                                                 aTGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTTCCTGAAAGCTCAGTTTGaaaaa  <  8:194329/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CCCATTGAACAATGTACTTTTGATGTAGAAACTAAGAATGACCTTCCCTCCACGGTTCAGCTATTTTATGAGGAATCTAAGGATAATGAACCTTTAGCGAATATACATTTTCAAGCAATAGGTTCTGGTTTTTTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAA  >  NC_000913/4260740‑4261033

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: ‑