Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,495,653 C→A 5.9% L28M (CTG→ATG) ‡ mdeA → methionine gamma‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,495,6530CA 5.9% 244.0 / 3.5 77L28M (CTG→ATG) ‡mdeAmethionine gamma‑lyase
Reads supporting (aligned to +/- strand):  ref base C (36/36);  new base A (3/2);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

TATTTTGGAAAAAATCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATGGCTTCG  >  NC_002947/1495507‑1495799
                                                                                                                                                     |                                                                                                                                                  
tATTTTGGAAAAAATCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTg                                                                                                                                                  >  2:129276/1‑149 (MQ=255)
tATTTTGGAAAAAATCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTg                                                                                                                                                  >  6:7370/1‑149 (MQ=255)
tATTTTGGAAAAAATCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTg                                                                                                                                                  >  2:202802/1‑149 (MQ=255)
  ttttGGAAAAAATCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAATCACAGGCTGACCCCAACGGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGccc                                                                                                                                                    <  1:334894/148‑1 (MQ=255)
     tGGAAAAAATCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGc                                                                                                                                              <  7:350753/148‑1 (MQ=255)
       gAAAAAATCGGGCAGCGCGGTGGCATAAGCTTATAGCAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCa                                                                                                                                            <  4:379343/148‑1 (MQ=255)
           aaaTCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCACCAGGGCGCCAccatg                                                                                                                                       <  7:299320/149‑4 (MQ=255)
           aaaTCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCACCAGGGCGCCAccatg                                                                                                                                       <  5:409247/149‑4 (MQ=255)
              tCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACA‑‑‑AGACGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTa                                                                                                                                    <  1:233836/149‑1 (MQ=255)
                gtgCAGCTCGGTGGCATAAGCTTATAGCAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCAACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTAc                                                                                                                                   <  7:12824/146‑1 (MQ=255)
                gggCAGCTCGGTGGCATAAGCTTATAGCAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTAc                                                                                                                                   <  1:176484/148‑1 (MQ=255)
                 ggCAGCTCGGTGGCATAAGCTTATAGCAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCa                                                                                                                                 <  8:83066/149‑1 (MQ=255)
                    aGCTCGGTGGCATAAGCTTATAGCAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGAc                                                                                                                              >  3:49079/1‑149 (MQ=255)
                       tCGGTGGCATAAGCTTATAGCAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTAcc                                                                                                                                  <  5:248809/142‑1 (MQ=255)
                       tCGGTGGCATAAGCTTATAGCAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTAcc                                                                                                                                  >  6:248809/1‑142 (MQ=255)
                           tGGCATAAGCTTATAGCAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTAcc                                                                                                                       >  3:16296/1‑149 (MQ=255)
                              cATAAGCTTATAGCAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTAt                                                                                                                    <  5:302450/149‑1 (MQ=255)
                                  aGCTTATAGCAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGcc                                                                                                                 <  4:49079/148‑1 (MQ=255)
                                    cTTATAGCAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGTGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTc                                                                                                              <  5:415859/149‑1 (MQ=255)
                                        taGCAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCa                                                                                                          >  1:397134/1‑149 (MQ=255)
                                           cAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTg                                                                                                       >  8:370979/1‑149 (MQ=255)
                                                  ca‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTa                                                                                                >  6:332905/1‑149 (MQ=255)
                                                   a‑‑‑aGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTAc                                                                                               >  8:145238/1‑149 (MQ=255)
                                                        gagGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACgg                                                                                             <  2:392786/149‑1 (MQ=255)
                                                          ggCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTc                                                                                                     >  6:207821/1‑139 (MQ=255)
                                                          ggCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTc                                                                                                     <  5:207821/139‑1 (MQ=255)
                                                                    aTGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCtt                                                                                 >  5:149310/1‑149 (MQ=255)
                                                                       cgcgACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGc                                                                              >  3:221877/1‑149 (MQ=255)
                                                                        gcgACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGcc                                                                             >  2:16376/1‑149 (MQ=255)
                                                                           aCTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCggg                                                                          <  8:251765/149‑1 (MQ=255)
                                                                            cTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCggg                                                                          >  6:25124/1‑148 (MQ=255)
                                                                               ccaccaCAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGgga                                                                        <  1:290892/147‑1 (MQ=255)
                                                                               ccaccaCAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGgga                                                                        <  5:332905/147‑1 (MQ=255)
                                                                               ccaccaCAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGgga                                                                        <  5:25124/147‑1 (MQ=255)
                                                                                caccacAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGggagga                                                                     <  5:7370/149‑1 (MQ=255)
                                                                                caccacAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGggag                                                                       >  4:390659/1‑147 (MQ=255)
                                                                                      accacCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTgg                                                               <  3:321408/147‑1 (MQ=255)
                                                                                      aacaCCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTgg                                                               <  4:221877/149‑1 (MQ=255)
                                                                                      aacaCCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTgg                                                               <  6:149310/149‑1 (MQ=255)
                                                                                       acacCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTggg                                                              >  2:250698/1‑149 (MQ=255)
                                                                                         acCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGTAGGAGCCTgggaa                                                            >  5:244250/1‑147 (MQ=255)
                                                                                         acCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTCGGCa                                                            >  6:173058/1‑149 (MQ=255)
                                                                                         acCGGTATTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCa                                                            >  6:410281/1‑149 (MQ=255)
                                                                                          ccGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCtt                                                                                 >  8:275617/1‑127 (MQ=255)
                                                                                          ccGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCtt                                                                                 <  7:275617/127‑1 (MQ=255)
                                                                                          ccGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCAt                                                           >  7:7351/1‑149 (MQ=255)
                                                                                                  tCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACa                                                    <  2:65146/148‑1 (MQ=255)
                                                                                                    cACGCCGGCCATTCATCCTGGCTAGGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACa                                                    <  8:313755/146‑1 (MQ=255)
                                                                                                      cgcgGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTggg                                                              <  4:58295/134‑1 (MQ=255)
                                                                                                      cgcgGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTggg                                                              >  3:58295/1‑134 (MQ=255)
                                                                                                         gggCCATTCATGATGTCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCAtt                                            <  4:16296/149‑1 (MQ=255)
                                                                                                          ggCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCAttt                                           >  3:326949/1‑149 (MQ=255)
                                                                                                          ggCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCAttt                                           >  7:219124/1‑149 (MQ=255)
                                                                                                           gCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTc                                          >  2:321400/1‑149 (MQ=255)
                                                                                                            ccATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTcc                                         <  6:400774/149‑1 (MQ=255)
                                                                                                              aTTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCaa                                       >  5:2930/1‑149 (MQ=255)
                                                                                                                 catcatGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAAccc                                    >  6:31184/1‑149 (MQ=255)
                                                                                                                        gCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTg                                                                <  4:278043/114‑1 (MQ=255)
                                                                                                                        gCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTg                                                                >  3:278043/1‑114 (MQ=255)
                                                                                                                            tGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGc                           <  3:390659/147‑1 (MQ=255)
                                                                                                                             gACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGctgcct                                                                                     <  7:219482/88‑1 (MQ=255)
                                                                                                                             gACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGctgcct                                                                                     >  8:219482/1‑88 (MQ=255)
                                                                                                                              aCCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTg                       >  1:156351/1‑149 (MQ=255)
                                                                                                                               cccGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGc                                                  >  4:138466/1‑121 (MQ=255)
                                                                                                                               cccGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGc                                                  <  3:138466/121‑1 (MQ=255)
                                                                                                                                 cGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTg                    <  7:370979/149‑1 (MQ=255)
                                                                                                                                  gCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTgg                   >  3:24474/1‑149 (MQ=255)
                                                                                                                                   cTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGa                  <  1:202802/149‑1 (MQ=255)
                                                                                                                                   cTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGa                  <  7:145238/149‑1 (MQ=255)
                                                                                                                                        ccATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGagcagc             >  5:404987/1‑149 (MQ=255)
                                                                                                                                        ccATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGagcagc             <  2:156351/149‑1 (MQ=255)
                                                                                                                                        ccATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGagca               <  5:366930/147‑1 (MQ=255)
                                                                                                                                        ccATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGagca               >  6:366930/1‑147 (MQ=255)
                                                                                                                                        caccggtggCACCAGGGCGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGagcagc             >  6:156441/5‑149 (MQ=255)
                                                                                                                                         cATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGa                  <  7:331442/143‑1 (MQ=255)
                                                                                                                                         cATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGa                  >  8:331442/1‑143 (MQ=255)
                                                                                                                                         accggtggCACCAGGGCGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCCCCCTGGCCTTGCTGGACCcgcg            >  8:60616/4‑149 (MQ=255)
                                                                                                                                             gtggCACCAGGGCGCCACCGGTGTACAAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATg        >  6:386017/1‑149 (MQ=255)
                                                                                                                                                   cccTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATGGCTTCg  <  8:7351/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TATTTTGGAAAAAATCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACA‑‑‑AGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATGGCTTCG  >  NC_002947/1495507‑1495799

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