Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,348,803 C→T 5.1% intergenic (+35/+19) yfaE → / ← inaA ferredoxin involved with ribonucleotide reductase diferric‑tyrosyl radical (Y*) cofactor maintenance/acid‑inducible Kdo/WaaP family putative kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,348,8030CT 5.1% 257.8 / 2.2 78intergenic (+35/+19)yfaE/inaAferredoxin involved with ribonucleotide reductase diferric‑tyrosyl radical (Y*) cofactor maintenance/acid‑inducible Kdo/WaaP family putative kinase
Reads supporting (aligned to +/- strand):  ref base C (36/38);  new base T (2/2);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

CGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAGCGGCGACGACTTTTC  >  NC_000913/2348655‑2348948
                                                                                                                                                    |                                                                                                                                                  
cGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTc                                                                                                                                                    <  7:81333/149‑1 (MQ=255)
 gCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCt                                                                                                                                                   <  5:257959/149‑1 (MQ=255)
   cTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTa                                                                                                                                                 >  4:198107/1‑149 (MQ=255)
      gTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATc                                                                                                                                               <  3:152229/148‑1 (MQ=255)
       ttGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTg                                                                                                                                             >  5:204766/1‑149 (MQ=255)
         gCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGc                                                                                                                                           >  2:140402/1‑149 (MQ=255)
          cAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGcc                                                                                                                                         <  6:129008/149‑1 (MQ=255)
          cAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGc                                                                                                                                           <  8:113152/148‑1 (MQ=255)
          cAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGc                                                                                                                                           <  3:303761/148‑1 (MQ=255)
          cAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATAAGACGCATAAGCGTTGTATCTGGca                                                                                                                                          <  6:279444/149‑2 (MQ=39)
           aGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGcc                                                                                                                                         >  7:177992/1‑148 (MQ=255)
             gTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTa                                                                                                                                       >  7:153868/1‑148 (MQ=255)
              tCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTAc                                                                                                                                      >  1:303669/1‑148 (MQ=255)
                aaGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAAc                                                                                                                                   >  4:258303/1‑149 (MQ=255)
                       cTGGATTGCCGAACCGTTAGCCTTTATTAAGCCGGGGGAAATTTTGCCCAGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACa                                                                                                                            >  2:315340/1‑149 (MQ=255)
                        tGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACAt                                                                                                                           >  8:97375/1‑149 (MQ=255)
                          gATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACAtcg                                                                                                                         >  4:350701/1‑147 (MQ=255)
                          gATTGCCGAACCGTTAGCCTTTATTCAGCCGCGGTAAATTTTGCCCTGTTGTGGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACAt                                                                                                                           <  3:191521/147‑1 (MQ=255)
                               ccGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgtag                                                                                                                    <  8:153868/149‑1 (MQ=255)
                                cGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgtagt                                                                                                                   <  5:219329/149‑1 (MQ=255)
                                  aaCCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTaa                                                                                                                 >  5:238281/1‑149 (MQ=255)
                                    ccGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCACTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTTGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATa                                                                                                               <  4:202190/149‑1 (MQ=255)
                                     cGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTCCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATAc                                                                                                              <  5:337568/149‑1 (MQ=255)
                                     cGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGGGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATAc                                                                                                              <  7:336726/149‑1 (MQ=255)
                                       ttAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGc                                                                                                            <  8:177992/149‑1 (MQ=255)
                                            cTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATgg                                                                                                                         >  3:278150/1‑131 (MQ=255)
                                            cTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATgg                                                                                                                         <  4:278150/131‑1 (MQ=255)
                                            cTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCa                                                                                                       >  8:17725/1‑149 (MQ=255)
                                             tttATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGAtt                                                                                                                               <  4:221025/124‑1 (MQ=255)
                                             tttATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGAtt                                                                                                                               >  3:221025/1‑124 (MQ=255)
                                                                ttttGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTg                                                                                   <  3:350701/149‑1 (MQ=255)
                                                                  ttGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTgg                                                                                  <  1:140402/148‑1 (MQ=255)
                                                                  ttGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTgg                                                                                  <  8:215823/148‑1 (MQ=255)
                                                                      ccTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGc                                                                                                                                           >  6:372862/1‑88 (MQ=255)
                                                                      ccTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGc                                                                                                                                           <  5:372862/88‑1 (MQ=255)
                                                                       cTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGGTGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGAAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTg                                                                            >  3:246196/1‑149 (MQ=255)
                                                                       cTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgta                                                                                                                     <  5:178583/108‑1 (MQ=255)
                                                                       cTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgta                                                                                                                     >  6:178583/1‑108 (MQ=255)
                                                                        tgttgtTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTgg                                                                                  <  7:261431/142‑1 (MQ=255)
                                                                        tgttgtTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTgg                                                                                  >  8:261431/1‑142 (MQ=255)
                                                                             ttgCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCa                                                                      <  1:315340/149‑1 (MQ=255)
                                                                               gCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTg                                                                            >  6:248652/1‑141 (MQ=255)
                                                                               gCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTg                                                                            <  5:248652/141‑1 (MQ=255)
                                                                                       aaaGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTc                                                            <  3:198107/149‑1 (MQ=255)
                                                                                       aaaGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTc                                                            <  2:303669/149‑1 (MQ=255)
                                                                                        aaGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTc                                                            >  3:195907/1‑148 (MQ=255)
                                                                                          ggCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAg                                                         >  7:132313/1‑149 (MQ=255)
                                                                                              atatTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTgg                                                     >  8:232946/1‑149 (MQ=255)
                                                                                                 ttGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATAtt                                                                >  5:76360/1‑135 (MQ=255)
                                                                                                 ttGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATAtt                                                                <  6:76360/135‑1 (MQ=255)
                                                                                                   gAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgta                                                                                                                     >  5:253255/1‑80 (MQ=255)
                                                                                                   gAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgta                                                                                                                     <  6:253255/80‑1 (MQ=255)
                                                                                                    aaaTCGAGGTGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGAAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGaa                                                <  4:246196/148‑1 (MQ=255)
                                                                                                    aaaTCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAg                                               >  7:26950/1‑149 (MQ=255)
                                                                                                         gagaTGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCAt                                           <  8:26950/148‑1 (MQ=255)
                                                                                                         gagaTGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATg                                          <  3:258303/149‑1 (MQ=255)
                                                                                                          agaTGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGa                                         <  8:132313/149‑1 (MQ=255)
                                                                                                          agaTGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGa                                         >  3:228125/1‑149 (MQ=255)
                                                                                                                 aaTTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCg                                  >  5:289951/1‑149 (MQ=255)
                                                                                                                   ttGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTg                                                                                   >  2:240747/1‑98 (MQ=255)
                                                                                                                   ttGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTg                                                                                   <  1:240747/98‑1 (MQ=255)
                                                                                                                      gTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGAt                                        >  1:381969/1‑138 (MQ=255)
                                                                                                                      gTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCACCTGGCGGAAGTCATGAt                                        <  2:381969/138‑1 (MQ=255)
                                                                                                                      gTTGTAGGCCAGATAAGACGCATAAGCGTTGTATCTGGCACTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCa                                                                      >  5:100892/1‑108 (MQ=18)
                                                                                                                      gTTGTAGGCCAGATAAGACGCATAAGCGTTGTATCTGGCACTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCa                                                                      <  6:100892/108‑1 (MQ=18)
                                                                                                                       ttGTAGGCCAGATAAGACGCATAAGCGTTGTATCTGGCACTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTAt                            >  4:281283/1‑149 (MQ=37)
                                                                                                                            gTGCCAGATACAGCGGGTATGCGTCTTAGCTGGGCTACAACGATTACATGGCGTAGTAATACGCTTTCGCCTGCTGCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGa                       >  8:317997/1‑149 (MQ=255)
                                                                                                                            gTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGa                       <  7:232946/149‑1 (MQ=255)
                                                                                                                                cAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGc                    >  8:275128/1‑148 (MQ=255)
                                                                                                                                       aaCGCTTATGCGTCTTATCTGGCCTACAACGATTACCTGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAgcggcg            >  2:49177/1‑149 (MQ=255)
                                                                                                                                       aaCGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgta                                                                                                                     <  8:12372/44‑1 (MQ=255)
                                                                                                                                       aaCGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGCCGCAAgcggcg            >  6:52729/1‑149 (MQ=255)
                                                                                                                                         cGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAGCGGcgac          <  7:275128/149‑1 (MQ=255)
                                                                                                                                           cTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAGCGGcgacga        <  7:17725/149‑1 (MQ=255)
                                                                                                                                            ttATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCaa                  <  7:86741/138‑1 (MQ=255)
                                                                                                                                            ttATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCaa                  >  8:86741/1‑138 (MQ=255)
                                                                                                                                                gCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGc                                                                                                  <  6:79089/54‑1 (MQ=255)
                                                                                                                                                gCGTCTTATCTGCCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGc                                                                                                  >  5:79089/1‑54 (MQ=255)
                                                                                                                                                 cGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAGCGGCGACGACTTTtc  >  6:133464/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTCACCTGCTCCCAGTCCGCTTTGGGGATTGGCTCCAGATATTTTTCCAGCTGGCGGAAGTCATGATTAATCGCTTTATCGCGACGCAAGCGGCGACGACTTTTC  >  NC_000913/2348655‑2348948

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