Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,499,500 G→A *302* (TAG→TAA insD1 → IS2 transposase TnpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,499,5000GA100.0% 265.0 / NA 87*302* (TAG→TAAinsD1IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (40/47);  total (40/47)

GCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAG  >  NC_000913/4499375‑4499624
                                                                                                                               |                                                                                                                            
gCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                 <  1:452487/139‑1 (MQ=11)
gCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                 <  2:1183603/139‑1 (MQ=11)
     ggCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                            <  2:963536/139‑1 (MQ=21)
          tCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGat                                                                                                       <  1:1026517/139‑1 (MQ=255)
          tCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGat                                                                                                       <  1:627629/139‑1 (MQ=255)
              gCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGa                                                                                                   >  2:1403807/1‑139 (MQ=255)
                    taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  1:1339031/139‑1 (MQ=255)
                    taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  1:1118388/139‑1 (MQ=255)
                         aaTGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATAtt                                                                                        >  2:1185438/1‑139 (MQ=255)
                          aTGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                 >  1:665961/1‑113 (MQ=11)
                          aTGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                 <  2:665961/113‑1 (MQ=11)
                          aTGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCtt                                                                                                               <  1:815916/115‑1 (MQ=17)
                          aTGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCtt                                                                                                               >  2:815916/1‑115 (MQ=17)
                              cATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTtata                                                                                   >  2:914599/1‑139 (MQ=255)
                                tCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATg                                                                                                    <  1:420105/120‑1 (MQ=255)
                                tCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATg                                                                                                    >  2:420105/1‑120 (MQ=255)
                                   gCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  1:91722/124‑1 (MQ=255)
                                   gCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             >  2:91722/1‑124 (MQ=255)
                                      tAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                 >  2:223282/1‑103 (MQ=11)
                                      tAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                 <  1:223282/103‑1 (MQ=11)
                                       aGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACt                                                                          >  1:750671/1‑139 (MQ=255)
                                         tGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACttt                                                                        >  1:137006/1‑139 (MQ=255)
                                           cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  2:434941/116‑1 (MQ=255)
                                           cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             >  1:434941/1‑116 (MQ=255)
                                            gcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCtt                                                                                                               >  1:574104/1‑97 (MQ=18)
                                            gcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCtt                                                                                                               <  2:574104/97‑1 (MQ=18)
                                               gggTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGat                                                                                                       <  2:293599/102‑1 (MQ=255)
                                               gggTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGat                                                                                                       >  1:293599/1‑102 (MQ=255)
                                                 gTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACa                                                                <  1:1429251/139‑1 (MQ=255)
                                                   tATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAg                                                              >  1:1387784/1‑139 (MQ=255)
                                                    aTCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                             <  2:137006/139‑1 (MQ=255)
                                                     tcgctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCg                                                            >  1:603734/1‑139 (MQ=255)
                                                      cgctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                           <  1:1185438/139‑1 (MQ=255)
                                                        ctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCtt                                                                                                               >  1:1230072/1‑85 (MQ=18)
                                                        ctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCtt                                                                                                               <  2:1230072/85‑1 (MQ=18)
                                                        ctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCtt                                                                                                               >  1:1055067/1‑85 (MQ=18)
                                                        ctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCtt                                                                                                               <  2:1055067/85‑1 (MQ=18)
                                                        ctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACt                                                         <  1:807943/139‑1 (MQ=255)
                                                          cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                             <  1:451355/133‑1 (MQ=255)
                                                          cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                             >  2:451355/1‑133 (MQ=255)
                                                                gggAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACa                                                 <  2:1387784/139‑1 (MQ=255)
                                                                   aaTATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTc                                              >  1:1240647/1‑139 (MQ=255)
                                                                    atatCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCt                                             <  2:388897/139‑1 (MQ=255)
                                                                     tatCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTg                                            <  2:750671/139‑1 (MQ=255)
                                                                       tCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                           >  2:989249/1‑122 (MQ=255)
                                                                       tCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                           <  1:989249/122‑1 (MQ=255)
                                                                           cGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCatat                                      <  2:603734/139‑1 (MQ=255)
                                                                            ggCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                            <  1:1230196/68‑1 (MQ=255)
                                                                            ggCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                            >  2:1230196/1‑68 (MQ=255)
                                                                             gCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACt                                                         <  1:633835/118‑1 (MQ=255)
                                                                             gCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACt                                                         >  2:633835/1‑118 (MQ=255)
                                                                                 cGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGat                                                                                                       >  2:1009609/1‑68 (MQ=255)
                                                                                 cGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGat                                                                                                       <  1:1009609/68‑1 (MQ=255)
                                                                                 cGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGAt                                 >  1:60803/1‑138 (MQ=255)
                                                                                 cGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGAt                                 <  2:60803/138‑1 (MQ=255)
                                                                                  gggCTTGTATTGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTa                               <  1:914599/139‑1 (MQ=255)
                                                                                   ggCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                 <  2:1102832/56‑1 (MQ=9)
                                                                                   ggCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                 >  1:1102832/1‑56 (MQ=9)
                                                                                   ggCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTAc                              <  1:460148/139‑1 (MQ=255)
                                                                                             gggTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAg                    <  1:1403807/139‑1 (MQ=255)
                                                                                              ggTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCat                                        >  2:637446/1‑118 (MQ=255)
                                                                                              ggTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCat                                        <  1:637446/118‑1 (MQ=255)
                                                                                              ggTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGAt                       >  1:1331680/1‑135 (MQ=255)
                                                                                              ggTTAA‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGAt                       <  2:1331680/135‑1 (MQ=255)
                                                                                                   a‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTa                                                   <  2:327908/102‑1 (MQ=255)
                                                                                                   a‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTa                                                   >  1:327908/1‑102 (MQ=255)
                                                                                                   a‑‑GTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGgaagaa              >  1:3087/1‑139 (MQ=255)
                                                                                                          tAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                           >  2:731602/1‑89 (MQ=255)
                                                                                                          tAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                           <  1:731602/89‑1 (MQ=255)
                                                                                                          tAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAg                                                <  2:1337067/100‑1 (MQ=255)
                                                                                                          tAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAg                                                >  1:1337067/1‑100 (MQ=255)
                                                                                                          tAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGt                                               <  2:725790/101‑1 (MQ=255)
                                                                                                          tAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGt                                               >  1:725790/1‑101 (MQ=255)
                                                                                                                aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             >  1:799677/1‑49 (MQ=255)
                                                                                                                aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  2:799677/49‑1 (MQ=255)
                                                                                                                aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGcc                                          <  1:11498/100‑1 (MQ=255)
                                                                                                                aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGcc                                          >  2:11498/1‑100 (MQ=255)
                                                                                                                aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCtgtg                         >  2:294411/1‑117 (MQ=255)
                                                                                                                aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCtgtg                         <  1:294411/117‑1 (MQ=255)
                                                                                                                 tGTCTGGAAGTATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCATGAAGAATAAGTGGCTGAg  >  1:960736/1‑139 (MQ=255)
                                                                                                                 tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             >  1:481014/1‑48 (MQ=255)
                                                                                                                 tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  2:481014/48‑1 (MQ=255)
                                                                                                                 tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                           <  2:1284353/82‑1 (MQ=255)
                                                                                                                 tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                           >  1:1284353/1‑82 (MQ=255)
                                                                                                                 tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAg  >  1:867189/1‑139 (MQ=255)
                                                                                                                   tCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTtata                                                                                   >  1:990900/1‑56 (MQ=255)
                                                                                                                   tCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTtata                                                                                   <  2:990900/56‑1 (MQ=255)
                                                                                                                               |                                                                                                                            
GCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAA‑‑GTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAG  >  NC_000913/4499375‑4499624

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