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,1540TC86.9% 626.6 / 80.3 252S83S (AGT→AGCrstBsensor histidine protein kinase
Reads supporting (aligned to +/- strand):  ref base T (23/10);  new base C (115/104);  total (138/114)
Fisher's exact test for biased strand distribution p-value = 9.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-02

Only 100 of 253 total aligned reads displayed.

GAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTAT  >  NC_000913/1681062‑1681243
                                                                                            |                                                                                         
gAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGt                                                                                           >  1:592791/1‑93 (MQ=255)
 aaTTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTa                                                                                          >  2:6148389/1‑93 (MQ=255)
        tGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCa                                                                                   <  2:6701539/93‑1 (MQ=255)
        tGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCa                                                                                   >  1:4127914/1‑93 (MQ=255)
        tGAGACCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGTTCTGCGTGTCGAGCCACTGAGCAAATACCa                                                                                   <  1:5593115/93‑1 (MQ=255)
            aTCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                               <  1:6807641/93‑1 (MQ=255)
             tCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTtg                                                                              >  2:2101412/1‑93 (MQ=255)
              ccccccACACGACTGGGGTAAAACTCTGAAAGAGATGGAGTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTtga                                                                             >  1:7501799/1‑93 (MQ=255)
                 cccACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTtgatga                                                                          >  1:7107088/1‑93 (MQ=255)
                 cccACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTtgatga                                                                          >  1:6408912/1‑93 (MQ=255)
                  ccACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTgatgat                                                                         >  1:1134942/1‑93 (MQ=255)
                  ccACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTgatgat                                                                         >  2:810039/1‑93 (MQ=255)
                  ccACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTgatgat                                                                         <  2:7520902/93‑1 (MQ=255)
                  ccACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTgatgat                                                                         >  1:618860/1‑93 (MQ=255)
                   cacaCGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGata                                                                        <  2:691469/93‑1 (MQ=255)
                    acacGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGatat                                                                       <  1:4009031/93‑1 (MQ=255)
                      acGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATAttt                                                                     <  2:4141629/93‑1 (MQ=255)
                       cGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTc                                                                    <  2:3715353/93‑1 (MQ=255)
                       cGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTc                                                                    <  2:1439619/93‑1 (MQ=255)
                          cTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCAt                                                                 <  2:4034313/93‑1 (MQ=255)
                          cTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCAt                                                                 >  1:3446567/1‑93 (MQ=255)
                           tGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATg                                                                >  1:7473839/1‑93 (MQ=255)
                            ggggTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGc                                                               >  1:4303730/1‑93 (MQ=255)
                            ggggTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGc                                                               >  1:252574/1‑93 (MQ=255)
                            ggggTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGc                                                               <  2:6296607/93‑1 (MQ=255)
                             gggTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCa                                                              >  1:4733395/1‑93 (MQ=255)
                             gggTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCa                                                              <  2:1779508/93‑1 (MQ=255)
                              ggTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCAc                                                             <  2:2427892/93‑1 (MQ=255)
                               gTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCAcc                                                            <  2:465211/93‑1 (MQ=255)
                                  aaaCTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGAc                                                         >  2:5592450/1‑93 (MQ=255)
                                  aaaCTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGAc                                                         >  1:5525341/1‑93 (MQ=255)
                                   aaCTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACt                                                        <  1:165200/93‑1 (MQ=255)
                                    aCTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTg                                                       <  2:3796374/93‑1 (MQ=255)
                                     ctctGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGc                                                      <  1:1990279/93‑1 (MQ=255)
                                     ctctGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGc                                                      >  1:2323844/1‑93 (MQ=255)
                                      tctGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCg                                                     >  1:6791194/1‑93 (MQ=255)
                                      tctGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCg                                                     >  2:1294867/1‑93 (MQ=255)
                                       ctGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGt                                                    <  1:5367897/93‑1 (MQ=255)
                                       ctGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGt                                                    <  2:6079251/93‑1 (MQ=255)
                                          aaaGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTggc                                                 <  2:1063847/93‑1 (MQ=255)
                                          aaaGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTggc                                                 >  1:5215923/1‑93 (MQ=255)
                                          aaaGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTggc                                                 >  2:2890911/1‑93 (MQ=255)
                                           aaGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTggcg                                                >  1:3919531/1‑93 (MQ=255)
                                             gagaTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTggcggc                                              >  2:3111468/1‑93 (MQ=255)
                                               gaTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGa                                            >  1:4524282/1‑93 (MQ=255)
                                                aTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGaa                                           >  2:2994484/1‑93 (MQ=255)
                                                aTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGaa                                           <  1:2494618/93‑1 (MQ=255)
                                                 tGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGaaa                                          >  2:3945114/1‑93 (MQ=255)
                                                 tGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGaaa                                          >  1:5686528/1‑93 (MQ=255)
                                                   gATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAAtt                                        <  1:3544593/93‑1 (MQ=255)
                                                    aTTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTg                                       >  1:48091/1‑93 (MQ=255)
                                                     tttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGt                                      <  2:6539987/93‑1 (MQ=255)
                                                     tttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGt                                      <  2:6327417/93‑1 (MQ=255)
                                                     tttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGt                                      >  1:716678/1‑93 (MQ=255)
                                                     tttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGt                                      >  1:2308099/1‑93 (MQ=255)
                                                       tAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCg                                    >  1:6643510/1‑93 (MQ=255)
                                                        aaaTCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGc                                   <  2:4784796/93‑1 (MQ=255)
                                                        aaaTCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGc                                   <  2:5072314/93‑1 (MQ=255)
                                                         aaTCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGcc                                  <  2:5936780/93‑1 (MQ=255)
                                                         aaTCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGcc                                  <  1:1926221/93‑1 (MQ=255)
                                                          aTCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGccc                                 >  1:4238574/1‑93 (MQ=255)
                                                          aTCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGccc                                 >  1:2104338/1‑93 (MQ=255)
                                                           tctctcTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCt                                <  1:882563/93‑1 (MQ=255)
                                                           tctctcTTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCt                                >  1:1052097/1‑93 (MQ=255)
                                                            ctctctTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTg                               >  2:7492574/1‑93 (MQ=255)
                                                             tctctTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTgg                              >  2:7173680/1‑93 (MQ=255)
                                                             tctctTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTgg                              >  2:4315685/1‑93 (MQ=255)
                                                               tctTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGAc                            <  1:58254/93‑1 (MQ=255)
                                                                ctTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACg                           >  2:5658133/1‑93 (MQ=255)
                                                                ctTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACg                           >  2:5437081/1‑93 (MQ=255)
                                                                ctTTCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACg                           >  1:6279474/1‑93 (MQ=255)
                                                                 tttCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGa                          >  2:3539334/1‑93 (MQ=255)
                                                                  ttCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGAt                         >  1:6053336/1‑93 (MQ=255)
                                                                   tCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATc                        <  1:3201562/93‑1 (MQ=255)
                                                                   tCGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGACATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATc                        <  1:7299950/93‑1 (MQ=255)
                                                                    cGATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCa                       <  1:529188/93‑1 (MQ=255)
                                                                     gATCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAg                      >  1:285508/1‑93 (MQ=255)
                                                                      aTCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGt                     >  2:6330614/1‑93 (MQ=255)
                                                                      aTCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGt                     <  1:2583362/93‑1 (MQ=255)
                                                                      aTCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGt                     >  1:5158646/1‑93 (MQ=255)
                                                                       tCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTa                    >  1:6335725/1‑93 (MQ=255)
                                                                       tCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTa                    >  1:2146501/1‑93 (MQ=255)
                                                                       tCTGCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTa                    <  1:4114346/93‑1 (MQ=255)
                                                                         tGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTaca                  >  2:3307376/1‑93 (MQ=255)
                                                                          gCGTGTCGAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTacac                 >  1:1313335/1‑93 (MQ=255)
                                                                            gtgtCGAGCCCCTTAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTTGACGATCAGGACACgg               >  1:3656160/1‑92 (MQ=255)
                                                                                 gAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTg          >  2:1604589/1‑93 (MQ=255)
                                                                                 gAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTg          >  2:1983340/1‑93 (MQ=255)
                                                                                 gAGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTg          >  1:4610179/1‑93 (MQ=255)
                                                                                  aGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGc         <  1:6900983/93‑1 (MQ=255)
                                                                                  aGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGc         <  2:526329/93‑1 (MQ=255)
                                                                                  aGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGc         <  2:4964727/93‑1 (MQ=255)
                                                                                  aGCCACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGc         >  1:3280589/1‑93 (MQ=255)
                                                                                   gCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCa        >  2:3411158/1‑93 (MQ=255)
                                                                                     cACTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGc      <  2:7204150/93‑1 (MQ=255)
                                                                                      aCTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCg     >  2:6110447/1‑93 (MQ=255)
                                                                                      aCTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCg     >  2:226242/1‑93 (MQ=255)
                                                                                      aCTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCg     >  2:70117/1‑93 (MQ=255)
                                                                                       cTGAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGt    >  1:6841615/1‑93 (MQ=255)
                                                                                         gAGCAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTAt  >  2:3181147/1‑93 (MQ=255)
                                                                                            |                                                                                         
GAATTGCGTGAGATCCCCCCACACGACTGGGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTGCAGCGTAT  >  NC_000913/1681062‑1681243

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