Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 144,977 (A)7→8 coding (401/441 nt) yadI → putative PTS enzyme IIA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913144,9701.A99.1% 810.1 / ‑0.9 220coding (394/441 nt)yadIputative PTS enzyme IIA component
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base A (99/119);  total (99/121)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

Only 100 of 267 total aligned reads displayed.

TGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATAAAAGCGTT  >  NC_000913/144878‑145057
                                                                                             |                                                                                       
tgGCTTGCCGTGAAACCATGACCAGTTCGGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACcaacaac                                                                                          >  2:3678608/1‑93 (MQ=255)
tgGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACcaacaac                                                                                          <  2:5853763/93‑1 (MQ=255)
tgGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACcaacaac                                                                                          <  2:3487500/93‑1 (MQ=255)
tgGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACcaacaac                                                                                          >  1:2028533/1‑93 (MQ=255)
 ggCTTGCCGTGAAACCATGACCAGTTCAGAGTTGCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACcaacaaca                                                                                        <  1:7617079/93‑1 (MQ=255)
  gCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACcaacaacaa                                                                                       <  2:6081640/93‑1 (MQ=255)
   gttGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACaaa                                                                                      <  1:3521772/92‑1 (MQ=255)
    ttGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaa                                                                                     >  1:935044/1‑93 (MQ=255)
     tGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaaa                                                                                    <  2:5873683/93‑1 (MQ=255)
       ccGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaaaaa                                                                                  <  1:7598869/93‑1 (MQ=255)
        cGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGGGAGTGGTCTTTGGCACCAACAACaaaaaaag                                                                                 >  1:2490602/1‑92 (MQ=255)
         gTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaaaaaat                                                                                >  1:2655090/1‑91 (MQ=255)
          tGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaaaaaatc                                                                               >  1:5524141/1‑90 (MQ=255)
           gAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaaaaaatcc                                                                              >  1:899958/1‑89 (MQ=255)
           gAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaaaaaatcc                                                                              >  1:7119955/1‑89 (MQ=255)
           gAAACCATGACCAGTTCAGAGTTTCGCGAGCGGATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaaaaaatcc                                                                              <  1:1089008/93‑5 (MQ=255)
             aaCCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGc                                                                             >  2:6761933/1‑93 (MQ=255)
             aaCCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGc                                                                             <  1:4082495/93‑1 (MQ=255)
               ccATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCt                                                                           <  1:783420/93‑1 (MQ=255)
               ccATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCt                                                                           <  1:1444760/93‑1 (MQ=255)
                cATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCtt                                                                          <  1:5838484/93‑1 (MQ=255)
                 aTGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCttt                                                                         <  2:2688478/93‑1 (MQ=255)
                  tGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTc                                                                        <  2:5437180/93‑1 (MQ=255)
                   gACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTAGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCg                                                                       >  2:7433786/1‑93 (MQ=255)
                        gTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCta                                                                  >  1:1666551/1‑92 (MQ=255)
                        gTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTc                                                                  <  1:3399549/93‑1 (MQ=255)
                        gTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGGGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTc                                                                  >  2:1588432/1‑93 (MQ=255)
                         ttCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCa                                                                 <  1:4873000/93‑1 (MQ=255)
                         ttCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCa                                                                 >  2:2298605/1‑93 (MQ=255)
                         ttCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCa                                                                 >  1:3125541/1‑93 (MQ=255)
                         ttCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCa                                                                 >  1:6194329/1‑93 (MQ=255)
                         ttCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTAGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCa                                                                 >  1:6086925/1‑93 (MQ=255)
                          tCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGCGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCaa                                                                >  2:1645700/1‑93 (MQ=255)
                           cAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCaaa                                                               <  1:1138341/93‑1 (MQ=255)
                           cAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCaaa                                                               <  2:2051650/93‑1 (MQ=255)
                            agagTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAAc                                                              >  1:2958973/1‑93 (MQ=255)
                             gagTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACa                                                             <  2:2339704/93‑1 (MQ=255)
                                tttCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATaa                                                          <  2:1111833/93‑1 (MQ=255)
                                 ttCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAAt                                                         <  2:6953888/93‑1 (MQ=255)
                                   cgcgAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAAttt                                                       <  1:417043/93‑1 (MQ=255)
                                      gAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGtt                                                    >  1:6308307/1‑92 (MQ=255)
                                      gAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTa                                                    <  1:7160552/93‑1 (MQ=255)
                                      gAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGGAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTa                                                    >  1:4257322/1‑93 (MQ=255)
                                       aGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTAt                                                   <  1:5369625/93‑1 (MQ=255)
                                       aGCGTAATGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTAt                                                   <  1:6407091/93‑1 (MQ=255)
                                        gCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATg                                                  >  2:7288115/1‑93 (MQ=255)
                                         cGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGa                                                 <  1:1510725/93‑1 (MQ=255)
                                         cGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGa                                                 <  1:7918699/93‑1 (MQ=255)
                                          gTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAg                                                <  2:1206003/93‑1 (MQ=255)
                                          gTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAg                                                <  2:6562932/93‑1 (MQ=255)
                                          gTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAg                                                >  2:497920/1‑93 (MQ=255)
                                           tATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGt                                               <  2:6189047/93‑1 (MQ=255)
                                            aTTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTa                                              <  1:524775/93‑1 (MQ=255)
                                             ttGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTAt                                             >  1:4972591/1‑93 (MQ=255)
                                             ttGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTAt                                             <  2:1707776/93‑1 (MQ=255)
                                              tGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTAtt                                            >  2:5475067/1‑93 (MQ=255)
                                              tGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTAtt                                            >  2:1946948/1‑93 (MQ=255)
                                                tCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTaa                                          <  1:1918333/93‑1 (MQ=255)
                                                 cGAACTGGGTGCGCCGAAGGTGAGTAGTCTTTGGCAACAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAAc                                         <  2:4788705/93‑1 (MQ=255)
                                                    aCTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCcg                                      >  1:1027637/1‑93 (MQ=255)
                                                     cTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCcgc                                     <  2:932540/93‑1 (MQ=255)
                                                      tGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCcgcg                                    >  1:3587718/1‑93 (MQ=255)
                                                       gggTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGa                                   <  2:7768713/93‑1 (MQ=255)
                                                       gggTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGa                                   <  2:4042274/93‑1 (MQ=255)
                                                        ggTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGAt                                  <  2:1418897/93‑1 (MQ=255)
                                                        ggTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGAt                                  <  1:6380541/93‑1 (MQ=255)
                                                              ccGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGAt                            <  2:272271/93‑1 (MQ=255)
                                                              ccGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTACTAACCCGCGATTCTGAt                            <  1:1595224/93‑1 (MQ=255)
                                                               cGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATg                           <  1:3354112/93‑1 (MQ=255)
                                                                ggAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATgg                          <  1:6979590/93‑1 (MQ=255)
                                                                ggAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATgg                          <  2:4638777/93‑1 (MQ=255)
                                                                 gAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGgc                         >  1:6817195/1‑93 (MQ=255)
                                                                 gAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGgc                         >  1:6631378/1‑93 (MQ=255)
                                                                 gAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGgc                         <  1:7357601/93‑1 (MQ=255)
                                                                   ggTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGgcgc                       >  1:7697239/1‑93 (MQ=255)
                                                                   ggTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGgcgc                       <  2:2870739/93‑1 (MQ=255)
                                                                    gTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCt                      >  2:5558315/1‑93 (MQ=255)
                                                                     tGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCtt                     >  2:3020101/1‑93 (MQ=255)
                                                                     tGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCtt                     >  1:7246609/1‑93 (MQ=255)
                                                                     tGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCtt                     <  1:2297861/93‑1 (MQ=255)
                                                                      gAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCttt                    <  1:1881706/93‑1 (MQ=255)
                                                                       agtagtCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCtttt                   >  2:752050/1‑93 (MQ=255)
                                                                       agtagtCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCtttt                   >  2:1334026/1‑93 (MQ=255)
                                                                        gtagtCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTg                  <  2:7564516/93‑1 (MQ=255)
                                                                          agtCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCt                >  1:3297260/1‑93 (MQ=255)
                                                                           gtCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTa               >  1:679182/1‑93 (MQ=255)
                                                                           gtCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTa               >  1:3933842/1‑93 (MQ=255)
                                                                            tCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTAc              <  1:5197414/93‑1 (MQ=255)
                                                                               ttGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAAt           <  1:4613773/93‑1 (MQ=255)
                                                                                 ggCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATaa         >  1:3217254/1‑93 (MQ=255)
                                                                                  gCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATaaa        <  2:5618968/93‑1 (MQ=255)
                                                                                    aCCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATAAAAg      >  2:219185/1‑93 (MQ=255)
                                                                                     ccAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATAAAAGc     >  2:7249240/1‑93 (MQ=255)
                                                                                     ccAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATAAAAGc     >  2:4002724/1‑93 (MQ=255)
                                                                                     ccAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATAAAAGc     <  2:4264909/93‑1 (MQ=255)
                                                                                       aacaacaaAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATAAAAGCgt   <  1:5783360/93‑1 (MQ=255)
                                                                                        acaacaaAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATAAAAGCgtt  >  2:247443/1‑93 (MQ=255)
                                                                                        acaacaaAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATAAAAGCgtt  <  1:2231990/93‑1 (MQ=255)
                                                                                        acaacaaAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATAAAAGCgtt  <  2:5349519/93‑1 (MQ=255)
                                                                                        acaacaaAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATAAAAGCgtt  >  1:5974634/1‑93 (MQ=255)
                                                                                             |                                                                                       
TGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTCTGATGGCGCTTTTGCTACAATAAAAGCGTT  >  NC_000913/144878‑145057

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: