Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,040,218 G→A A127T (GCA→ACA)  appA → phosphoanhydride phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,040,2180GA99.2% 438.6 / ‑5.4 127A127T (GCA→ACA) appAphosphoanhydride phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (57/69);  total (58/69)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 130 total aligned reads displayed.

AAGGGCTGCCCGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGGCACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGTTTGCCAA  >  NC_000913/1040128‑1040307
                                                                                          |                                                                                         
  gggCTGCCCGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACAcc                                                                                       <  1:755982/93‑1 (MQ=255)
  gggCTGCCCGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACAcc                                                                                       >  2:575549/1‑93 (MQ=255)
  gggCTGCCCGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACAcc                                                                                       <  2:4645515/93‑1 (MQ=255)
  gggCTGCCCGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACAcc                                                                                       <  1:2379864/93‑1 (MQ=255)
    gCTGCCCGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTg                                                                                     <  2:3173889/93‑1 (MQ=255)
     cTGCCCGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGGCGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGa                                                                                    <  1:900275/93‑1 (MQ=255)
     cTGCCCGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGa                                                                                    <  2:1775285/93‑1 (MQ=255)
       gCCCGCAGTCGGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACt                                                                                  <  1:972598/93‑1 (MQ=255)
        cccGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACtg                                                                                 <  1:1331745/93‑1 (MQ=255)
           gCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGc                                                                              >  1:1151569/1‑93 (MQ=255)
           gCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGc                                                                              >  2:1948080/1‑93 (MQ=255)
             aGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCaa                                                                            >  2:1034758/1‑93 (MQ=255)
              gTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAAt                                                                           >  2:3703085/1‑93 (MQ=255)
              gTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAAt                                                                           >  2:2791944/1‑93 (MQ=255)
              gTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAAt                                                                           >  1:787723/1‑93 (MQ=255)
               tCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATa                                                                          >  1:3377455/1‑93 (MQ=255)
                cTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATaa                                                                         >  2:2380207/1‑93 (MQ=255)
                cTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATaa                                                                         <  2:4138986/93‑1 (MQ=255)
                 tGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAAc                                                                        >  1:598483/1‑93 (MQ=255)
                 tGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAAc                                                                        >  2:3276157/1‑93 (MQ=255)
                 tGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAAc                                                                        >  1:3562273/1‑93 (MQ=255)
                 tGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGACGCCGGGCTGACACCTGACTGTGCAATAAc                                                                        >  2:341830/1‑93 (MQ=255)
                  ggTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAAcc                                                                       <  1:4657487/93‑1 (MQ=255)
                    tCAGGTCGCGATTATTGCTTATGTCGACGAGCGTACCCGTAAAATAGGCGAAGTCTTCGCCGTTTTGCTGACACTTGACTGTGCAATTACCGt                                                                     >  2:4679836/1‑93 (MQ=255)
                    tCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGt                                                                     >  2:234746/1‑93 (MQ=255)
                     cAGGTCGCTATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTa                                                                    <  2:1551148/93‑1 (MQ=255)
                     cAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTa                                                                    <  1:709520/93‑1 (MQ=255)
                     cAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTa                                                                    <  2:2248702/93‑1 (MQ=255)
                      aGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTAc                                                                   >  1:951991/1‑93 (MQ=255)
                      aGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTAc                                                                   <  2:1411190/93‑1 (MQ=255)
                       ggTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGGCACCTGACTGTGCAATAACCGTACa                                                                  >  1:2291917/1‑93 (MQ=255)
                       ggTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACa                                                                  >  2:897609/1‑93 (MQ=255)
                        gTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGGACAt                                                                 >  1:2904250/1‑93 (MQ=255)
                         tCGCGGTTATTGGTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATa                                                                <  2:1086762/93‑1 (MQ=255)
                         tCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGACTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATa                                                                >  1:2843793/1‑93 (MQ=255)
                          cgcgATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATAc                                                               <  1:2153278/93‑1 (MQ=255)
                            cgATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATAccc                                                             <  1:4275046/93‑1 (MQ=255)
                             tattatgGCTGATGTCGCCGAGCGGTCCAGTAAAGCAGGCTAAGCCTTCGCCG‑CGGGCTGACACCTGACTGTGCAATAACCGTACATCCCCAg                                                           <  1:1582304/92‑1 (MQ=255)
                             gATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCa                                                            <  1:3638166/93‑1 (MQ=255)
                              attattGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAg                                                           >  1:2227965/1‑93 (MQ=255)
                              attagtGCTGATGTCGACGAGCGTTCCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTTCAATAACCGTACATACCCAg                                                           <  1:1234423/93‑1 (MQ=255)
                                tattGCTGATGTCGCCGAGCTTCCCCGTAAAACGGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGc                                                         <  1:6433/93‑1 (MQ=255)
                                tattGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGc                                                         >  1:3566357/1‑93 (MQ=255)
                                 tttGCTGAGGTCGACGAGCGTCCCCGAAAAACAGGCGAGGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCa                                                        <  2:4040769/92‑1 (MQ=255)
                                 attGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCa                                                        <  1:274862/93‑1 (MQ=255)
                                 attGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCa                                                        >  1:3788398/1‑93 (MQ=255)
                                 attGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCa                                                        >  2:1372744/1‑93 (MQ=255)
                                   tGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACCCCTGACTGTGCAATAACCGTACATACCCAGGCAGa                                                      <  1:1516466/93‑1 (MQ=255)
                                   tGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGa                                                      <  1:2019762/93‑1 (MQ=255)
                                   tGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGa                                                      <  1:658026/93‑1 (MQ=255)
                                    gCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGAt                                                     >  2:3831269/1‑93 (MQ=255)
                                     cTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATa                                                    <  1:2663306/93‑1 (MQ=255)
                                       gATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACg                                                  <  1:3064973/93‑1 (MQ=255)
                                         tGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTc                                                >  2:3646403/1‑93 (MQ=255)
                                          gTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCACATACGtaa                                               >  2:227682/1‑91 (MQ=255)
                                            cgacgaGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAg                                             <  1:3787929/93‑1 (MQ=255)
                                            cgacgaGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAg                                             <  1:2980007/93‑1 (MQ=255)
                                              acgaGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTc                                           <  1:3858665/93‑1 (MQ=255)
                                               cgaGCGTCCCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGTTACGTCCAGTcc                                          <  1:2533520/93‑1 (MQ=255)
                                                     tccccGTAAAACAGGCGACCCTTGCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATcc                                    <  2:3544960/91‑1 (MQ=255)
                                                      aCCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCg                                   <  1:3214768/93‑1 (MQ=255)
                                                       cccGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGt                                  <  1:1101594/93‑1 (MQ=255)
                                                        ccGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGtt                                 >  2:4074953/1‑93 (MQ=255)
                                                        ccGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGtt                                 >  1:88669/1‑93 (MQ=255)
                                                        ccGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGtt                                 <  1:2623057/93‑1 (MQ=255)
                                                         cGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTa                                >  2:3533527/1‑93 (MQ=255)
                                                         cGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTa                                <  1:4209655/93‑1 (MQ=255)
                                                          gTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTAt                               <  2:2998688/93‑1 (MQ=255)
                                                          gTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCATACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTAt                               <  1:3406013/93‑1 (MQ=255)
                                                           tAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTAtt                              >  1:2040656/1‑93 (MQ=255)
                                                            aaaaCAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTAttt                             >  1:2241826/1‑93 (MQ=255)
                                                              aaCAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTaa                           >  1:3627782/1‑93 (MQ=255)
                                                              aaCAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTaa                           >  1:3492686/1‑93 (MQ=255)
                                                               aCAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAAt                          <  2:3290077/93‑1 (MQ=255)
                                                               aCAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAAt                          >  2:3552213/1‑93 (MQ=255)
                                                               aCAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAAt                          >  1:3981507/1‑93 (MQ=255)
                                                                cAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATc                         <  2:1562909/93‑1 (MQ=255)
                                                                cAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATc                         >  1:3155632/1‑93 (MQ=255)
                                                                 aGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATcc                        <  2:1420952/93‑1 (MQ=255)
                                                                  ggCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCct                       >  2:568272/1‑93 (MQ=255)
                                                                   gCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCctc                      <  1:1798832/93‑1 (MQ=255)
                                                                    cGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCctct                     <  2:4610752/93‑1 (MQ=255)
                                                                    cGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCctct                     <  1:4363863/93‑1 (MQ=255)
                                                                      aaGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTaa                   <  1:4588188/93‑1 (MQ=255)
                                                                      aaGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTaa                   <  1:4620178/93‑1 (MQ=255)
                                                                       aGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTaaa                  <  2:4621841/93‑1 (MQ=255)
                                                                        gCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTaaaa                 >  2:1052916/1‑93 (MQ=255)
                                                                        gCCGTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTaaaa                 <  2:980178/93‑1 (MQ=255)
                                                                         ccTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTaaaaa                >  1:3564848/1‑93 (MQ=255)
                                                                          cTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAAc               <  2:670614/93‑1 (MQ=255)
                                                                           ttCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACt              >  1:4585539/1‑93 (MQ=255)
                                                                            tCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTg             <  2:2741487/93‑1 (MQ=255)
                                                                              gccgccGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATTCTCTAAAGACTggg           >  2:1659441/1‑92 (MQ=255)
                                                                                cgccgGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGt         <  2:441782/93‑1 (MQ=255)
                                                                                 gccgGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGtt        <  1:4038044/93‑1 (MQ=255)
                                                                                  ccgGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGttt       <  2:1946803/93‑1 (MQ=255)
                                                                                    gggCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGTTTGc     <  2:4490416/93‑1 (MQ=255)
                                                                                      gCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGTTTGCCa   >  2:4439690/1‑93 (MQ=255)
                                                                                      gCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGTTTGCCa   <  2:3948878/93‑1 (MQ=255)
                                                                                      gCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGTTTGCCa   <  2:966372/93‑1 (MQ=255)
                                                                                          |                                                                                         
AAGGGCTGCCCGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGGCACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGTTTGCCAA  >  NC_000913/1040128‑1040307

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑