Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,871,779 A→G D79D (GAT→GAC cysC ← adenosine 5'‑phosphosulfate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,871,7790AG86.9% 720.5 / 95.5 297D79D (GAT→GACcysCadenosine 5'‑phosphosulfate kinase
Reads supporting (aligned to +/- strand):  ref base A (23/16);  new base G (136/122);  total (159/138)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.64e-01

Only 100 of 298 total aligned reads displayed.

GGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTC  >  NC_000913/2871687‑2871868
                                                                                            |                                                                                         
ggCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCa                                                                                           <  1:6628219/93‑1 (MQ=255)
        aaaTGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTaaa                                                                                   <  2:2740494/93‑1 (MQ=255)
         aaTGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTaaaa                                                                                  <  1:108080/93‑1 (MQ=255)
         aaTGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTaaaa                                                                                  <  2:3348208/93‑1 (MQ=255)
          aTGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAAc                                                                                 <  1:34954/93‑1 (MQ=255)
          aTGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAAc                                                                                 >  1:4539071/1‑93 (MQ=255)
          aTGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAAc                                                                                 <  2:5909860/93‑1 (MQ=255)
           tGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAAcc                                                                                <  2:2933261/93‑1 (MQ=255)
             cGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCga                                                                              <  2:1599701/93‑1 (MQ=255)
               gTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCgaga                                                                            >  1:4790430/1‑93 (MQ=255)
                tCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGAt                                                                           <  1:6452688/93‑1 (MQ=255)
                tCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGAt                                                                           <  2:482877/93‑1 (MQ=255)
                   gCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGc                                                                        <  1:6293463/93‑1 (MQ=255)
                   gCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGc                                                                        <  2:6084734/93‑1 (MQ=255)
                   gCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGc                                                                        >  2:6532737/1‑93 (MQ=255)
                      ccaccaGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGc                                                                     >  1:3439869/1‑93 (MQ=255)
                      ccaccaGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGc                                                                     >  2:123899/1‑93 (MQ=255)
                       caccaGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCa                                                                    <  2:3098264/93‑1 (MQ=255)
                         ccaGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATa                                                                  >  2:365469/1‑93 (MQ=255)
                         ccaGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATa                                                                  <  1:4300247/93‑1 (MQ=255)
                         ccaGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATa                                                                  >  1:3635671/1‑93 (MQ=255)
                          caGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATaa                                                                 <  1:4661953/93‑1 (MQ=255)
                            gTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATc                                                               <  1:3943155/93‑1 (MQ=255)
                            gTCCGACTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATc                                                               >  1:4700637/1‑93 (MQ=255)
                             tCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATcc                                                              <  1:5916765/93‑1 (MQ=255)
                             tCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATcc                                                              <  1:1801697/93‑1 (MQ=255)
                               cGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGt                                                            <  1:1567701/93‑1 (MQ=255)
                                ggCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTg                                                           <  1:3723574/93‑1 (MQ=255)
                                ggCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTg                                                           <  1:1200349/93‑1 (MQ=255)
                                ggCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTg                                                           >  2:602364/1‑93 (MQ=255)
                                 gCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTgg                                                          >  1:6506555/1‑93 (MQ=255)
                                  cTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGc                                                         <  1:3027746/93‑1 (MQ=255)
                                   ttCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCg                                                        >  2:7119130/1‑93 (MQ=255)
                                   ttCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCg                                                        >  2:3580821/1‑93 (MQ=255)
                                    tCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGa                                                       <  2:5756661/93‑1 (MQ=255)
                                     cAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGaa                                                      >  2:3894617/1‑93 (MQ=255)
                                     cAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGaa                                                      <  1:956578/93‑1 (MQ=255)
                                      aaCCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAAc                                                     >  2:1160639/1‑93 (MQ=255)
                                       aCCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACa                                                    <  1:1446100/93‑1 (MQ=255)
                                       aCCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACa                                                    >  2:6191985/1‑93 (MQ=255)
                                         cATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACAtt                                                  <  1:3343433/93‑1 (MQ=255)
                                         cATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACAtt                                                  >  1:7489626/1‑93 (MQ=255)
                                           tCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGt                                                >  1:6770603/1‑93 (MQ=255)
                                           tCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGt                                                >  1:208194/1‑93 (MQ=255)
                                           tCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGt                                                <  2:5227705/93‑1 (MQ=255)
                                            cAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTc                                               <  1:1679518/93‑1 (MQ=255)
                                            cAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTc                                               <  2:2191436/93‑1 (MQ=255)
                                                    ccACTTCACCGACGCGACGGATATTCTGTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCca                                       >  2:6415721/1‑93 (MQ=255)
                                                    ccACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCca                                       >  2:3498248/1‑93 (MQ=255)
                                                    ccACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCca                                       <  1:1958711/93‑1 (MQ=255)
                                                    ccACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCca                                       >  1:1476840/1‑93 (MQ=255)
                                                     cACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcag                                      >  2:1401250/1‑93 (MQ=255)
                                                     cACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCCGCATAATCCGTGGCGAACATTGTCGCCATCcag                                      >  2:6931920/1‑93 (MQ=255)
                                                        ttCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcagcag                                   >  1:5048239/1‑93 (MQ=255)
                                                        ttCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcagcag                                   >  2:673068/1‑93 (MQ=255)
                                                         tCACCGACGCGACGGATATTCTCTTTATGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGa                                  <  2:5837942/93‑1 (MQ=255)
                                                         tCACCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGa                                  <  1:4102298/93‑1 (MQ=255)
                                                           aCCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa                                >  1:1018405/1‑93 (MQ=255)
                                                           aCCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa                                <  1:2784526/93‑1 (MQ=255)
                                                           aCCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa                                <  1:1495960/93‑1 (MQ=255)
                                                           aCCGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa                                <  2:262062/93‑1 (MQ=255)
                                                            ccGACGCGACGGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATAc                               >  1:4565904/1‑93 (MQ=255)
                                                             cgacgcgacgGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACg                              <  2:1140049/93‑1 (MQ=255)
                                                               acgcgacgGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTa                            <  1:5795497/93‑1 (MQ=255)
                                                                cgcgacgGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTAc                           <  2:3764078/93‑1 (MQ=255)
                                                                 gcgacgGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACt                          >  1:5310974/1‑93 (MQ=255)
                                                                  cgacgGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTg                         <  1:6195524/93‑1 (MQ=255)
                                                                   gacgGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGa                        <  2:105495/93‑1 (MQ=255)
                                                                   gacgGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGa                        <  2:4954682/93‑1 (MQ=255)
                                                                    acgGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGAc                       <  1:717237/93‑1 (MQ=255)
                                                                    acgGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGAc                       >  1:275964/1‑93 (MQ=255)
                                                                    acgGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGAc                       >  2:1885419/1‑93 (MQ=255)
                                                                    acgGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGAc                       >  2:1491580/1‑93 (MQ=255)
                                                                    acgGATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGAc                       >  1:3184415/1‑93 (MQ=255)
                                                                       gATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGcc                    >  2:363244/1‑93 (MQ=255)
                                                                       gATATTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGcc                    <  1:5716836/93‑1 (MQ=255)
                                                                        atatTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCg                   <  1:7025847/93‑1 (MQ=255)
                                                                        atatTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCg                   >  2:1277401/1‑93 (MQ=255)
                                                                         tatTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGa                  >  1:5473176/1‑93 (MQ=255)
                                                                         tatTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGa                  >  2:7178786/1‑93 (MQ=255)
                                                                          atTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAg                 >  1:5594932/1‑93 (MQ=255)
                                                                          atTCTCTTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAg                 <  1:2628389/93‑1 (MQ=255)
                                                                              tctTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTa             >  2:2493949/1‑93 (MQ=255)
                                                                               ctTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTAt            >  2:751533/1‑93 (MQ=255)
                                                                               ctTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTAt            <  2:5103745/93‑1 (MQ=255)
                                                                               ctTTACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTAt            >  2:4543070/1‑93 (MQ=255)
                                                                                tttACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATg           <  1:85177/93‑1 (MQ=255)
                                                                                 ttACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGt          >  1:2259760/1‑93 (MQ=255)
                                                                                 ttACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGt          >  2:6972630/1‑93 (MQ=255)
                                                                                 ttACGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGt          <  2:3217808/93‑1 (MQ=255)
                                                                                  tACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTa         <  2:5642451/93‑1 (MQ=255)
                                                                                    cGATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAAc       <  2:4430543/93‑1 (MQ=255)
                                                                                     gATCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACg      >  2:5021228/1‑93 (MQ=255)
                                                                                      aTCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGc     >  2:2573429/1‑93 (MQ=255)
                                                                                      aTCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGc     >  2:7347948/1‑93 (MQ=255)
                                                                                      aTCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGc     >  2:5419280/1‑93 (MQ=255)
                                                                                       tCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGcc    >  2:4616651/1‑93 (MQ=255)
                                                                                       tCGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGcc    <  2:5783641/93‑1 (MQ=255)
                                                                                       tCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGcc    >  2:234235/1‑93 (MQ=255)
                                                                                        cGGCGTCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGcct   <  1:738945/93‑1 (MQ=255)
                                                                                            |                                                                                         
GGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTC  >  NC_000913/2871687‑2871868

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