Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,661,821 C→G C104S (TGC→TCC)  iscR ← DNA‑binding transcriptional dual regulator IscR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,661,8210CG96.8% 186.1 / NA 64C104S (TGC→TCC) iscRDNA‑binding transcriptional dual regulator IscR
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base G (24/37);  total (25/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTAC  >  NC_000913/2661681‑2661960
                                                                                                                                            |                                                                                                                                           
gggTGCGGGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGAt                                                                                                                                           <  1:384996/143‑1 (MQ=255)
        ggCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATcgcc                                                                                                                                   >  1:660104/1‑143 (MQ=255)
              acgggAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAAACGGAGAGACGGTCGCTCAAATCACGCCACAGAGCGTGGGCGAGAGATATATCGACGCCCtc                                                                                                                             >  2:667940/5‑142 (MQ=255)
                 tGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGCGTCAGGCCTCTATCACCGCCCCGGcc                                                                                                                          >  2:378305/1‑142 (MQ=255)
                    gTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACCGCGCGTCAGTCAGGGATTTATCGCCGCCCTGGCAgcc                                                                                                                       <  2:818386/143‑1 (MQ=255)
                     tATGCTGACGACCAAACACATCCAGCACTTCATGGTTATTAACCAGTTCGCCTAAAGTAATGTTATTGAGAAAACAGATGAGACGGTCGCTTAAACTACACCACAGCGCGAAGGTCAGCGATATATCGCAGACCAGGAAgccg                                                                                                                      >  2:2713/1‑143 (MQ=255)
                       tGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAgccgcc                                                                                                                    >  1:549297/1‑143 (MQ=255)
                         cTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACAGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGACTTATCGCCGCCCTGGCAgccgccgt                                                                                                                  >  2:1000554/1‑141 (MQ=255)
                           gacgacCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTa                                                                                                                <  2:695584/143‑1 (MQ=255)
                           gacgacCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTa                                                                                                                >  1:949465/1‑143 (MQ=255)
                            acgtctaGACACATCCAGCACTTCCTGGGTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTAc                                                                                                               <  1:983435/137‑1 (MQ=255)
                             cgacCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTAcc                                                                                                              <  2:665721/143‑1 (MQ=255)
                                 cAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTCACGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGGATTTATCGCCGCGCTGGCAGCCGCCTTTCCCCTGa                                                                                                          >  2:30342/1‑143 (MQ=255)
                                  atacacGTCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCGAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGAc                                                                                                         <  1:424379/141‑1 (MQ=255)
                                  aGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTg                                                                                                                             >  1:349764/1‑123 (MQ=255)
                                  aGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTg                                                                                                                             <  2:349764/123‑1 (MQ=255)
                                   gACACATCCATGACTTCCTTGTTATGAACCAGGTCTCTGAAAGTAATGTTGTTGAGAGAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACa                                                                                                        <  2:804150/143‑1 (MQ=255)
                                     cacaTCCAGCACTTCCTGGTTATGAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCACCTTTACCCTGACAAc                                                                                                      <  2:265019/143‑1 (MQ=255)
                                      acaTCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACg                                                                                                     <  1:547734/143‑1 (MQ=255)
                                         tCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGt                                                                                                  <  1:387926/143‑1 (MQ=255)
                                          ccAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCGCTGGCAGCCGCCTTTACCCTGACAACAGGTg                                                                                                 >  1:478618/1‑143 (MQ=255)
                                          ccAGCACTTCCTGGTCATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTg                                                                                                 <  1:595863/143‑1 (MQ=255)
                                           cAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTgg                                                                                                >  1:391638/1‑143 (MQ=255)
                                              cACTTCCTGGTTGTTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCAt                                                                                             <  2:553846/143‑1 (MQ=255)
                                              cACGTCCTGGTTATTAACCAGTTCGCCTAAAGTGATGGTGTTGAGAAAGCCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCAt                                                                                             <  2:58500/143‑1 (MQ=255)
                                                   ccTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACa                                                                                        >  1:154201/1‑143 (MQ=255)
                                                     tGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTCGCATCTACAGa                                                                                      >  2:850232/1‑143 (MQ=255)
                                                          attaACAAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGTTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGt                                                                                 <  1:408245/143‑1 (MQ=255)
                                                            taACCAGTTCGCCTAAAGTAATGTTGTTGAGAAGACCGGTGGGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGGGGCATCTACAGATTCGTCa                                                                               >  1:737334/1‑143 (MQ=255)
                                                            taACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACCACGGGTGGCATCTCCAGACTCGTcc                                                                               >  1:173405/1‑142 (MQ=255)
                                                                cAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACgg                                                                           <  2:549297/143‑1 (MQ=255)
                                                                 aGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCATCGgc                                                                          >  1:472720/1‑143 (MQ=255)
                                                                 aGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGgc                                                                          <  2:949465/143‑1 (MQ=255)
                                                                      gCCTAAAGTAATGTTGGTGGGAAAAACGGTGAGACGGTTGCTCCCAGCACGCCACAGCACGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTaa                                                                     <  2:860846/143‑1 (MQ=255)
                                                                       ccTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAGATCACGCCACAGCGCGTCGGTCCGGGGTCTATCGCCGCCCTGGCAGCCGCCCGTACCCTGACAACGGTTGGCCTCTACCGAATCCTCAACGgcacgaaa                                                                    >  2:341518/1‑137 (MQ=255)
                                                                       ccTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAAt                                                                    <  2:306477/143‑1 (MQ=255)
                                                                         tAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAAttc                                                                  >  1:981325/1‑142 (MQ=255)
                                                                           aaGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGACAGGGATTTATCGCCGCCCTGGCAGCCGCCTGTACCCTGACAACGGGTGGCATCTACATAGTCGTCAACGCCGCTAATTACt                                                                >  2:735487/1‑143 (MQ=255)
                                                                              tAAGGTTGTCGATAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCg                                                             <  1:970594/143‑1 (MQ=255)
                                                                                     gttgAGAAAACCGGTGAGACGGTCGATTAAAAATACAGATAAGGCGCGGGTCAAGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGAGGCATCCACAGACTCGACAACGGCGCTAGTCACTCCGCCAACgg                                                      >  2:551629/1‑143 (MQ=255)
                                                                                      ttgAGAAAACCGGTGAGACGGTCGCTCAAATCACGCGACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCGACGGCGCTAATTACTCAGCCAACGGc                                                     >  2:255076/1‑143 (MQ=255)
                                                                                        gagaAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTGAACGGCGCTAATTACTTCGCCAGCGGCGa                                                   >  2:547924/1‑143 (MQ=255)
                                                                                            aaaCCTGTGAGAGGGTCGCTCAAACCACGCCACAGCGCGTGGGTCACGGATTTATCGCCGCCCTGGCAGCCGCCTTCACCCTGACAACGGGTGGCATCTACAGATTCGTCTACGGCGCTAATTACTTCGCCAACGGCGAtgct                                               <  2:526566/143‑1 (MQ=255)
                                                                                               cccgcgaggCGGTCGGTCACATCACGGCGCAGCGCGTGGGTCAGGGATTTATCGCGGCCCTGGCAGCCGGCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATgctgct                                            <  2:55158/138‑1 (MQ=255)
                                                                                               ccGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATgctgct                                            <  1:647269/143‑1 (MQ=255)
                                                                                               ccGGTGAGACGGTCGCTCAAACCACGCCACAGCGCGCGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTCCAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATgctgct                                            <  1:514125/143‑1 (MQ=255)
                                                                                                        cGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGTGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCttt                                   <  1:738126/143‑1 (MQ=255)
                                                                                                          ggcGCTCAAATCACGCCGCAGCGCGGGGGGCAGGGAGTGATCGCCGCCCTGGCAGCCGCCGTTACCCTGACAACGGGTGGCACCTACAGAGTCGTCAACGGCGCTAATGACTTCGCCAACGGCGATGCTGCTGGCATCTTTGc                                 <  1:337501/141‑1 (MQ=255)
                                                                                                          gTCGCTCATATCACGCCACGGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGc                                 <  2:105593/143‑1 (MQ=255)
                                                                                                          gTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGc                                 <  1:909971/143‑1 (MQ=255)
                                                                                                          gTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGc                                 <  2:774070/143‑1 (MQ=255)
                                                                                                          gTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGc                                 <  2:154201/143‑1 (MQ=255)
                                                                                                              cTCAAATCACGCCACGGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTGTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTaa                             <  1:325580/143‑1 (MQ=255)
                                                                                                                     cACGCCACAGCGCGTGGGTGAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATaa                      <  1:3947/143‑1 (MQ=255)
                                                                                                                     cACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTAGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTGACCGATaa                      >  2:501715/1‑143 (MQ=255)
                                                                                                                       cGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCCACGGCGACGCTGCTGGCATCTTTGCCTAACAGATAacc                    <  1:490338/143‑1 (MQ=255)
                                                                                                                         acacaGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAaccac                  <  1:255739/142‑1 (MQ=255)
                                                                                                                          cacaGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCTCTGCTGCTGGCATCTTTGCCTAACAGATAgccacc                 >  1:600821/1‑143 (MQ=255)
                                                                                                                             aGCGCGTGGGCCAGGTATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGcc              <  2:410210/143‑1 (MQ=255)
                                                                                                                               cgcgTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTg            >  1:475447/1‑143 (MQ=255)
                                                                                                                                gcgTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTgg           <  1:134144/143‑1 (MQ=255)
                                                                                                                                  tggggTCCGGGCTTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTCGCATCTACGGATTCGTCAACGGCGCTAATTACGTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTc         <  2:456314/141‑1 (MQ=255)
                                                                                                                                        cAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTa   <  1:667768/143‑1 (MQ=255)
                                                                                                                                         aGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTAc  <  1:667940/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                           
GGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTAC  >  NC_000913/2661681‑2661960

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