Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,511,306 A→T 7.0% L51H (CTT→CAT)  fecE ← ferric citrate ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,511,3060AT 7.0% 193.5 / 5.2 57L51H (CTT→CAT) fecEferric citrate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (28/25);  new base T (2/2);  total (30/27)
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.78e-01

ACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAA  >  NC_000913/4511158‑4511452
                                                                                                                                                    |                                                                                                                                                  
aCCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTaaaa                                                                                                                                                    <  3:466704/149‑1 (MQ=255)
    tACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCcg                                                                                                                                                <  4:221917/149‑1 (MQ=255)
       gAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGa                                                                                                                                             >  1:526699/1‑149 (MQ=255)
           ccAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCTTTTAATGCcgcagtctg                                                                                                                                         <  2:552876/149‑7 (MQ=255)
             aGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACa                                                                                                                                       <  1:484491/149‑1 (MQ=255)
              gCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAg                                                                                                                                      >  2:516887/1‑149 (MQ=255)
              gCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAg                                                                                                                                      <  3:379848/149‑1 (MQ=255)
                tCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGtt                                                                                                                                    >  3:539453/1‑149 (MQ=255)
                  cTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTa                                                                                                                                  >  1:182497/1‑149 (MQ=255)
                       cTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGc                                                                                                                             >  4:127469/1‑149 (MQ=255)
                                 ctctGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTcc                                                                                                                   <  1:398544/149‑1 (MQ=255)
                                 ctctGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTcc                                                                                                                   >  4:198332/1‑149 (MQ=255)
                                 ctctGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTcc                                                                                                                   >  1:307098/1‑149 (MQ=255)
                                  tctGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTccc                                                                                                                  <  1:285301/149‑1 (MQ=255)
                                      gCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCaa                                                                                                              <  2:307098/149‑1 (MQ=255)
                                        gTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAAcc                                                                                                            >  2:45012/1‑149 (MQ=255)
                                        gTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAAcc                                                                                                            >  3:81404/1‑149 (MQ=255)
                                         ttAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCg                                                                                                           >  3:238183/1‑149 (MQ=255)
                                          tAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGt                                                                                                          <  1:45012/149‑1 (MQ=255)
                                             aTGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGTCGCGGAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAg                                                                                                       >  4:45127/1‑149 (MQ=255)
                                              tGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAgg                                                                                                      >  2:589/1‑149 (MQ=255)
                                               ggTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGa                                                                                                     <  4:539453/149‑1 (MQ=255)
                                                  gCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCg                                                                                                  >  1:574141/1‑149 (MQ=255)
                                                    tGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGAt                                                                                                >  1:294369/1‑149 (MQ=255)
                                                      aGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTATAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCa                                                                                              <  3:236647/149‑1 (MQ=255)
                                                      aGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCa                                                                                              <  3:198332/149‑1 (MQ=255)
                                                              gcGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTg                                                                                      <  3:524336/149‑1 (MQ=255)
                                                               cGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGa                                                                                     <  1:391311/149‑1 (MQ=255)
                                                                          gggCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGtt                                                                          >  4:214543/1‑149 (MQ=255)
                                                                            gCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTa                                                                                                        >  3:7496/1‑117 (MQ=255)
                                                                            gCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTa                                                                                                        <  4:7496/117‑1 (MQ=255)
                                                                             ccAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGc                                                                       <  2:238562/149‑1 (MQ=255)
                                                                                        gATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGa                                                            >  1:339127/1‑149 (MQ=255)
                                                                                           gagaGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAAc                                                         <  3:127469/149‑1 (MQ=255)
                                                                                              agCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAAcgtc                                                      <  4:121099/149‑1 (MQ=255)
                                                                                                cATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAAcgtcgt                                                    <  2:116070/149‑1 (MQ=255)
                                                                                                    tttATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAg                                                >  3:285504/1‑149 (MQ=255)
                                                                                                       aTGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTAc                                             <  1:365837/149‑1 (MQ=255)
                                                                                                         gggATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCt                                           >  4:246834/1‑149 (MQ=255)
                                                                                                          ggATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCtt                                          <  1:450153/149‑1 (MQ=255)
                                                                                                           gATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTtg                                         >  4:16180/1‑149 (MQ=255)
                                                                                                               aTCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTtgtctg                                     >  1:493016/1‑149 (MQ=255)
                                                                                                                tCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCTGCGTAAAACAGTTTAACAGCGTAGATTTCCCGCAACAGCTAGGACCGATCAGGGCGGTGATCTTCCCAGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTtgtctgt                                    >  3:361386/1‑149 (MQ=255)
                                                                                                                 cGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTtgtctgtc                                   <  3:246834/149‑1 (MQ=255)
                                                                                                                           tACGGTGCCAGTCTGCGGCATTAAAAGCCGTGAACAACAGTTTAACAGCGTCGATTGCCCGCAACCGTTAGGACCGATCTGGGCGGTGATCTTCCCCGTAGGCAGTGGGAGTGAAACGTCGTTAAGTACCTCGTCTGTCCCGTAACTGa                         <  4:361386/149‑1 (MQ=255)
                                                                                                                           tACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGa                         <  2:493016/149‑1 (MQ=255)
                                                                                                                           tACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGa                         >  2:282482/1‑149 (MQ=255)
                                                                                                                            aCGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGAc                        <  2:574141/149‑1 (MQ=255)
                                                                                                                            aCGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGAc                        >  4:330569/1‑149 (MQ=255)
                                                                                                                              ggTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCg                      >  1:22786/1‑149 (MQ=255)
                                                                                                                                  gtttttcgCGGCTTTTAATGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAg                  >  2:4754/8‑149 (MQ=255)
                                                                                                                                    ttttcgCGGCTTTTAATGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGAt                >  4:105338/6‑149 (MQ=255)
                                                                                                                                    aGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGAt                >  1:180740/1‑149 (MQ=255)
                                                                                                                                        tGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCa                   >  2:75005/1‑142 (MQ=255)
                                                                                                                                        tGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCa                   <  1:75005/142‑1 (MQ=255)
                                                                                                                                         gCGGCTTTTAATGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCa           <  3:105338/149‑1 (MQ=255)
                                                                                                                                              aTTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGGCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGACTTTCAGTTCg      >  1:114028/1‑149 (MQ=255)
                                                                                                                                               ttAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGt     <  1:495181/149‑1 (MQ=255)
                                                                                                                                                  aaaGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTaaa  <  2:542318/149‑1 (MQ=255)
                                                                                                                                                  aaaGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTaaa  <  4:91127/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAA  >  NC_000913/4511158‑4511452

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