Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 642,964 A→G H112R (CAC→CGC) 
V76A (GTG→GCG) 
b0609 →
rnk ←
hypothetical protein
regulator of nucleoside diphosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913642,9640AG88.7% 751.8 / 76.8 293H112R (CAC→CGC) 
V76A (GTG→GCG) 
b0609
rnk
hypothetical protein
regulator of nucleoside diphosphate kinase
Reads supporting (aligned to +/- strand):  ref base A (23/10);  new base G (134/126);  total (157/136)
Fisher's exact test for biased strand distribution p-value = 6.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

Only 100 of 297 total aligned reads displayed.

TGGAATCGCCAACGCGCAGCCCCAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACACCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCATCTCTTCTG  >  NC_000913/642872‑643054
                                                                                            |                                                                                          
  gAATCGCCAACGCGCAGCCCCAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGcc                                                                                          <  1:4190566/93‑1 (MQ=255)
    aTCGCCAACGCGCACCCCCAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAg                                                                                        <  1:3725906/93‑1 (MQ=255)
     tCGCCAACGCGCAGCCCCAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGc                                                                                       <  1:3104461/93‑1 (MQ=255)
     tCGCCAACGCGCAGCCCCAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGc                                                                                       >  1:4047238/1‑93 (MQ=255)
       gCCAACGCGCAGCCCCAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGt                                                                                     <  1:6749681/93‑1 (MQ=255)
        ccAACGCGCAGCCCCAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTg                                                                                    <  2:4269259/93‑1 (MQ=255)
          aaCGCGCAGCCCCAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACACCAGCGTgcg                                                                                  >  2:4623356/1‑93 (MQ=255)
           aCGCGCAGCCCCAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTgcgc                                                                                 >  2:5953234/1‑93 (MQ=255)
             gcgcAGCCCCAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGcac                                                                               >  1:4915518/1‑93 (MQ=255)
             gcgcAGCCCCAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGcac                                                                               >  1:5507852/1‑93 (MQ=255)
             gcgcAGCCCCAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGcac                                                                               >  1:3755764/1‑93 (MQ=255)
             gcgcAGCCCCAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACACCAGCGTGCGcac                                                                               >  2:3206516/1‑93 (MQ=255)
              cgcAGCCCCAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACACCAGCGTGCGcaca                                                                              <  2:1466518/93‑1 (MQ=255)
               gcAGCCCCAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCacac                                                                             <  2:2356400/93‑1 (MQ=255)
                cAGCCCCAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACg                                                                            <  1:3879460/93‑1 (MQ=255)
                  gCCCCAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCa                                                                          >  2:2592306/1‑93 (MQ=255)
                  gCCCCAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACACCAGCGTGCGCACACGCa                                                                          >  2:6808358/1‑93 (MQ=255)
                   ccccAGCAGTGCGGCGCCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCAc                                                                         <  1:3051354/93‑1 (MQ=255)
                   ccccAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCAc                                                                         >  1:4353978/1‑93 (MQ=255)
                     ccAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACtt                                                                       >  1:5181385/1‑93 (MQ=255)
                     ccAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACtt                                                                       >  2:7263485/1‑93 (MQ=255)
                     ccAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACtt                                                                       >  2:3098454/1‑93 (MQ=255)
                     ccAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACtt                                                                       >  1:2448935/1‑93 (MQ=255)
                      cagcagTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTc                                                                      >  2:4437341/1‑93 (MQ=255)
                      cagcagTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTc                                                                      >  2:2526280/1‑93 (MQ=255)
                      cagcagTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTc                                                                      >  1:2270796/1‑93 (MQ=255)
                      cagcagTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTc                                                                      <  1:5931825/93‑1 (MQ=255)
                      cagcagTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTc                                                                      <  2:6260764/93‑1 (MQ=255)
                      cagcagTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTc                                                                      >  2:2032689/1‑93 (MQ=255)
                       agcagTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCg                                                                     >  1:6384734/1‑93 (MQ=255)
                        gcagTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGc                                                                    >  2:5649103/1‑93 (MQ=255)
                        gcagTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGc                                                                    >  2:1145459/1‑93 (MQ=255)
                            tGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATc                                                                >  2:2214815/1‑93 (MQ=255)
                            tGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATc                                                                >  1:1392606/1‑93 (MQ=255)
                            tGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATc                                                                <  2:7054039/93‑1 (MQ=255)
                              cGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGc                                                              <  2:5555135/93‑1 (MQ=255)
                              cGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGc                                                              <  1:5386311/93‑1 (MQ=255)
                              cGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGc                                                              <  2:5414368/93‑1 (MQ=255)
                                gCACCTACCGGAGCCATAACCGAAAGCTGTCCATTGCTAACGGTCATTTTTGCCGGATCCGCCACCTTGCGGACACGCTCTTCGCCATAGCTa                                                            >  1:7187304/1‑93 (MQ=255)
                                 cACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTaa                                                           <  2:5336912/93‑1 (MQ=255)
                                  aCCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAg                                                          >  2:942710/1‑93 (MQ=255)
                                   ccTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGa                                                         <  1:1184508/93‑1 (MQ=255)
                                   ccTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGa                                                         >  2:1773349/1‑93 (MQ=255)
                                     tACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGAtt                                                       <  1:6089546/93‑1 (MQ=255)
                                     tACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGAtt                                                       <  2:3337442/93‑1 (MQ=255)
                                        cGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCg                                                    <  1:1571050/93‑1 (MQ=255)
                                        cGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCg                                                    >  2:5714062/1‑93 (MQ=255)
                                          gAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGa                                                  <  1:5675678/93‑1 (MQ=255)
                                            gCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAAt                                                <  1:5768171/93‑1 (MQ=255)
                                            gCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAAt                                                <  2:7417507/93‑1 (MQ=255)
                                            gCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAAt                                                >  1:558997/1‑93 (MQ=255)
                                            gCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACACCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAAt                                                <  2:4052129/93‑1 (MQ=255)
                                             ccATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAAtt                                               >  2:166614/1‑93 (MQ=255)
                                              cATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAAttt                                              >  1:295807/1‑93 (MQ=255)
                                              cATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAAttt                                              <  1:4740767/93‑1 (MQ=255)
                                              cATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAAttt                                              >  1:4593556/1‑93 (MQ=255)
                                               aTAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTa                                             >  2:659964/1‑93 (MQ=255)
                                                tAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTaa                                            <  2:1337607/93‑1 (MQ=255)
                                                tAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTaa                                            <  1:3518362/93‑1 (MQ=255)
                                                   cGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAAccc                                         <  1:1060206/93‑1 (MQ=255)
                                                    ggAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCg                                        >  2:4750803/1‑93 (MQ=255)
                                                     gAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACACCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCgg                                       >  2:2934576/1‑93 (MQ=255)
                                                       aaGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCt                                     >  2:7150504/1‑93 (MQ=255)
                                                       aaGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACACCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCt                                     <  1:3387059/93‑1 (MQ=255)
                                                         gCTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGt                                   >  1:6997125/1‑93 (MQ=255)
                                                          cTGAGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGtt                                  <  1:7446063/93‑1 (MQ=255)
                                                             aGTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCAt                               <  1:752739/93‑1 (MQ=255)
                                                              gTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCAtt                              >  2:6181986/1‑93 (MQ=255)
                                                              gTATTGCTATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCAtt                              <  2:834985/93‑1 (MQ=255)
                                                                     tATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTcacca                       >  2:6702584/1‑93 (MQ=255)
                                                                     tATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTcacca                       >  1:280803/1‑93 (MQ=255)
                                                                     tATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTcacca                       <  2:4302352/93‑1 (MQ=255)
                                                                     tATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTcacca                       <  1:934490/93‑1 (MQ=255)
                                                                     tATCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTcacca                       >  1:4017367/1‑93 (MQ=255)
                                                                       tCGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCAcg                     >  1:1885642/1‑93 (MQ=255)
                                                                        cGGTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCAcgt                    <  1:1765222/93‑1 (MQ=255)
                                                                         ggTCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCAcgtc                   <  2:6888052/93‑1 (MQ=255)
                                                                           tCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCAcgtcgt                 >  2:7320547/1‑93 (MQ=255)
                                                                           tCATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCAcgtcgt                 >  2:6831564/1‑93 (MQ=255)
                                                                            cATTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCgtg                >  2:3712644/1‑93 (MQ=255)
                                                                             aTTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCgtgt               >  1:3189893/1‑93 (MQ=255)
                                                                             aTTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCgtgt               <  2:6458395/93‑1 (MQ=255)
                                                                             aTTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCgtgt               <  1:1645983/93‑1 (MQ=255)
                                                                             aTTTTTGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCgtgt               >  2:3835944/1‑93 (MQ=255)
                                                                              tttttGCCGGATACGCCAGCGTGCGCACACGCACTTCGCGATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGtgtg              >  1:6211771/1‑93 (MQ=255)
                                                                              tttttGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGtgtg              <  1:2594586/93‑1 (MQ=255)
                                                                                 ttGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCa           <  1:1299897/93‑1 (MQ=255)
                                                                                 ttGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCa           <  1:6273166/93‑1 (MQ=255)
                                                                                  tGCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCAt          >  2:2264489/1‑93 (MQ=255)
                                                                                   gCCGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCAtc         <  1:6504897/93‑1 (MQ=255)
                                                                                    ccGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCAtct        <  2:5361239/93‑1 (MQ=255)
                                                                                     cGGATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCAtctc       >  1:270729/1‑93 (MQ=255)
                                                                                       gATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCATCtctt     >  1:180501/1‑93 (MQ=255)
                                                                                       gATACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCATCtctt     >  1:2053676/1‑93 (MQ=255)
                                                                                        aTACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCATCtcttc    <  2:5172521/93‑1 (MQ=255)
                                                                                        aTACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCATCtcttc    <  2:512172/93‑1 (MQ=255)
                                                                                        aTACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCATCtcttc    >  1:262848/1‑93 (MQ=255)
                                                                                         tACGCCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCATCtcttct   <  2:3468077/93‑1 (MQ=255)
                                                                                          acacCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCATCTCTTCTg  >  2:5989293/1‑93 (MQ=255)
                                                                                          acacCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCATCTCTTCTg  <  2:3662103/93‑1 (MQ=255)
                                                                                            |                                                                                          
TGGAATCGCCAACGCGCAGCCCCAGCAGTGCGGCACCTACCGGAGCCATAACGGAAAGCTGAGTATTGCTATCGGTCATTTTTGCCGGATACACCAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCATCTCTTCTG  >  NC_000913/642872‑643054

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