Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,681,154 T→C S83S (AGT→AGC rstB → sensor histidine protein kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,681,1540TC99.7% 1079.3 / NA 310S83S (AGT→AGCrstBsensor histidine protein kinase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (164/145);  total (165/145)
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.99e-01

Only 100 of 318 total aligned reads displayed.

ATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTAT  >  NC_000913/1681064‑1681243
                                                                                          |                                                                                         
aTTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCaa                                                                                         <  2:4304485/93‑1 (MQ=255)
 ttGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAATTCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCaaa                                                                                        <  2:180028/93‑1 (MQ=255)
      tGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCa                                                                                   <  2:751604/93‑1 (MQ=255)
      tGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCa                                                                                   <  2:630711/93‑1 (MQ=255)
      tGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCa                                                                                   >  2:4838836/1‑93 (MQ=255)
          aTCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCtt                                                                               <  2:6971852/93‑1 (MQ=255)
             cccccACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTtgat                                                                            <  1:7034180/93‑1 (MQ=255)
              ccccACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTtgatg                                                                           <  2:8423157/93‑1 (MQ=255)
                 cacaCGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGata                                                                        >  2:4097785/1‑93 (MQ=255)
                   aacGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATAtt                                                                      <  2:2456685/92‑1 (MQ=255)
                      gACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTcc                                                                   <  2:2747233/93‑1 (MQ=255)
                      gACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTcc                                                                   <  2:3023535/93‑1 (MQ=255)
                       aCTGTGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCa                                                                  <  2:3759842/93‑1 (MQ=255)
                       aCTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCa                                                                  <  2:2113509/93‑1 (MQ=255)
                       aCTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCa                                                                  <  2:5159316/93‑1 (MQ=255)
                        cTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCAt                                                                 <  2:2257271/93‑1 (MQ=255)
                         tGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATg                                                                <  2:1643656/93‑1 (MQ=255)
                         tGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATg                                                                >  2:1307389/1‑93 (MQ=255)
                             gaaaaaCTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCAcc                                                            <  2:3139235/91‑1 (MQ=255)
                             gTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCAcc                                                            <  2:7241859/93‑1 (MQ=255)
                              tAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCg                                                           >  1:84163/1‑93 (MQ=255)
                               aaaaCTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGa                                                          >  1:336242/1‑93 (MQ=255)
                                aaaCTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGAc                                                         <  1:3803180/93‑1 (MQ=255)
                                aaaCTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGAc                                                         <  1:6489344/93‑1 (MQ=255)
                                aaaCTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGAc                                                         <  1:4979567/93‑1 (MQ=255)
                                 aaCTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACt                                                        <  1:3629925/93‑1 (MQ=255)
                                 aaCTCTGAAAGAGAGGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACt                                                        <  2:8692793/93‑1 (MQ=255)
                                    tctGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCg                                                     >  2:202749/1‑93 (MQ=255)
                                      tGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTg                                                   <  1:8189372/93‑1 (MQ=255)
                                       gAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACTATCTTGATGATATTTCCATGCACCGACTGCGTgg                                                  <  2:5498669/93‑1 (MQ=255)
                                       gAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTgg                                                  >  2:4028133/1‑93 (MQ=255)
                                       gAAAGAGATGGATTGAAATCTCTCTTTCGATCCGCGTGGCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTgg                                                  <  1:2375575/93‑1 (MQ=255)
                                        aaaGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTggg                                                 >  1:1888745/1‑92 (MQ=255)
                                        aaaGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTggc                                                 <  1:3218239/93‑1 (MQ=255)
                                        aaaGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTggc                                                 >  2:7184120/1‑93 (MQ=255)
                                        aaaGAGATGCATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTAATAAAATACCATCTTGCTGATATTTCCATGCTCCGACTGCGTggc                                                 >  2:5815999/1‑93 (MQ=255)
                                         aaGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTggcg                                                >  2:150510/1‑93 (MQ=255)
                                         aaGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTggcg                                                <  1:1216082/93‑1 (MQ=255)
                                          agagATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTggcgg                                               >  2:4165554/1‑93 (MQ=255)
                                          agagATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTggcgg                                               >  1:746117/1‑93 (MQ=255)
                                          agagATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTggcgg                                               >  1:2761687/1‑93 (MQ=255)
                                           gagaTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTggcggc                                              >  1:4924584/1‑93 (MQ=255)
                                            agaTTGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGgcggcg                                             >  2:5604439/1‑93 (MQ=255)
                                            agaTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGgcggcg                                             >  1:6396613/1‑93 (MQ=255)
                                            agaTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGgcggcg                                             <  1:8534530/93‑1 (MQ=255)
                                              aTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGaa                                           >  1:4796789/1‑93 (MQ=255)
                                              aTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGaa                                           >  1:7717691/1‑93 (MQ=255)
                                                ggATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAAt                                         >  1:106888/1‑93 (MQ=255)
                                                ggATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAAt                                         <  2:6079154/93‑1 (MQ=255)
                                                   tttAAATCGCTCTTTCTATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGt                                      <  1:3008261/93‑1 (MQ=255)
                                                    ttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTc                                     <  2:3447167/93‑1 (MQ=255)
                                                     tAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCg                                    <  1:6652858/93‑1 (MQ=255)
                                                     tAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCg                                    >  2:8040367/1‑93 (MQ=255)
                                                      aaaTCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGc                                   >  2:5624607/1‑93 (MQ=255)
                                                      aaaTCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGc                                   >  1:2620426/1‑93 (MQ=255)
                                                       aaTCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGcc                                  >  2:4781767/1‑93 (MQ=255)
                                                        aTCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGccc                                 >  1:406841/1‑93 (MQ=255)
                                                         tctctcTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCt                                >  2:3474336/1‑93 (MQ=255)
                                                          ctctttTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTg                               <  2:1955837/93‑1 (MQ=255)
                                                           tctctTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTgg                              >  1:4105701/1‑93 (MQ=255)
                                                            ctctTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGa                             <  2:1234697/93‑1 (MQ=255)
                                                             tctTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGAc                            <  1:2258816/93‑1 (MQ=255)
                                                             tctTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGAc                            >  1:1362646/1‑93 (MQ=255)
                                                             tctTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGAc                            >  1:522957/1‑93 (MQ=255)
                                                             tctTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGAc                            >  1:6932942/1‑93 (MQ=255)
                                                               tttCGATCTGCGTGTGGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGa                          <  2:6788771/93‑1 (MQ=255)
                                                                ttCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGTTGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCTCTCTGGACGAt                         >  2:7838719/1‑93 (MQ=255)
                                                                  cGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCa                       <  1:4158921/93‑1 (MQ=255)
                                                                  cGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCa                       >  2:3906080/1‑93 (MQ=255)
                                                                   gATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAg                      >  2:8540561/1‑93 (MQ=255)
                                                                    aTCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAga                     >  1:4377979/1‑92 (MQ=255)
                                                                    aTCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGt                     >  1:5517044/1‑93 (MQ=255)
                                                                    aTCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGt                     >  1:976941/1‑93 (MQ=255)
                                                                     tCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTa                    <  1:1549655/93‑1 (MQ=255)
                                                                     tCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTa                    >  2:8228990/1‑93 (MQ=255)
                                                                     tCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTa                    >  1:5431739/1‑93 (MQ=255)
                                                                      cTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTac                   <  1:5643084/93‑1 (MQ=255)
                                                                      cTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTac                   <  1:5781812/93‑1 (MQ=255)
                                                                      cTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTac                   <  2:5710468/93‑1 (MQ=255)
                                                                      cTGCGTGTCGAGCCACTAAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTac                   <  2:4521185/93‑1 (MQ=255)
                                                                       tGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTaca                  >  1:3072452/1‑93 (MQ=255)
                                                                        gCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTacac                 <  1:3450083/93‑1 (MQ=255)
                                                                          gtgtCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGt               >  2:176696/1‑93 (MQ=255)
                                                                           tgtCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGtt              >  1:6539381/1‑93 (MQ=255)
                                                                            gtCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGttt             >  2:4504464/1‑93 (MQ=255)
                                                                             tCGAGCCACTGAGCAAATACCATCTTGATTATATTTCCATGCTCCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGtttt            <  1:2105357/93‑1 (MQ=255)
                                                                             tCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGtttt            >  1:5121477/1‑93 (MQ=255)
                                                                               gAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTg          >  2:332728/1‑93 (MQ=255)
                                                                               gAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTg          <  1:3302177/93‑1 (MQ=255)
                                                                                 gCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCa        <  1:2433977/93‑1 (MQ=255)
                                                                                 gCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCa        >  1:6629195/1‑93 (MQ=255)
                                                                                 gCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCa        <  2:3543373/93‑1 (MQ=255)
                                                                                   cACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGc      >  1:3859107/1‑93 (MQ=255)
                                                                                   cACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGc      <  2:5943335/93‑1 (MQ=255)
                                                                                    aCTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCg     <  1:6124345/93‑1 (MQ=255)
                                                                                    aCTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCg     >  1:8643206/1‑93 (MQ=255)
                                                                                     cTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGt    <  2:6412231/93‑1 (MQ=255)
                                                                                      tGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTa   >  1:7224351/1‑93 (MQ=255)
                                                                                       gAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTAt  >  2:7040504/1‑93 (MQ=255)
                                                                                       gAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTAt  <  2:5281402/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTAT  >  NC_000913/1681064‑1681243

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: