Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,018,760 G→A *476* (TAG→TAA rmuC → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,018,7600GA98.9% 285.2 / NA 89*476* (TAG→TAArmuCDNA recombination protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (41/47);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAGCCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACTTTG  >  NC_000913/4018627‑4018888
                                                                                                                                     |                                                                                                                                
aaCGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAt                                                                                                                             <  1:1247115/139‑1 (MQ=255)
 aCGCGAGGTTAATCCTGAATTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAAACAGCAGTCGCGCTAACCCAtt                                                                                                                            <  1:375824/139‑1 (MQ=255)
 aCGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAtt                                                                                                                            <  1:1020994/139‑1 (MQ=255)
 aCGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAtt                                                                                                                            <  2:659016/139‑1 (MQ=255)
 aCGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAtt                                                                                                                            >  2:1366434/1‑139 (MQ=255)
    cgAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTggg                                                                                                                         <  2:1691405/139‑1 (MQ=255)
    cgAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTggg                                                                                                                         >  1:1679045/1‑139 (MQ=255)
     gagaTTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGa                                                                                                                        <  1:1594880/139‑1 (MQ=255)
      agaTTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGag                                                                                                                       >  2:1136025/1‑139 (MQ=255)
               cGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGc                                                                                                              <  2:1590827/139‑1 (MQ=255)
                 gATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCg                                                                                                            <  1:1136025/139‑1 (MQ=255)
                 gATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCg                                                                                                            <  1:1216480/139‑1 (MQ=255)
                      ggCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAg                                                                                                       <  1:1420615/139‑1 (MQ=255)
                                     gAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCg                                                                                        <  2:792784/139‑1 (MQ=255)
                                          aGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                   <  2:313062/139‑1 (MQ=255)
                                          aGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                   <  2:1230237/139‑1 (MQ=255)
                                                     tATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGTAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                   >  2:874129/1‑128 (MQ=255)
                                                     tATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGTAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                   <  1:874129/128‑1 (MQ=255)
                                                     tATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTac                                                                        >  1:1418732/1‑139 (MQ=255)
                                                      tgtataagagacaggtTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTaca                                                                       >  1:1319576/14‑139 (MQ=255)
                                                      aTCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAtt                                                                                                                            >  2:763330/1‑86 (MQ=255)
                                                      aTCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAtt                                                                                                                            <  1:763330/86‑1 (MQ=255)
                                                       tCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTacac                                                                      >  2:1121467/1‑139 (MQ=255)
                                                               ggtcggtTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAg                                                                                                       <  1:37769/98‑1 (MQ=255)
                                                               ggtcggtTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAg                                                                                                       >  2:37769/1‑98 (MQ=255)
                                                               ggtcggtTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGaa                                                              >  1:1372637/1‑139 (MQ=255)
                                                                 tcggtTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGTCCAATCTGTTACACTTCTGGAACa                                                            >  2:1002055/1‑139 (MQ=255)
                                                                              agcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                               <  2:1319576/139‑1 (MQ=255)
                                                                               gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTa                                                                                                        <  2:669353/81‑1 (MQ=255)
                                                                               gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTa                                                                                                        >  1:669353/1‑81 (MQ=255)
                                                                               gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAg                                                                                                       <  1:1460834/82‑1 (MQ=255)
                                                                               gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAg                                                                                                       >  2:1460834/1‑82 (MQ=255)
                                                                               gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATc                                                                                                  <  1:476185/87‑1 (MQ=255)
                                                                               gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATc                                                                                                  >  2:476185/1‑87 (MQ=255)
                                                                               gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTgg                                                                >  1:1054148/1‑121 (MQ=255)
                                                                               gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTgg                                                                <  2:1054148/121‑1 (MQ=255)
                                                                               gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                              <  2:1514929/139‑1 (MQ=255)
                                                                                       atgaAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTg                                      <  1:1121467/139‑1 (MQ=255)
                                                                                         gaAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                   <  2:1382201/92‑1 (MQ=255)
                                                                                         gaAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                   >  1:1382201/1‑92 (MQ=255)
                                                                                          aaGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTgaga                                   >  2:1082814/1‑139 (MQ=255)
                                                                                           aGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGAt                                  <  1:1082814/139‑1 (MQ=255)
                                                                                             cTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGAtgg                                >  1:429576/1‑139 (MQ=255)
                                                                                                aTCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAg                                                                                                       <  1:53549/65‑1 (MQ=255)
                                                                                                aTCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAg                                                                                                       >  2:53549/1‑65 (MQ=255)
                                                                                                 tCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACt                                                                     >  2:97605/1‑98 (MQ=255)
                                                                                                 tCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACt                                                                     <  1:97605/98‑1 (MQ=255)
                                                                                                  cAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGAt                           <  2:1679045/139‑1 (MQ=255)
                                                                                                     cgcgATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGt                                                                           >  1:1309749/1‑88 (MQ=255)
                                                                                                     cgcgATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGt                                                                           <  2:1309749/88‑1 (MQ=255)
                                                                                                      gcgATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGag                                                                                                                       >  1:585709/1‑43 (MQ=255)
                                                                                                      gcgATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGag                                                                                                                       <  2:585709/43‑1 (MQ=255)
                                                                                                      gcgATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                               <  2:363706/115‑1 (MQ=255)
                                                                                                      gcgATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                               >  1:363706/1‑115 (MQ=255)
                                                                                                      gcgATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGt                       >  2:552083/1‑139 (MQ=255)
                                                                                                         atgatgAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTcac                    >  1:818624/1‑139 (MQ=255)
                                                                                                                atatAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACa                                                            >  2:740305/1‑92 (MQ=255)
                                                                                                                atatAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACa                                                            <  1:740305/92‑1 (MQ=255)
                                                                                                                   taATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAacgac          <  2:429576/139‑1 (MQ=255)
                                                                                                                    aaTCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAAt                                                                               <  1:678674/69‑1 (MQ=255)
                                                                                                                    aaTCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAAt                                                                               >  2:678674/1‑69 (MQ=255)
                                                                                                                    aaTCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAacgacg         >  1:726180/1‑139 (MQ=255)
                                                                                                                      tCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGcc                                                                                                             >  1:283278/1‑37 (MQ=255)
                                                                                                                      tCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGcc                                                                                                             <  2:283278/37‑1 (MQ=255)
                                                                                                                      tCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGcc                                                                                                             >  2:1510628/1‑37 (MQ=255)
                                                                                                                      tCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGcc                                                                                                             <  1:1510628/37‑1 (MQ=255)
                                                                                                                      tCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGAtgg                                <  2:444922/114‑1 (MQ=255)
                                                                                                                      tCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGAtgg                                >  1:444922/1‑114 (MQ=255)
                                                                                                                      tCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCa       <  1:242012/139‑1 (MQ=255)
                                                                                                                         gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTa                                                                                                        >  1:1712313/1‑39 (MQ=255)
                                                                                                                         gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTa                                                                                                        <  2:1712313/39‑1 (MQ=255)
                                                                                                                         gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGaa                                                                                                     >  2:402036/1‑42 (MQ=255)
                                                                                                                         gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGaa                                                                                                     <  1:402036/42‑1 (MQ=255)
                                                                                                                         gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                               >  2:888134/1‑96 (MQ=255)
                                                                                                                         gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                               <  1:888134/96‑1 (MQ=255)
                                                                                                                         gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                              <  2:1586678/97‑1 (MQ=255)
                                                                                                                         gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                              >  1:1586678/1‑97 (MQ=255)
                                                                                                                         gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACtt    >  2:1612848/1‑139 (MQ=255)
                                                                                                                         gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACtt    >  1:626604/1‑139 (MQ=255)
                                                                                                                         gcagTCGCACTAGCCCATCGGGAGTAGTTAAGCCGGGTAGAGAACTAGGGCACCGATGCCCGACCTGTTACACTTCTGGTACAATTTTTTGAAGAGCAGCCATTGAGAAGGAGGAAAAGCCACAAGAAaagacgacact    >  1:57344/1‑134 (MQ=255)
                                                                                                                          cagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                   >  2:1061965/1‑59 (MQ=255)
                                                                                                                          cagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                   <  1:1061965/59‑1 (MQ=255)
                                                                                                                          cagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACttt   <  2:1372637/139‑1 (MQ=255)
                                                                                                                           agTCGCGCTAACCCATTGGGTGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCa                                                                                 <  2:1619933/60‑1 (MQ=255)
                                                                                                                           agTCGCGCTAACCCATTGGGTGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCa                                                                                 >  1:1619933/1‑60 (MQ=255)
                                                                                                                           agTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACTTTg  <  1:1002055/139‑1 (MQ=255)
                                                                                                                           agTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACTTTg  >  1:302968/1‑139 (MQ=255)
                                                                                                                              cgcgCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCa                                                                                 <  1:1451443/57‑1 (MQ=255)
                                                                                                                              cgcgCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCa                                                                                 >  2:1451443/1‑57 (MQ=255)
                                                                                                                                     |                                                                                                                                
AACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAGCCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACTTTG  >  NC_000913/4018627‑4018888

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