Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 983,581 T→A 7.9% intergenic (+39/‑50) iscX → / → ndk cysteine desulfurase activity regulator/nucleoside diphosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947983,5810TA 7.9% 259.5 / 9.1 88intergenic (+39/‑50)iscX/ndkcysteine desulfurase activity regulator/nucleoside diphosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (40/41);  new base A (2/5);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01

ACCTGCACCGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCGGCCT  >  NC_002947/983433‑983725
                                                                                                                                                       |                                                                                                                                                
aCCTGCACCGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGt                                                                                                                                                  >  5:218608/1‑149 (MQ=255)
  cTGCACCGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGta                                                                                                                                                 >  2:151748/1‑148 (MQ=255)
     cACCGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATaa                                                                                                                                              <  2:86291/148‑1 (MQ=255)
     cACCGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAAt                                                                                                                                             >  2:34490/1‑149 (MQ=255)
     cACCGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAAt                                                                                                                                             >  2:28397/1‑149 (MQ=255)
        cGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTcg                                                                                                                                          >  1:392828/1‑149 (MQ=255)
             ggTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCgg                                                                                                                                       <  3:238782/147‑1 (MQ=255)
             ggTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGAATTATACGCggg                                                                                                                                      <  5:209671/148‑1 (MQ=255)
              gTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGt                                                                                                                                     <  1:282971/149‑1 (MQ=255)
              actcTGGCAATGCCTGAGTTCAGCGACGATTCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGt                                                                                                                                     >  1:237347/4‑149 (MQ=255)
               tGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGAATTATACGCGGGttc                                                                                                                                   <  1:43327/149‑2 (MQ=255)
               tGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGAAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGtt                                                                                                                                    >  7:149670/1‑148 (MQ=255)
                 cTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTaa                                                                                                                                 <  7:159734/149‑1 (MQ=255)
                    gctttGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGAATTATACGCGGGttccgg                                                                                                                                <  1:33379/144‑5 (MQ=255)
                    gCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGAATTATACGCGGGttccgg                                                                                                                                <  8:37446/147‑5 (MQ=255)
                    gCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGAATTATACGCGGGttccgg                                                                                                                                <  1:81180/147‑5 (MQ=255)
                         gCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCg                                                                                                                          <  2:126004/148‑1 (MQ=255)
                           cAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCgg                                                                                                                                       >  5:261580/1‑133 (MQ=255)
                           cAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCgg                                                                                                                                       <  6:261580/133‑1 (MQ=255)
                           cAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTaa                                                                                                                       <  6:283832/149‑1 (MQ=255)
                               gTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCTGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAacac                                                                                                                    <  4:105034/148‑1 (MQ=255)
                                  cAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCa                                                                                                                 >  6:271416/1‑148 (MQ=255)
                                   aGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCata                                                                                                               >  2:357633/1‑149 (MQ=255)
                                     cgacgaTCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATAtt                                                                                                             >  7:63058/1‑149 (MQ=255)
                                        cgaTCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTg                                                                                                          <  1:303112/149‑1 (MQ=255)
                                          aTCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGa                                                                                                        <  6:122004/149‑1 (MQ=255)
                                          aTCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGa                                                                                                        <  3:346284/149‑1 (MQ=255)
                                          aTCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGa                                                                                                        <  6:226602/149‑1 (MQ=255)
                                            tCCTTCACTCGGCGGTGATAAAGTGCTCTTGGCAATCCAGGCCGCCTGTATGGAGTAAGCGGACTATGCTTCACTTGAAGGTTAGGCAATCCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGa                                                                                                        <  1:320142/149‑1 (MQ=255)
                                              cGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGt                                                                                                      <  8:278114/148‑1 (MQ=255)
                                              cGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGt                                                                                                      <  8:192891/148‑1 (MQ=255)
                                               gTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGttt                                                                                                    >  6:31103/1‑149 (MQ=255)
                                                tCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGtttt                                                                                                   >  1:268729/1‑149 (MQ=255)
                                                 cACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGa                                                                                                        >  7:179117/1‑143 (MQ=255)
                                                    gcggcgGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCAt                                                                                               <  6:218608/149‑1 (MQ=255)
                                                      ggcggTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATgg                                                                                             <  1:151748/149‑1 (MQ=255)
                                                       gcggTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGc                                                                                            <  1:357633/149‑1 (MQ=255)
                                                          gTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGt                                                                                         <  7:363240/149‑1 (MQ=255)
                                                               aaaGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTNATATTCTGGAGTTTTCCATGGCTGTTCAAc                                                                                    <  1:34490/149‑1 (MQ=255)
                                                               aaaGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAAc                                                                                    <  1:28397/149‑1 (MQ=255)
                                                                     cTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACt                                                                               <  2:392828/148‑1 (MQ=255)
                                                                      tCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACttt                                                                             >  5:372800/1‑149 (MQ=255)
                                                                      tCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACttt                                                                             >  1:182817/1‑149 (MQ=255)
                                                                      tCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACttt                                                                             >  6:70999/1‑149 (MQ=255)
                                                                      tCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACtt                                                                              >  8:191174/1‑148 (MQ=255)
                                                                        gAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCGCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAGCACTCATATTCTGGAGTTTTCCATGGCTGTTCAGCGTACTTtct                                                                           >  4:368843/1‑149 (MQ=255)
                                                                        gAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTtct                                                                           >  2:115440/1‑149 (MQ=255)
                                                                         aGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTtctc                                                                          <  7:246200/149‑1 (MQ=255)
                                                                         aGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTtctc                                                                          <  2:300302/149‑1 (MQ=255)
                                                                         aGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTtctc                                                                          <  2:184189/149‑1 (MQ=255)
                                                                            cAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCcat                                                                       <  5:307228/149‑1 (MQ=255)
                                                                                  aGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCatcatc                                                                   >  6:97162/1‑147 (MQ=255)
                                                                                  aGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCatcatc                                                                   >  2:74402/1‑147 (MQ=255)
                                                                                  aGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCaa                                                                 <  2:237347/149‑1 (MQ=255)
                                                                                  aGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCaa                                                                 <  5:31103/149‑1 (MQ=255)
                                                                                    gccgccTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTc                                                                                                                  >  7:87714/1‑98 (MQ=255)
                                                                                    gccgccTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTc                                                                                                                  <  8:87714/98‑1 (MQ=255)
                                                                                         cTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTg                                                                                          <  3:44814/117‑1 (MQ=255)
                                                                                         cTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTg                                                                                          >  4:44814/1‑117 (MQ=255)
                                                                                          tGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACg                                                         >  3:368258/1‑149 (MQ=255)
                                                                                           ggATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGc                                                        >  1:381745/1‑149 (MQ=255)
                                                                                           ggATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGc                                                        >  4:351403/1‑149 (MQ=255)
                                                                                             aTCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTg                                                      >  3:97396/1‑149 (MQ=255)
                                                                                                gacgaAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTAc                                                                                >  7:317247/1‑120 (MQ=255)
                                                                                                gacgaAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTAc                                                                                <  8:317247/120‑1 (MQ=255)
                                                                                                 acgaAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGc                                                  >  4:387516/1‑149 (MQ=255)
                                                                                                 acgaAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGc                                                  >  6:358161/1‑149 (MQ=255)
                                                                                                    aaGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCa                                                >  4:123160/1‑148 (MQ=255)
                                                                                                        cGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAAc                                           <  5:97162/149‑1 (MQ=255)
                                                                                                         gACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACg                                          >  5:232800/1‑149 (MQ=255)
                                                                                                         gACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACg                                          <  5:358161/149‑1 (MQ=255)
                                                                                                          aCTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGt                                         >  2:246989/1‑149 (MQ=255)
                                                                                                          aCTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGt                                         <  5:271416/149‑1 (MQ=255)
                                                                                                                 gTCACTTGGAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGc                                  <  8:149670/149‑1 (MQ=255)
                                                                                                                        gCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATca                           >  2:362200/1‑149 (MQ=255)
                                                                                                                          cAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATcac                          <  2:268729/149‑1 (MQ=255)
                                                                                                                            tgTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATcacca                        <  5:190879/148‑1 (MQ=255)
                                                                                                                            ggTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATcacca                        <  3:123160/149‑1 (MQ=255)
                                                                                                                                   cAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCtt                 >  1:319367/1‑149 (MQ=255)
                                                                                                                                   cAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCtt                 <  4:368258/149‑1 (MQ=255)
                                                                                                                                         ccccGGAACCCGCGTATAATTCGCAGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCATGAAGCTGTTGCCAAAAACGTCATAGGCAAGATAACCACTCGCTTCtaa             >  7:187989/1‑148 (MQ=255)
                                                                                                                                           ccGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCgaagaa          <  3:351403/149‑1 (MQ=255)
                                                                                                                                            ttaaaCCCGCGAATTATACGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGaagaag         >  5:254936/4‑149 (MQ=255)
                                                                                                                                             ggAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGc        <  6:232800/149‑1 (MQ=255)
                                                                                                                                              aaaCCCGCGAATTATACGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGcc       >  8:227558/2‑149 (MQ=255)
                                                                                                                                                  ccGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCGGcc   >  5:167999/1‑149 (MQ=255)
                                                                                                                                                   cgcgTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCGGCCt  >  5:5910/1‑149 (MQ=255)
                                                                                                                                                   agcgTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCGGCCt  >  6:118853/2‑149 (MQ=255)
                                                                                                                                                    gcgTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCGGCCt  >  3:135413/1‑148 (MQ=255)
                                                                                                                                                    gcgTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTACATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGTCAAGATCACCACTCGCTTCGAAGAAGCCGGCCt  >  4:107819/1‑148 (MQ=255)
                                                                                                                                                       |                                                                                                                                                
ACCTGCACCGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCGGCCT  >  NC_002947/983433‑983725

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: