Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,961,588 C→G R86T (AGG→ACG)  gppA ← guanosine pentaphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,961,5880CG99.8% 1613.2 / ‑5.2 500R86T (AGG→ACG) gppAguanosine pentaphosphatase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (267/232);  total (268/232)
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.95e-01

Only 100 of 518 total aligned reads displayed.

TCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACCTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATT  >  NC_000913/3961498‑3961677
                                                                                          |                                                                                         
tCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTc                                                                                         >  1:1205916/1‑93 (MQ=255)
tCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTc                                                                                         >  1:8890073/1‑93 (MQ=255)
tCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTc                                                                                         <  1:9627079/93‑1 (MQ=255)
tCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTc                                                                                         >  2:129470/1‑93 (MQ=255)
     ttCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGccct                                                                                    >  2:7196099/1‑91 (MQ=255)
      tCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTa                                                                                   >  1:4601113/1‑93 (MQ=255)
          cGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGa                                                                               <  2:4938324/93‑1 (MQ=255)
            cTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACa                                                                             >  2:9081857/1‑93 (MQ=255)
            cTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACa                                                                             >  1:472018/1‑93 (MQ=255)
             tGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACaa                                                                            >  2:3677269/1‑93 (MQ=255)
             tGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACaa                                                                            >  1:3454274/1‑93 (MQ=255)
              gATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCGCCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAc                                                                           <  2:1181111/93‑1 (MQ=255)
              gATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAc                                                                           >  2:453555/1‑93 (MQ=255)
              gATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAc                                                                           >  1:9985949/1‑93 (MQ=255)
               aTCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAca                                                                          <  2:5416724/93‑2 (MQ=255)
               aTCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAca                                                                          >  1:9065774/1‑92 (MQ=255)
                tCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcac                                                                         >  2:499534/1‑91 (MQ=255)
                tCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcac                                                                         <  1:4420705/93‑3 (MQ=255)
                 cACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcacg                                                                        >  1:9893420/1‑93 (MQ=255)
                  aCCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGa                                                                       >  1:2884185/1‑93 (MQ=255)
                   ccTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGaa                                                                      >  2:9581175/1‑93 (MQ=255)
                   ccTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGaa                                                                      <  1:135615/93‑1 (MQ=255)
                   ccTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGaa                                                                      <  2:6430807/93‑1 (MQ=255)
                     tGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAAtt                                                                    >  2:5936259/1‑93 (MQ=255)
                     tGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAAtt                                                                    <  2:7919044/93‑1 (MQ=255)
                      gTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAAttt                                                                   <  2:1508332/93‑1 (MQ=255)
                       tACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTg                                                                  >  1:7095833/1‑93 (MQ=255)
                        aCCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGc                                                                 <  1:8278418/93‑1 (MQ=255)
                        aCCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGc                                                                 <  1:1708152/93‑1 (MQ=255)
                           ggACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgag                                                              <  2:683761/93‑1 (MQ=255)
                           ggACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgag                                                              <  2:6684342/93‑1 (MQ=255)
                            gACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCCCCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagg                                                             <  1:311853/93‑1 (MQ=255)
                            gACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagg                                                             >  1:5185991/1‑93 (MQ=255)
                             aCAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgaggg                                                            <  1:4892801/93‑1 (MQ=255)
                               aaCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagggag                                                          >  2:7473303/1‑93 (MQ=255)
                               aaCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagggag                                                          <  1:3240468/93‑1 (MQ=255)
                                 ccGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGggagga                                                        >  2:4425253/1‑93 (MQ=255)
                                  cGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGaa                                                       <  1:7417573/93‑1 (MQ=255)
                                  cGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGaa                                                       >  1:3661763/1‑93 (MQ=255)
                                  cGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGaa                                                       <  1:858308/93‑1 (MQ=255)
                                   gagaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAat                                                      >  1:4173269/1‑93 (MQ=255)
                                   gagaATTTCCTGCGCTTGTGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAat                                                      <  1:7928682/93‑1 (MQ=255)
                                    agaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAata                                                     <  1:9331594/93‑1 (MQ=255)
                                     gaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAatat                                                    <  2:6199911/93‑1 (MQ=255)
                                     gaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAatat                                                    >  2:2833755/1‑93 (MQ=255)
                                      aaTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATc                                                   <  1:4523401/93‑1 (MQ=255)
                                      aaTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATc                                                   >  1:3953200/1‑93 (MQ=255)
                                       aTTTCCTGCGCTTTGGCAATAAAATCCCCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCt                                                  <  2:7758754/93‑1 (MQ=255)
                                       aTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCt                                                  <  1:8125408/93‑1 (MQ=255)
                                       aTTTCCTGCGCTGTGGCAATAAAATCCCCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCt                                                  <  1:1947826/93‑1 (MQ=255)
                                           ccTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCa                                              <  2:5763066/93‑1 (MQ=255)
                                           ccTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCa                                              <  2:873617/93‑1 (MQ=255)
                                             tGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGa                                            >  2:2014974/1‑93 (MQ=255)
                                             tGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGa                                            >  1:2129043/1‑93 (MQ=255)
                                             tGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGa                                            <  2:8049897/93‑1 (MQ=255)
                                                 cTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGagctc                                        >  2:1450543/1‑89 (MQ=255)
                                                 cTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGtt                                        <  2:7126691/93‑1 (MQ=255)
                                                 cTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGtt                                        >  1:7780447/1‑93 (MQ=255)
                                                  tttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTc                                       >  1:10008069/1‑93 (MQ=255)
                                                   ttGGCAATAAAATCTCCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGATAACACGAATTTGCGAGGCAGGAATATCTTGCAGACGTTCa                                      <  1:1779190/93‑1 (MQ=255)
                                                      gCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCa                                   <  2:408793/93‑1 (MQ=255)
                                                        aaTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCaaa                                 >  1:2417327/1‑93 (MQ=255)
                                                         aTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAAc                                <  1:6938872/93‑1 (MQ=255)
                                                         aTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAAc                                >  1:7578946/1‑93 (MQ=255)
                                                          tAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa                               <  2:8122502/93‑1 (MQ=255)
                                                          tAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa                               >  1:6102559/1‑93 (MQ=255)
                                                          tAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa                               >  2:2484349/1‑93 (MQ=255)
                                                           aaaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAg                              <  1:5866822/93‑1 (MQ=255)
                                                            aaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAgg                             <  2:3481948/93‑1 (MQ=255)
                                                            aaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAgg                             >  2:6026255/1‑93 (MQ=255)
                                                             aaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgc                            >  2:1868360/1‑93 (MQ=255)
                                                             aaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgc                            <  1:2716149/93‑1 (MQ=255)
                                                             aaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgc                            <  1:3764098/93‑1 (MQ=255)
                                                              aTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgcg                           >  2:3218693/1‑93 (MQ=255)
                                                                  gccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaaa                       >  1:2537475/2‑93 (MQ=255)
                                                                  cccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaaa                       >  1:6794390/1‑93 (MQ=255)
                                                                  cccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaaa                       >  1:6626023/1‑93 (MQ=255)
                                                                   ccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAAc                      <  1:1493614/93‑1 (MQ=255)
                                                                   ccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAAc                      >  2:37978/1‑93 (MQ=255)
                                                                    cGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACa                     >  2:266207/1‑93 (MQ=255)
                                                                     gCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAt                    >  1:9490558/1‑93 (MQ=255)
                                                                     gCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAt                    <  1:8385713/93‑1 (MQ=255)
                                                                      cATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAtt                   <  1:6337928/93‑1 (MQ=255)
                                                                       aTTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTg                  >  2:2353323/1‑93 (MQ=255)
                                                                        ttGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGc                 <  1:5497629/93‑1 (MQ=255)
                                                                         tGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGcc                >  2:5321375/1‑93 (MQ=255)
                                                                         tGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGcc                >  1:3011760/1‑93 (MQ=255)
                                                                          gACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCa               >  2:1757235/1‑93 (MQ=255)
                                                                            cGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAAc             <  2:9783662/93‑1 (MQ=255)
                                                                             ggCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAAcc            <  1:5935987/93‑1 (MQ=255)
                                                                              gCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACcg           >  2:8920188/1‑93 (MQ=255)
                                                                               cAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACcgc          >  2:55733/1‑93 (MQ=255)
                                                                               cAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACcgc          <  1:5649315/93‑1 (MQ=255)
                                                                                   gCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTc      >  2:3095251/1‑93 (MQ=255)
                                                                                     gTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCa    >  2:8517980/1‑93 (MQ=255)
                                                                                     gTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCa    <  2:8385591/93‑1 (MQ=255)
                                                                                     gTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCa    <  2:9299357/93‑1 (MQ=255)
                                                                                     gTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCa    >  2:2217948/1‑93 (MQ=255)
                                                                                      tAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCAt   <  1:9114369/93‑1 (MQ=255)
                                                                                      tAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCAt   >  2:1673532/1‑93 (MQ=255)
                                                                                          |                                                                                         
TCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACCTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATT  >  NC_000913/3961498‑3961677

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