Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,468 A→T 5.8% E517V (GAA→GTA) ‡ rsxC → SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,708,4680AT 5.8% 246.4 / 5.2 69E517V (GAA→GTA) ‡rsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base A (37/28);  new base T (2/2);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

GCGAAAAAGCGGCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGC  >  NC_000913/1708324‑1708615
                                                                                                                                                |                                                                                                                                                   
gcgAAAAAGCGGCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACg                                                                                                                                                  <  2:21054/148‑1 (MQ=255)
  gAAAAAGCGGCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGccc                                                                                                                                               <  4:682628/149‑1 (MQ=255)
  gAAAAAGCGGCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGcc                                                                                                                                                <  2:499421/148‑1 (MQ=255)
   aaaaaGCGGCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGccc                                                                                                                                               >  2:582450/1‑148 (MQ=255)
   aaaaaGCGGCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGccc                                                                                                                                               >  1:741036/1‑148 (MQ=255)
    aaaaGCGGCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGccc                                                                                                                                               >  4:222331/1‑147 (MQ=255)
    aaaaGCGGCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGccc                                                                                                                                               <  3:222331/147‑1 (MQ=255)
         cGGCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATaa                                                                                                                                         <  1:491618/148‑1 (MQ=255)
           gCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAg                                                                                                                                      >  3:619815/1‑149 (MQ=255)
                     aaCGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAattat                                                                                                                             >  2:441897/1‑148 (MQ=255)
                       cGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGTTCGCCCgcttta                                                                                                                                         >  1:102140/1‑129 (MQ=255)
                       cGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGTTCGCCCgcttta                                                                                                                                         <  2:102140/134‑6 (MQ=255)
                       cGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTgc                                                                                                                          <  3:197285/149‑1 (MQ=255)
                       cGACATAAGAGCGCAGCCGTTCAACCTGAAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTgc                                                                                                                          <  3:871656/149‑1 (MQ=255)
                        gACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTgca                                                                                                                         <  3:84100/149‑1 (MQ=255)
                           aTAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGTTCGCCCgct                                                                                                                                            <  3:832549/127‑3 (MQ=255)
                           aTAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGTTCGCCCgct                                                                                                                                            >  4:832549/1‑125 (MQ=255)
                             aaGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCTCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTTCAATTATTGCAGCAc                                                                                                                     >  2:263887/1‑148 (MQ=255)
                                   gcAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGc                                                                                                              >  4:849396/1‑149 (MQ=255)
                                     aGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGccc                                                                                                            >  3:758323/1‑149 (MQ=255)
                                         gTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTaa                                                                                                        >  4:710164/1‑149 (MQ=255)
                                            cAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCTCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGCATAACAGTGCAATTATTGCAGTACGGGAAGCCCGTAAAg                                                                                                      <  1:263887/148‑1 (MQ=255)
                                            cAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAg                                                                                                      <  1:411786/148‑1 (MQ=255)
                                                cTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCaa                                                                                                 <  4:394328/149‑1 (MQ=255)
                                                 tGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAat                                                                                                                                >  4:415865/1‑117 (MQ=255)
                                                 tGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAat                                                                                                                                <  3:415865/117‑1 (MQ=255)
                                                      ccAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCaga                                                                                           >  4:427108/1‑149 (MQ=255)
                                                       cAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCagag                                                                                          <  4:323419/149‑1 (MQ=255)
                                                         aagataaagatGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGc                                                                                         >  2:67008/1‑148 (MQ=255)
                                                             taaagatGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGaaa                                                                                     >  1:907802/1‑148 (MQ=255)
                                                               aagatGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAg                                                                                  <  4:215851/149‑1 (MQ=255)
                                                                   tGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCa                                                                               >  2:129954/1‑148 (MQ=255)
                                                                     cGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGa                                                                             <  2:614175/148‑1 (MQ=255)
                                                                     cGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGa                                                                             <  1:441897/148‑1 (MQ=255)
                                                                          gctgcCGCTCTGGCGCGGGTGGAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTgca                                                                       <  4:171564/149‑1 (MQ=255)
                                                                          gctgcCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTgc                                                                        <  1:129954/148‑1 (MQ=255)
                                                                             gcCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTgcagca                                                                    >  4:203469/1‑149 (MQ=255)
                                                                                 ctctGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACt                                                                >  3:54997/1‑149 (MQ=255)
                                                                                  tctGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACt                                                                >  2:39804/1‑148 (MQ=255)
                                                                                       cgcgGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAacga                                                           <  2:907802/148‑1 (MQ=255)
                                                                                        gcgGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAacgac                                                          <  1:2362/148‑1 (MQ=255)
                                                                                          gggTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACgc                                                        <  1:39804/148‑1 (MQ=255)
                                                                                            gTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACgcagc                                                     >  4:610151/1‑149 (MQ=255)
                                                                                                aaGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAAc                                                  >  2:219237/1‑148 (MQ=255)
                                                                                                  gagaAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCg                                                >  2:192323/1‑148 (MQ=255)
                                                                                                    gaAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGtt                                              >  1:201621/1‑148 (MQ=255)
                                                                                                    gaAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGtt                                              >  2:317293/1‑148 (MQ=255)
                                                                                                     aaaCAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTg                                             <  1:582450/148‑1 (MQ=255)
                                                                                                      aaCAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCt                                           >  3:428698/1‑149 (MQ=255)
                                                                                                        cAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGa                                         <  4:54997/149‑1 (MQ=255)
                                                                                                         aGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGAt                                        <  3:849396/149‑1 (MQ=255)
                                                                                                               aGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTa                                                               >  2:809931/1‑120 (MQ=255)
                                                                                                               aGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTa                                                               <  1:809931/120‑1 (MQ=255)
                                                                                                               aGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGt                                  >  3:442382/1‑149 (MQ=255)
                                                                                                                    acacAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAAc                             >  3:875617/1‑149 (MQ=255)
                                                                                                                           cTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGt                       <  1:67008/148‑1 (MQ=255)
                                                                                                                              ttGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGaa                   >  4:530218/1‑149 (MQ=255)
                                                                                                                               tgtgATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGaa                   >  1:803587/1‑148 (MQ=255)
                                                                                                                                gtgATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAag                  >  1:653908/1‑148 (MQ=255)
                                                                                                                                 tgATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAagca                >  3:731316/1‑149 (MQ=255)
                                                                                                                                  gATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAagcag               <  3:710164/149‑1 (MQ=255)
                                                                                                                                    ttAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTAGAagcagc              >  2:786391/1‑148 (MQ=255)
                                                                                                                                    ttAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGACAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCt             >  4:536145/1‑149 (MQ=255)
                                                                                                                                     tAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCt             >  2:546894/1‑148 (MQ=255)
                                                                                                                                     tAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCt             >  2:409655/1‑148 (MQ=255)
                                                                                                                                         gCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGc        >  4:902633/1‑149 (MQ=255)
                                                                                                                                            ggCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTaaaa                              >  4:221602/1‑124 (MQ=255)
                                                                                                                                            ggCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTaaaa                              <  3:221602/124‑1 (MQ=255)
                                                                                                                                               gAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCgcgc  >  4:577830/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCGAAAAAGCGGCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGC  >  NC_000913/1708324‑1708615

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: