Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,440,519 A→G E187E (GAA→GAG ytfM → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,440,5190AG100.0% 743.8 / NA 230E187E (GAA→GAGytfMhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (105/125);  total (105/125)

Only 100 of 234 total aligned reads displayed.

TTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATGA  >  NC_000913/4440429‑4440608
                                                                                          |                                                                                         
ttATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCg                                                                                         >  1:716867/1‑93 (MQ=255)
  aTTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGtt                                                                                       <  1:1945052/93‑1 (MQ=255)
     tCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTAcc                                                                                    >  2:3081693/1‑93 (MQ=255)
         tAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCtt                                                                                >  2:5635773/1‑93 (MQ=255)
          aGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCttt                                                                               <  2:552153/93‑1 (MQ=255)
          aGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCttt                                                                               <  2:3587998/93‑1 (MQ=255)
           gCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTg                                                                              >  2:1991025/1‑93 (MQ=255)
           gCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTg                                                                              >  2:120517/1‑93 (MQ=255)
             gAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTggg                                                                            >  1:5661259/1‑93 (MQ=255)
             gAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTggg                                                                            >  2:2175951/1‑93 (MQ=255)
              aaTTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGc                                                                           >  1:1087957/1‑93 (MQ=255)
              aaTTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGc                                                                           <  1:2688862/93‑1 (MQ=255)
               aTTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCa                                                                          <  2:1610897/93‑1 (MQ=255)
               aTTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCa                                                                          >  1:4866928/1‑93 (MQ=255)
                tttACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCAt                                                                         <  1:5574217/93‑1 (MQ=255)
                 ttACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCAtg                                                                        <  2:550057/93‑1 (MQ=255)
                 ttACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCAtg                                                                        <  1:5930298/93‑1 (MQ=255)
                   aCCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCAtgtg                                                                      <  1:6290374/93‑1 (MQ=255)
                   aCCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCAtgtg                                                                      >  1:4786584/1‑93 (MQ=255)
                    ccAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGa                                                                     <  2:2139723/93‑1 (MQ=255)
                      aaaGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGAcc                                                                   >  2:1082315/1‑93 (MQ=255)
                       aaGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCt                                                                  <  2:5281224/93‑1 (MQ=255)
                       aaGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCt                                                                  >  1:2062125/1‑93 (MQ=255)
                       aaGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCt                                                                  >  1:4613995/1‑93 (MQ=255)
                        aGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCtt                                                                 <  1:6173536/93‑1 (MQ=255)
                        aGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCtt                                                                 >  1:4026334/1‑93 (MQ=255)
                         gcgcAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCttt                                                                <  2:356218/93‑1 (MQ=255)
                         gcgcAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCttt                                                                >  2:4689336/1‑93 (MQ=255)
                           gcAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGa                                                              >  2:1250082/1‑93 (MQ=255)
                           gcAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGa                                                              <  1:2660802/93‑1 (MQ=255)
                           gcAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGa                                                              >  2:859195/1‑93 (MQ=255)
                            cAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGaa                                                             >  1:2531692/1‑93 (MQ=255)
                            cAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGaa                                                             <  2:4925541/93‑1 (MQ=255)
                            cAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGaa                                                             <  2:2355790/93‑1 (MQ=255)
                             aGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAg                                                            <  2:1505128/93‑1 (MQ=255)
                              gCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAgg                                                           >  1:4663771/1‑93 (MQ=255)
                              gCTGGGCAATGCGCTCGGCTTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAgg                                                           >  2:6165030/1‑93 (MQ=255)
                                 gggCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATc                                                        >  2:4754775/1‑93 (MQ=255)
                                 gggCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATc                                                        <  2:3936650/93‑1 (MQ=255)
                                 gggCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATc                                                        <  1:1481970/93‑1 (MQ=255)
                                 gggCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATc                                                        <  2:2701514/93‑1 (MQ=255)
                                 gggCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATc                                                        <  2:5798508/93‑1 (MQ=255)
                                  ggCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATca                                                       <  2:1718007/93‑1 (MQ=255)
                                  ggCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATca                                                       <  2:2272277/93‑1 (MQ=255)
                                   gCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATcac                                                      <  1:1580163/93‑1 (MQ=255)
                                    cATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATcaca                                                     <  2:3695296/93‑1 (MQ=255)
                                      ttGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACaaa                                                   <  1:3528330/93‑1 (MQ=255)
                                       tGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAAt                                                  <  2:587091/93‑1 (MQ=255)
                                           cTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCcgc                                              >  1:1714874/1‑93 (MQ=255)
                                            tCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCcgcg                                             >  1:5960900/1‑93 (MQ=255)
                                            tCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCcgcg                                             >  1:575486/1‑93 (MQ=255)
                                              ggCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGAt                                           >  1:3928365/1‑93 (MQ=255)
                                                 cTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGaa                                        >  1:3034130/1‑93 (MQ=255)
                                                 cTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGaa                                        <  2:4653070/93‑1 (MQ=255)
                                                  tGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAAt                                       >  1:2122086/1‑93 (MQ=255)
                                                  tGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAAt                                       >  1:704193/1‑93 (MQ=255)
                                                  tGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAAt                                       <  2:588487/93‑1 (MQ=255)
                                                   gCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATa                                      >  2:1864334/1‑93 (MQ=255)
                                                   gCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATa                                      <  2:5301121/93‑1 (MQ=255)
                                                    cATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATAc                                     <  2:4547692/93‑1 (MQ=255)
                                                    cATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATAc                                     <  1:5046835/93‑1 (MQ=255)
                                                     aTAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATAcc                                    <  1:1059010/93‑1 (MQ=255)
                                                     aTAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATAcc                                    <  1:5314812/93‑1 (MQ=255)
                                                     aTAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATAcc                                    <  2:3197791/93‑1 (MQ=255)
                                                      tAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCt                                   >  2:482467/1‑93 (MQ=255)
                                                       aaaGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTg                                  >  1:3755367/1‑93 (MQ=255)
                                                       aaaGCCTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTg                                  >  1:931925/1‑93 (MQ=255)
                                                            cTTCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCaaaa                             >  1:4526463/1‑93 (MQ=255)
                                                             ttCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAAt                            >  2:3393815/1‑93 (MQ=255)
                                                              tCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATc                           >  2:2499704/1‑93 (MQ=255)
                                                              tCTGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATc                           <  2:5715440/93‑1 (MQ=255)
                                                                tGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTg                         <  1:1109226/93‑1 (MQ=255)
                                                                tGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTg                         <  2:4219334/93‑1 (MQ=255)
                                                                tGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTg                         >  1:2970953/1‑93 (MQ=255)
                                                                tGGGATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTg                         >  2:4087322/1‑93 (MQ=255)
                                                                 gggATATTGATTATAACAGTGGCGAGCGTTACCGCTTTTGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTgg                        >  2:2680053/1‑93 (MQ=255)
                                                                 gggATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTgg                        <  2:733095/93‑1 (MQ=255)
                                                                 gggATATTGATTACAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGGCAAATCTgg                        <  2:6151668/93‑1 (MQ=255)
                                                                  ggATATTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGt                       >  1:793666/1‑93 (MQ=255)
                                                                    atatTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGc                     >  1:491359/1‑93 (MQ=255)
                                                                    atatTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGc                     <  1:3302760/93‑1 (MQ=255)
                                                                    atatTGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGc                     >  1:3813598/1‑93 (MQ=255)
                                                                       ttGATTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGt                  >  2:2816726/1‑93 (MQ=255)
                                                                          aTTATAACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTa               >  1:2882685/1‑93 (MQ=255)
                                                                            tataACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTaaa             <  2:4132683/93‑1 (MQ=255)
                                                                             ataACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGGGCCGTTTAAag            <  1:5797683/93‑1 (MQ=255)
                                                                              taACAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAaga           <  1:3412272/93‑1 (MQ=255)
                                                                               aaCAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAagag          <  1:1669841/93‑1 (MQ=255)
                                                                               aaCAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAagag          <  2:4052048/93‑1 (MQ=255)
                                                                               aaCAGTGGCGAGCGTGACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAagag          <  1:1285646/93‑1 (MQ=255)
                                                                                 cAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAggg        <  2:2031783/93‑1 (MQ=255)
                                                                                 cAGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAggg        <  1:4159629/93‑1 (MQ=255)
                                                                                  aGTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAgggg       >  1:3765445/1‑92 (MQ=255)
                                                                                   gTGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCg      >  2:508625/1‑93 (MQ=255)
                                                                                    tGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGa     <  2:1033662/93‑1 (MQ=255)
                                                                                    tGGCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGa     >  1:4291024/1‑93 (MQ=255)
                                                                                     ggCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGAt    <  2:5994406/93‑1 (MQ=255)
                                                                                     ggCGAGCGTGACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGAt    <  2:2578995/93‑1 (MQ=255)
                                                                                      gCGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATg   <  1:1125216/93‑1 (MQ=255)
                                                                                       cGAGCGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATGa  <  1:4446677/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATGA  >  NC_000913/4440429‑4440608

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