Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 983,587 T→A 8.2% intergenic (+45/‑44) 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,5870TA 8.2% 251.5 / 9.2 85intergenic (+45/‑44)iscX/ndkcysteine desulfurase activity regulator/nucleoside diphosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (38/40);  new base A (2/5);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

CGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGA‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCGGCCTGAAAATCGTT  >  NC_002947/983441‑983735
                                                                                                                                                     |                                                                                                                                                    
cGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGA‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTcg                                                                                                                                                    >  1:392828/1‑149 (MQ=255)
     ggTGCTGGCATTGCCAGAGTTCAGCGACGA‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCgg                                                                                                                                                 <  3:238782/147‑1 (MQ=255)
     ggTGCTGGCATTGCCAGAGTTCAGCGACGA‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGAATTATACGCggg                                                                                                                                                <  5:209671/148‑1 (MQ=255)
      gTGCTGGCATTGCCAGAGTTCAGCGACGA‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCCAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGt                                                                                                                                               <  1:282971/149‑1 (MQ=255)
      actcTGGCAATGCCTGAGTTCAGCGACGATTCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGt                                                                                                                                               >  1:237347/4‑149 (MQ=255)
       tGCTGGCATTGCCAGAGTTCAGCGACGA‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGAATTATACGCGGGttc                                                                                                                                             <  1:43327/149‑2 (MQ=255)
       tGCTGGCATTGCCAGAGTTCAGCGACGA‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑AAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGtt                                                                                                                                              >  7:149670/1‑148 (MQ=255)
         cTGGCATTGCCAGAGTTCAGCGACGA‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTaa                                                                                                                                           <  7:159734/149‑1 (MQ=255)
            gctttGCCAGAGTTCAGCGACGA‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGAATTATACGCGGGttccgg                                                                                                                                          <  1:33379/144‑5 (MQ=255)
            gCATTGCCAGAGTTCAGCGACGA‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGAATTATACGCGGGttccgg                                                                                                                                          <  1:81180/147‑5 (MQ=255)
            gCATTGCCAGAGTTCAGCGACGA‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGAATTATACGCGGGttccgg                                                                                                                                          <  8:37446/147‑5 (MQ=255)
                 gCCAGAGTTCAGCGACGA‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCg                                                                                                                                    <  2:126004/148‑1 (MQ=255)
                   cAGAGTTCAGCGACGA‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCgg                                                                                                                                                 >  5:261580/1‑133 (MQ=255)
                   cAGAGTTCAGCGACGA‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCgg                                                                                                                                                 <  6:261580/133‑1 (MQ=255)
                   cAGAGTTCAGCGACGA‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTaa                                                                                                                                 <  6:283832/149‑1 (MQ=255)
                       gTTCAGCGACGA‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCTGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAacac                                                                                                                              <  4:105034/148‑1 (MQ=255)
                          cAGCGACGA‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCa                                                                                                                           >  6:271416/1‑148 (MQ=255)
                           aGCGACGA‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCata                                                                                                                         >  2:357633/1‑149 (MQ=255)
                             cgacga‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATAtt                                                                                                                       >  7:63058/1‑149 (MQ=255)
                                cga‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTg                                                                                                                    <  1:303112/149‑1 (MQ=255)
                                  a‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGa                                                                                                                  <  3:346284/149‑1 (MQ=255)
                                  a‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGa                                                                                                                  <  6:122004/149‑1 (MQ=255)
                                  a‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGa                                                                                                                  <  6:226602/149‑1 (MQ=255)
                                    tCCTTCACTCGGCGGTGATAAAGTGCTCTTGGCAATCCAGGCCGCCTGTATGGAGTAAGCGGACTATGCTTCACTTG‑AAGGTTAGGCAATCCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGa                                                                                                                  <  1:320142/149‑1 (MQ=255)
                                      cGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGt                                                                                                                <  8:278114/148‑1 (MQ=255)
                                      cGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGt                                                                                                                <  8:192891/148‑1 (MQ=255)
                                       gTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGttt                                                                                                              >  6:31103/1‑149 (MQ=255)
                                        tCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGtttt                                                                                                             >  1:268729/1‑149 (MQ=255)
                                         cACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGa                                                                                                                  >  7:179117/1‑143 (MQ=255)
                                            gcggcgGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCAt                                                                                                         <  6:218608/149‑1 (MQ=255)
                                              ggcggTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATgg                                                                                                       <  1:151748/149‑1 (MQ=255)
                                               gcggTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGc                                                                                                      <  1:357633/149‑1 (MQ=255)
                                                  gTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGt                                                                                                   <  7:363240/149‑1 (MQ=255)
                                                       aaaGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTNATATTCTGGAGTTTTCCATGGCTGTTCAAc                                                                                              <  1:34490/149‑1 (MQ=255)
                                                       aaaGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAAc                                                                                              <  1:28397/149‑1 (MQ=255)
                                                             cTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACt                                                                                         <  2:392828/148‑1 (MQ=255)
                                                              tCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACttt                                                                                       >  1:182817/1‑149 (MQ=255)
                                                              tCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACttt                                                                                       >  5:372800/1‑149 (MQ=255)
                                                              tCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACttt                                                                                       >  6:70999/1‑149 (MQ=255)
                                                              tCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACtt                                                                                        >  8:191174/1‑148 (MQ=255)
                                                                gAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCGCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAGCACTCATATTCTGGAGTTTTCCATGGCTGTTCAGCGTACTTtct                                                                                     >  4:368843/1‑149 (MQ=255)
                                                                gAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTtct                                                                                     >  2:115440/1‑149 (MQ=255)
                                                                 aGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTtctc                                                                                    <  7:246200/149‑1 (MQ=255)
                                                                 aGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTtctc                                                                                    <  2:300302/149‑1 (MQ=255)
                                                                 aGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTtctc                                                                                    <  2:184189/149‑1 (MQ=255)
                                                                    cAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCcat                                                                                 <  5:307228/149‑1 (MQ=255)
                                                                          aGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCatcatc                                                                             >  6:97162/1‑147 (MQ=255)
                                                                          aGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCatcatc                                                                             >  2:74402/1‑147 (MQ=255)
                                                                          aGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCaa                                                                           <  2:237347/149‑1 (MQ=255)
                                                                          aGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCaa                                                                           <  5:31103/149‑1 (MQ=255)
                                                                            gccgccTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTc                                                                                                                            >  7:87714/1‑98 (MQ=255)
                                                                            gccgccTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTc                                                                                                                            <  8:87714/98‑1 (MQ=255)
                                                                                 cTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTg                                                                                                    <  3:44814/117‑1 (MQ=255)
                                                                                 cTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTg                                                                                                    >  4:44814/1‑117 (MQ=255)
                                                                                  tGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACg                                                                   >  3:368258/1‑149 (MQ=255)
                                                                                   ggATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGc                                                                  >  1:381745/1‑149 (MQ=255)
                                                                                   ggATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGc                                                                  >  4:351403/1‑149 (MQ=255)
                                                                                     aTCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTg                                                                >  3:97396/1‑149 (MQ=255)
                                                                                        gacgaAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTAc                                                                                          <  8:317247/120‑1 (MQ=255)
                                                                                        gacgaAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTAc                                                                                          >  7:317247/1‑120 (MQ=255)
                                                                                         acgaAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGc                                                            >  6:358161/1‑149 (MQ=255)
                                                                                         acgaAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGc                                                            >  4:387516/1‑149 (MQ=255)
                                                                                            aaGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCa                                                          >  4:123160/1‑148 (MQ=255)
                                                                                                cGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAAc                                                     <  5:97162/149‑1 (MQ=255)
                                                                                                 gACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACg                                                    <  5:358161/149‑1 (MQ=255)
                                                                                                 gACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACg                                                    >  5:232800/1‑149 (MQ=255)
                                                                                                  aCTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGt                                                   <  5:271416/149‑1 (MQ=255)
                                                                                                  aCTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGt                                                   >  2:246989/1‑149 (MQ=255)
                                                                                                         gTCACTTG‑GAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGc                                            <  8:149670/149‑1 (MQ=255)
                                                                                                                g‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATca                                     >  2:362200/1‑149 (MQ=255)
                                                                                                                  cAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATcac                                    <  2:268729/149‑1 (MQ=255)
                                                                                                                    tgTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATcacca                                  <  5:190879/148‑1 (MQ=255)
                                                                                                                    ggTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATcacca                                  <  3:123160/149‑1 (MQ=255)
                                                                                                                           cAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCtt                           <  4:368258/149‑1 (MQ=255)
                                                                                                                           cAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCtt                           >  1:319367/1‑149 (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)
                                                                                                                                                     tCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCGGCCTGAAAATCGtt  >  2:193298/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGA‑TCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTG‑CAGGTTAGGCAAT‑CCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCGGCCTGAAAATCGTT  >  NC_002947/983441‑983735

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: ‑