Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,950,559 (C)5→6 intergenic (‑37/‑273) aspS ← / → yecD aspartyl‑tRNA synthetase/isochorismatase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,950,5541.C98.5% 233.9 / ‑3.2 66intergenic (‑32/‑278)aspS/yecDaspartyl‑tRNA synthetase/isochorismatase family protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (34/31);  total (35/31)
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.98e-01

AGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑CCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTACC  >  NC_000913/1950425‑1950684
                                                                                                                                  |                                                                                                                                   
aGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑ccccc                                                                                                                                >  2:1212392/1‑135 (MQ=255)
aGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑ccccc                                                                                                                                <  1:1212392/135‑1 (MQ=255)
 ggCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑cccccctgct                                                                                                                           >  1:302330/1‑134 (MQ=255)
 ggCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑cccccctgct                                                                                                                           >  1:171227/1‑134 (MQ=255)
       cAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAg                                                                                                                     >  2:10147/1‑139 (MQ=255)
       cAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAg                                                                                                                     <  2:302330/139‑1 (MQ=255)
            tCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGAc                                                                                                                >  2:532320/1‑139 (MQ=255)
            tCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGAc                                                                                                                >  2:623255/1‑139 (MQ=255)
             cACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGAcc                                                                                                               >  1:915193/1‑139 (MQ=255)
              acgacgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCa                                                                                                              <  1:490796/139‑1 (MQ=255)
                gacgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAg                                                                                                            >  1:27372/1‑139 (MQ=255)
                  cgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTc                                                                                                          <  1:623255/139‑1 (MQ=255)
                       gTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGc                                                                                                     >  1:291008/1‑139 (MQ=255)
                         tGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGt                                                                                                   <  2:27372/139‑1 (MQ=255)
                            cccAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAg                                                                                                <  1:924617/139‑1 (MQ=255)
                            cccAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAg                                                                                                <  2:291008/139‑1 (MQ=255)
                               aaCCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                             <  1:532320/139‑1 (MQ=255)
                                     acaGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑cccccctg                                                                                                                             <  1:254756/101‑4 (MQ=255)
                                     acaGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑cccccctg                                                                                                                             >  2:254756/1‑98 (MQ=255)
                                         agTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACaa                                                                                   >  1:72510/1‑139 (MQ=255)
                                          gTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑ccc                                                                                                                                  <  1:1259511/91‑1 (MQ=255)
                                          gTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑ccc                                                                                                                                  >  2:1259511/1‑91 (MQ=255)
                                           tCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACt                                                                                 >  2:658340/1‑139 (MQ=255)
                                                 gctgcCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACtgaatga                                                                           >  2:572573/1‑139 (MQ=255)
                                                     ccccACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGaaaaa                                                                       >  1:1018390/1‑139 (MQ=255)
                                                        cACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAgg                                                                    >  1:437708/1‑139 (MQ=255)
                                                         aCGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                   >  2:835113/1‑139 (MQ=255)
                                                                aCAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCtatt                                                            <  2:1018390/139‑1 (MQ=255)
                                                                     cGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACt                                                       >  2:1191193/1‑139 (MQ=255)
                                                                            gTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAAtt                                                <  1:10147/139‑1 (MQ=255)
                                                                             tCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                             <  1:741521/93‑1 (MQ=255)
                                                                             tCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                             >  2:741521/1‑93 (MQ=255)
                                                                                acaATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACtgaa                                                                              >  1:874486/1‑105 (MQ=255)
                                                                                acaATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACtgaa                                                                              <  2:874486/105‑1 (MQ=255)
                                                                                acaATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGc                                            <  1:658340/139‑1 (MQ=255)
                                                                                acaATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGc                                            <  2:915193/139‑1 (MQ=255)
                                                                                       tCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCg                                     >  2:1412382/1‑139 (MQ=255)
                                                                                       tCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCg                                     >  2:891309/1‑139 (MQ=255)
                                                                                              gCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                   <  1:728382/102‑1 (MQ=255)
                                                                                              gCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                   >  2:728382/1‑102 (MQ=255)
                                                                                                aTGAGATATCCCTTAACTTAGCTGTCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTa                                   <  1:132182/132‑1 (MQ=255)
                                                                                                aTGAGATATCCCTTAACTTAGCTGTCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTa                                   >  2:132182/1‑132 (MQ=255)
                                                                                                aTGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGt                                                                                            <  1:148073/75‑1 (MQ=255)
                                                                                                aTGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGt                                                                                            >  2:148073/1‑75 (MQ=255)
                                                                                                aTGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                   <  1:398151/100‑1 (MQ=255)
                                                                                                aTGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                   >  2:398151/1‑100 (MQ=255)
                                                                                                aTGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCt                            >  1:1096030/1‑139 (MQ=255)
                                                                                                aTGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑cccccc                                                                                                                               <  1:1329486/40‑2 (MQ=255)
                                                                                                aTGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑cccccc                                                                                                                               >  2:1329486/1‑39 (MQ=255)
                                                                                                aTGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑CCCCCGTGCTTCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAACAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCt                            >  1:998974/1‑139 (MQ=255)
                                                                                                    gaTATCCCTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGcc                        <  2:437708/139‑1 (MQ=255)
                                                                                                           cTTAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTg                 >  2:1263300/1‑139 (MQ=255)
                                                                                                            ttAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAAc                                                                                  >  2:1129303/1‑73 (MQ=255)
                                                                                                            ttAACTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAAc                                                                                  <  1:1129303/73‑1 (MQ=255)
                                                                                                              aaCTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                   >  1:1106781/1‑86 (MQ=255)
                                                                                                              aaCTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                   <  2:1106781/86‑1 (MQ=255)
                                                                                                               aCTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAAc             >  1:1017290/1‑139 (MQ=255)
                                                                                                                cTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                             <  2:1459740/58‑1 (MQ=39)
                                                                                                                cTTAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                             >  1:1459740/1‑58 (MQ=39)
                                                                                                                  tAGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGcc          >  1:1013821/1‑139 (MQ=255)
                                                                                                                   aGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAg                                                                                                <  1:1332601/52‑1 (MQ=38)
                                                                                                                   aGCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAg                                                                                                >  2:1332601/1‑52 (MQ=38)
                                                                                                                    gCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCt         >  2:129056/1‑138 (MQ=255)
                                                                                                                    gCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCt         <  1:129056/138‑1 (MQ=255)
                                                                                                                    gCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTc        <  1:1244312/139‑1 (MQ=255)
                                                                                                                    gCTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTc        >  1:364255/1‑139 (MQ=255)
                                                                                                                     cTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                             <  2:931549/53‑1 (MQ=38)
                                                                                                                     cTGCCGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                             >  1:931549/1‑53 (MQ=38)
                                                                                                                        ccGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                   >  1:713404/1‑76 (MQ=255)
                                                                                                                        ccGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                   <  2:713404/76‑1 (MQ=255)
                                                                                                                        ccGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTa    <  1:835113/139‑1 (MQ=255)
                                                                                                                        ccGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTa    <  2:1013821/139‑1 (MQ=255)
                                                                                                                        ccGGCGGATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTa    <  2:1017290/139‑1 (MQ=255)
                                                                                                                         cggcggATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTAc   <  2:364255/139‑1 (MQ=255)
                                                                                                                          ggcggATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGc                                                                                                     <  1:317416/40‑1 (MQ=37)
                                                                                                                          ggcggATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGc                                                                                                     >  2:317416/1‑40 (MQ=37)
                                                                                                                          ggcggATGCCCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGTGCCTGGCCCGCGCTGGAACGCCTCGTTAcc  <  2:1096030/139‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                   
AGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑CCCCC‑TGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTACC  >  NC_000913/1950425‑1950684

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