Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,945,381 C→A 6.9% Q54H (CAG→CAT)  PP_5212 ← iron‑sulfur‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,945,3810CA 6.9% 223.1 / 4.1 70Q54H (CAG→CAT) PP_5212iron‑sulfur‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (35/30);  new base A (2/3);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

CGTCGACGCACTCACTGACGCTGCACGCCAGCTTGCGCACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCATCCAGGATCCTTTCCCCGGGTTCAAGCGCCAGCACCGCCCCGGACGGCTGC  >  NC_002947/5945236‑5945528
                                                                                                                                                 |                                                                                                                                                   
cGTCGACGCACTCACTGACGCTGCACGCCAGCTTGCGCACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTgg                                                                                                                                                  <  7:98505/149‑1 (MQ=255)
       gCACTCACTGACGCTGCACGCCAGCTTGCGCACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGG‑GTTCGCCAGGATGGCGaagtcc                                                                                                                                          <  5:70249/149‑5 (MQ=255)
             actGACGCTGCACGCCAGCTTGCGCACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCAGCCCtt                                                                                                                                     <  5:15513/149‑1 (MQ=255)
             actGACGCTGCACGCCAGCTTGCGCACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCAGCCCtt                                                                                                                                     >  6:15513/1‑149 (MQ=255)
              ctGACGCTGCACGCCAGCTTGCGCACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACACCTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTc                                                                                                                                    <  5:6456/149‑1 (MQ=255)
              ctGACGCTGCACGCCAGCTTGCGCACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGATCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTc                                                                                                                                    <  5:2148/149‑1 (MQ=255)
                 aCGCTGCACGCCAGCTTGCGCACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTcaac                                                                                                                                 >  1:47935/1‑149 (MQ=255)
                     tGCACGCCAGCTTGCGCACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTgg                                                                                                                                                  >  8:61212/1‑128 (MQ=255)
                     tGCACGCCAGCTTGCGCACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTgg                                                                                                                                                  <  7:61212/128‑1 (MQ=255)
                     tGCACGCCAGCTTGCGCACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTcaaccaa                                                                                                                              >  4:189297/1‑148 (MQ=255)
                        aCGCCAGCTTGCGCACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCaacaac                                                                                                                          <  7:149792/149‑1 (MQ=255)
                         cGCCAGCTTGCGCACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCAGGATGGCGaagt                                                                                                                                            <  8:237049/130‑3 (MQ=255)
                         cGCCAGCTTGCGCACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCAGGATGGCGaagt                                                                                                                                            >  7:237049/1‑128 (MQ=255)
                            cAGCTTGCGCACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACgcgg                                                                                                                      >  6:300365/1‑149 (MQ=255)
                             aGCTTGCGCACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACgcggc                                                                                                                     >  8:64395/1‑149 (MQ=255)
                                ttGCGCACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGACCCCGCCCTTAAACCAACAACGCGgcgc                                                                                                                   >  7:192285/1‑148 (MQ=255)
                                ttGCGCACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGgcgc                                                                                                                   >  2:277003/1‑148 (MQ=255)
                                 tGCTCACTGGCAGTTCTCCCAAGGCAAGAACAACTTCCCAGAGCAATAGTAATTCCTCCAGAGGTTGTGAAATGAATTTTAACAACTCGCCATGGTCGCGTAGTTCTCCATCATGTAGAACCAGTAGTTAAACCAACAAagcttcgta                                                                                                                  <  4:107944/148‑7 (MQ=255)
                                    gcACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAAcgcg                                                                                                                       <  1:189643/140‑1 (MQ=255)
                                    gcACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAAcgcg                                                                                                                       >  2:189643/1‑140 (MQ=255)
                                       cTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTAGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATgg                                                                                                           <  8:51360/149‑1 (MQ=255)
                                            aGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAAtt                                                                                                      <  2:281311/149‑1 (MQ=255)
                                            aGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAAtt                                                                                                      >  1:56663/1‑149 (MQ=255)
                                                cGCCCAAGGCAAGCACACCGTCCCAGAGCAACGCGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGc                                                                                                    <  3:7827/147‑1 (MQ=255)
                                                cGCCCAAGGCAAGCACACAGTCCCAGAGCAACGCGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGCCAATTGc                                                                                                    >  4:7827/1‑147 (MQ=255)
                                                   ccAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTg                                                                                               >  8:147481/1‑149 (MQ=255)
                                                    cAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGc                                                                                              >  4:298948/1‑149 (MQ=255)
                                                    cAAGGCAAGCACACCGGCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGc                                                                                              <  5:129810/149‑1 (MQ=255)
                                                         cAAGCACACCGTCCCAGAGCAACACGCAGGCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAg                                                                                         <  7:64395/149‑1 (MQ=255)
                                                          aaGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGc                                                                                        >  3:190324/1‑149 (MQ=255)
                                                             cacaCCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGt                                                                                     <  7:147481/149‑1 (MQ=255)
                                                                acCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTg                                                                                  >  7:29980/1‑149 (MQ=255)
                                                                acCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTg                                                                                  >  8:113093/1‑149 (MQ=255)
                                                                 ccGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTgg                                                                                 >  3:128003/1‑149 (MQ=255)
                                                                 ccGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTgg                                                                                 <  3:189297/149‑1 (MQ=255)
                                                                  cGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTAGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTggg                                                                                <  4:128003/149‑1 (MQ=255)
                                                                    tCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCa                                                                              <  7:113093/149‑1 (MQ=255)
                                                                          agCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATa                                                                        <  5:300365/149‑1 (MQ=255)
                                                                          agCAACACGCAGTCCTCCAGCGGTTAGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTTGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGATGTTTGGGCAGTCa                                                                          >  8:114118/1‑147 (MQ=255)
                                                                           gCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAg                                                                       <  6:309871/149‑1 (MQ=255)
                                                                                  gCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAgcc                                                                >  7:106400/1‑149 (MQ=255)
                                                                                         tccAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAgccg                                                               <  6:111290/143‑1 (MQ=255)
                                                                                         tccAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAgccg                                                               >  5:111290/1‑143 (MQ=255)
                                                                                         tccAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGc                                                         >  5:312423/1‑149 (MQ=255)
                                                                                         tccAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGc                                                         >  4:110684/1‑149 (MQ=255)
                                                                                           cAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGc                                                         <  8:106400/147‑1 (MQ=255)
                                                                                               ggTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCtccatc                                                   <  6:312423/149‑1 (MQ=255)
                                                                                                gTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCtccatcc                                                  >  7:243872/1‑149 (MQ=255)
                                                                                                         tGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCATCCAGGATCCtt                                         <  1:277003/149‑1 (MQ=255)
                                                                                                           cAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCATCCAGGATCCttt                                        <  3:298948/148‑1 (MQ=255)
                                                                                                            aGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGc                                                         <  7:67444/130‑1 (MQ=255)
                                                                                                            aGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGc                                                         >  8:67444/1‑130 (MQ=255)
                                                                                                              gTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCtccatc                                                   >  3:189752/1‑134 (MQ=255)
                                                                                                              gTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCtccatc                                                   <  4:189752/134‑1 (MQ=255)
                                                                                                                 aacaacTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCATCCAGGATCCTTTCCCCggg                                 >  8:84369/1‑149 (MQ=255)
                                                                                                                 aacaacTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCATCCAGGATCCTTTCCCCggg                                 >  4:124440/1‑149 (MQ=255)
                                                                                                                      cTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATATCCCAGCCGCCGTGCTCCATCCAGGATCCTTTCCCCGGGTGCa                             >  6:216847/1‑148 (MQ=255)
                                                                                                                       tCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCATCCAGGATCCTTTCCCCGGGTTCAAg                           >  7:22867/1‑149 (MQ=255)
                                                                                                                         gCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCATCCAGGATCCTTTCCCCGGGTTCAAgcg                         >  3:170779/1‑149 (MQ=255)
                                                                                                                         gCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCATCCAGGATCCTTTCCCCGGGTTCAAgc                          <  8:29980/148‑1 (MQ=255)
                                                                                                                           cATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCATCCAGGATCCTTTCCCCGGGTTCAAgcgc                        >  7:200991/1‑148 (MQ=255)
                                                                                                                           cATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCATCCAGGATCCTTTCCCCGGGTTCAAgcgc                        >  7:197452/1‑148 (MQ=255)
                                                                                                                           cATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCATCCAGGATCCTTTCCCCGGGTTCAAGCGcc                       >  3:271095/1‑149 (MQ=255)
                                                                                                                               gTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCATCCAGGATCCTTTCCCCGGGTTCAAGCGCCAGCa                   <  3:110684/149‑1 (MQ=255)
                                                                                                                                     gACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGAGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCATCCAGGATCCTTTCACCGGGTTCAAGCGACAGCACCGccc             >  8:175287/1‑149 (MQ=255)
                                                                                                                                       gttCGCCAGGATGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCATCCAGGATCCTTTCCCCGGGTTCAAGCGCCAGCACCGCCCCg           >  1:15176/2‑149 (MQ=255)
                                                                                                                                          cGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCATCCAGGATCCTTTCCCCGGGTTCAAGCGCCAGCACCGCCCCGGa         <  8:22867/148‑1 (MQ=255)
                                                                                                                                           gCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCATCCAGGATCCTTTCCCCGGGTTCAAGCGCCAGCACCGCCCCGGACg       >  4:77868/1‑149 (MQ=255)
                                                                                                                                            ccATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCATCCAGGATCCTTTCCCCGGGTTCAAGCGCCAGCACCGCCCCGGACgg      <  7:116284/149‑1 (MQ=255)
                                                                                                                                              aTCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCATCCAGTATCCTTTCCCCGGGTTCAAGCGCCAGCACCGCCCCGGACGGCt    <  4:190324/149‑1 (MQ=255)
                                                                                                                                                 cTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCATCCAGGATCCTTTCCCCGGGTTCAAGCGCCAGCACCGCCCCGGACGGCTGc  >  3:83722/1‑148 (MQ=255)
                                                                                                                                                 cTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCCGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCAGCCAGGATCCTTTCCCCGGGTTCAAGCGCCAGCACCGCCCCGGACCGCTGc  >  2:245916/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGTCGACGCACTCACTGACGCTGCACGCCAGCTTGCGCACTGGCAGTTCGCCCAAGGCAAGCACACCGTCCCAGAGCAACACGCAGTCCTCCAGCGGTTCGGCAATGCAGGTGAACAACTCGCCATGGTCGCGGACTTCGCCATCCTGGCGAACCCGCCCTTCAACCAACAACGCGGCGCAAACATGGCAATTGCCATTGCGGCAGCTGTTTGGGCAGTCATAGCCCAGCCGCCGTGCTCCATCCAGGATCCTTTCCCCGGGTTCAAGCGCCAGCACCGCCCCGGACGGCTGC  >  NC_002947/5945236‑5945528

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