Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,450,023 T→A 5.1% C33S (TGC→AGC)  nicA → nicotinate dehydrogenase subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,450,0230TA 5.1% 246.5 / 2.2 79C33S (TGC→AGC) nicAnicotinate dehydrogenase subunit A
Reads supporting (aligned to +/- strand):  ref base T (38/36);  new base A (2/2);  total (40/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 = 9.51e-01

GTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCA  >  NC_002947/4449877‑4450170
                                                                                                                                                  |                                                                                                                                                     
gTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGc                                                                                                                                                     <  1:124872/149‑1 (MQ=255)
gTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGc                                                                                                                                                     >  5:106389/1‑149 (MQ=255)
    tttGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTg                                                                                                                                                  >  8:154557/1‑148 (MQ=255)
     ttGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGAGCCTGaa                                                                                                                                                >  6:190064/1‑149 (MQ=255)
     ttGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGAGCCTGaa                                                                                                                                                >  6:184689/1‑149 (MQ=255)
        aGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACgg                                                                                                                                             >  2:197191/1‑149 (MQ=255)
         gTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGc                                                                                                                                            <  6:106389/149‑1 (MQ=255)
          tCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGcc                                                                                                                                           <  3:172874/149‑1 (MQ=255)
            cTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGccc                                                                                                                                          >  5:349753/1‑148 (MQ=255)
             tttGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGcccc                                                                                                                                         >  2:318840/1‑148 (MQ=255)
               tGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCa                                                                                                                                        <  4:305735/147‑1 (MQ=255)
                ggTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACACCCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTCAACGTCCCAAAg                                                                                                                                      >  7:361627/1‑148 (MQ=255)
                 gTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTa                                                                                                                                    >  3:131906/1‑149 (MQ=255)
                    gcgcAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGc                                                                                                                                                     <  1:332181/129‑1 (MQ=255)
                    gcgcAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGc                                                                                                                                                     >  2:332181/1‑129 (MQ=255)
                     cgcAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGc                                                                                                                                >  6:141577/1‑149 (MQ=255)
                          gcaaCCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCgg                                                                                                                           >  3:281325/1‑149 (MQ=255)
                              ccTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTg                                                                                                                       <  3:126908/149‑1 (MQ=255)
                               cTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTgg                                                                                                                      >  4:165600/1‑149 (MQ=255)
                               cTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTgg                                                                                                                      >  1:276403/1‑149 (MQ=255)
                                tCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTggg                                                                                                                     >  1:173451/1‑149 (MQ=255)
                                 ccAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTgg                                                                                                                      >  6:184603/1‑147 (MQ=255)
                                 ccAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTgg                                                                                                                      <  5:184603/147‑1 (MQ=255)
                                  cAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGc                                                                                                                    <  6:349753/148‑1 (MQ=255)
                                     agaTGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAAt                                                                                                                <  7:154557/149‑1 (MQ=255)
                                          aGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCgg                                                                                                            <  3:353475/148‑1 (MQ=255)
                                          aGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCgg                                                                                                            <  3:165600/148‑1 (MQ=255)
                                          aGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGCGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCgg                                                                                                            <  3:180327/148‑1 (MQ=255)
                                           gCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCgggg                                                                                                          >  1:255077/1‑149 (MQ=255)
                                           gCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCgggg                                                                                                          <  2:276403/149‑1 (MQ=255)
                                            ccAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCg                                                                                                        >  5:361976/1‑149 (MQ=255)
                                                    gCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACgg                                                                                                 >  5:175547/1‑149 (MQ=255)
                                                      aaaCAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGt                                                                                                >  5:7669/1‑148 (MQ=255)
                                                      aaaCAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTg                                                                                               <  6:7669/149‑1 (MQ=255)
                                                      aaaCAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGAGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTg                                                                                               <  5:190064/149‑1 (MQ=255)
                                                      aaaCAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGAGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTg                                                                                               <  5:184689/149‑1 (MQ=255)
                                                       aacaacCATCTCACTGCAGGTCAACGGGCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGa                                                                                              <  8:361627/149‑1 (MQ=255)
                                                                   ccTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACg                                                                                    <  4:300851/146‑1 (MQ=255)
                                                                   ccTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACg                                                                                    <  2:173451/146‑1 (MQ=255)
                                                                     tGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTg                                                                               >  3:242414/1‑149 (MQ=255)
                                                                          gTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAAc                                                                                                                                               <  8:90689/81‑1 (MQ=255)
                                                                          gTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAAc                                                                                                                                               >  7:90689/1‑81 (MQ=255)
                                                                            cAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGccgccc                                                                        >  8:169553/1‑149 (MQ=255)
                                                                               cGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTg                                                                                               >  6:384922/1‑124 (MQ=255)
                                                                               cGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTg                                                                                               <  5:384922/124‑1 (MQ=255)
                                                                               cGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCt                                                                    >  8:76486/1‑149 (MQ=255)
                                                                                ggCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTc                                                                   <  4:281325/149‑1 (MQ=255)
                                                                                   cAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTg                                                                >  6:215761/1‑149 (MQ=255)
                                                                                   cAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTg                                                                <  6:175547/149‑1 (MQ=255)
                                                                                    aGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCg                                                              <  6:361976/149‑1 (MQ=255)
                                                                                     gCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCACTCGACGGCGTGGCCGCCCGCTCCTGc                                                               <  1:197191/149‑1 (MQ=255)
                                                                                      cccGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGt                                                             >  5:280403/1‑149 (MQ=255)
                                                                                          tCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATc                                                         >  7:322784/1‑149 (MQ=255)
                                                                                               gTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTTCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCtt                                                    >  5:95722/1‑149 (MQ=255)
                                                                                               gTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCtt                                                    <  5:141577/149‑1 (MQ=255)
                                                                                                  aGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGAGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTgg                                                  <  1:318840/148‑1 (MQ=255)
                                                                                                    cgcCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCg                                               <  1:199011/149‑1 (MQ=255)
                                                                                                     gcCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCgg                                              <  4:131906/149‑1 (MQ=255)
                                                                                                      ccATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGgc                                             >  7:363544/1‑149 (MQ=255)
                                                                                                           ccTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGc                                                               >  7:94857/1‑126 (MQ=255)
                                                                                                           ccTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGc                                                               <  8:94857/126‑1 (MQ=255)
                                                                                                             tGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGccgccg                                      >  2:83494/1‑149 (MQ=255)
                                                                                                             tGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGccgccg                                      >  4:142859/1‑149 (MQ=255)
                                                                                                              gACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCgccgg                                     >  3:48690/1‑149 (MQ=255)
                                                                                                                 acCCCGCTGCTGCTGATTCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATc                                                         <  4:94300/126‑1 (MQ=255)
                                                                                                                 acCCCGCTGCTGCTGATTCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATc                                                         >  3:94300/1‑126 (MQ=255)
                                                                                                                   cccGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATcctc                                                      >  1:169837/1‑125 (MQ=255)
                                                                                                                   cccGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATcctc                                                      <  2:169837/127‑3 (MQ=255)
                                                                                                                   cccGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGc                                 >  4:349741/1‑148 (MQ=255)
                                                                                                                   cccGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCa                                <  6:280403/149‑1 (MQ=255)
                                                                                                                     cGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTg                                                   <  4:100764/128‑1 (MQ=255)
                                                                                                                     cGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTg                                                   >  3:100764/1‑128 (MQ=255)
                                                                                                                     cGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAat                              <  8:322784/149‑1 (MQ=255)
                                                                                                                          ctgctgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATcac                         <  4:242414/149‑1 (MQ=255)
                                                                                                                             ctgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCAGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATcac                         <  1:83494/147‑1 (MQ=255)
                                                                                                                                      cgcAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGaa               <  7:76486/147‑1 (MQ=255)
                                                                                                                                        cAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTGATCCCCTTGGCCGGCGCCGGCGGCCGCACTCTCCCCACCCTCGAAGGGc           >  8:297569/1‑149 (MQ=255)
                                                                                                                                         aaCGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCt          >  3:88084/1‑149 (MQ=255)
                                                                                                                                           cGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTgg        <  4:88084/149‑1 (MQ=255)
                                                                                                                                                  tgCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgcagca  <  8:363544/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTCATTTGAGTCCTTTGGTAGCGCAAGCAACCTCCAGAGATGAGCCAATAATGCAAACAACCATCTCCCTGCAGGTCAACGGCCAGCCCGTCGAAGTCAGCGCCATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCA  >  NC_002947/4449877‑4450170

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