Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,511,303 A→T 7.0% L52* (TTA→TAA)  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,3030AT 7.0% 192.2 / 4.8 57L52* (TTA→TAA) 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 = 9.89e-01

ACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGT  >  NC_000913/4511155‑4511449
                                                                                                                                                    |                                                                                                                                                  
aCGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTa                                                                                                                                                    >  1:542318/1‑149 (MQ=255)
  gACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTaaa                                                                                                                                                  >  1:116070/1‑149 (MQ=255)
  gACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTaaa                                                                                                                                                  >  4:466704/1‑149 (MQ=255)
   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:307098/1‑149 (MQ=255)
                                    ctctGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTcc                                                                                                                <  1:398544/149‑1 (MQ=255)
                                    ctctGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTcc                                                                                                                >  4:198332/1‑149 (MQ=255)
                                     tctGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTccc                                                                                                               <  1:285301/149‑1 (MQ=255)
                                         gCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCaa                                                                                                           <  2:307098/149‑1 (MQ=255)
                                           gTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAAcc                                                                                                         >  3:81404/1‑149 (MQ=255)
                                           gTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAAcc                                                                                                         >  2:45012/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:282482/1‑149 (MQ=255)
                                                                                                                              tACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGa                      <  2:493016/149‑1 (MQ=255)
                                                                                                                               aCGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGAc                     >  4:330569/1‑149 (MQ=255)
                                                                                                                               aCGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGAc                     <  2:574141/149‑1 (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                <  1:75005/142‑1 (MQ=255)
                                                                                                                                           tGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCa                >  2:75005/1‑142 (MQ=255)
                                                                                                                                            gCGGCTTTTAATGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCa        <  3:105338/149‑1 (MQ=255)
                                                                                                                                                 aTTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGGCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGACTTTCAGTTCg   >  1:114028/1‑149 (MQ=255)
                                                                                                                                                  ttAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGt  <  1:495181/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGT  >  NC_000913/4511155‑4511449

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