Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,082,474 G→C 6.5% R96G (CGC→GGC)  fdoH ← formate dehydrogenase O subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,082,4740GC 6.5% 174.1 / 3.8 69R96G (CGC→GGC) fdoHformate dehydrogenase O subunit beta
Reads supporting (aligned to +/- strand):  ref base G (37/27);  new base C (3/2);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

AGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAA  >  NC_000913/4082333‑4082613
                                                                                                                                             |                                                                                                                                           
aGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGAAGAAACGGTGTGTTTAGTGAATAATTGGAGCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGGTGGGAAGGGTGCATACAGCCGGCTGTGCgg                                                                                                                                            >  1:211555/1‑143 (MQ=255)
  cccGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGGTACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATCCAGCCGTCTTTGCGGAt                                                                                                                                          >  1:1075970/1‑143 (MQ=255)
  cccGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGAt                                                                                                                                          >  1:748038/1‑143 (MQ=255)
  cccGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAACCGTCTTTGCGGAt                                                                                                                                          >  1:1087200/1‑143 (MQ=255)
    ctccttgcaatAACCGGATCAAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCa                                                                                                                                        <  1:633552/138‑1 (MQ=255)
    cGCAAGGCAATAACCGCAGCCAATGCACTGCTCGGCCGGGAAGTCGACGATACCGTTGGCCTGTTGGTTCTTACCACGTTCTGTGGGGTCCACCCTCAGACAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCa                                                                                                                                        <  1:524239/143‑1 (MQ=255)
       aatggaatAACTGCAGCCAATGCAGTGCTCGGACTGGAAGGCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGcc                                                                                                                                     <  1:279383/143‑1 (MQ=255)
       aatgcaatAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCCCGCTTTCCGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCCGATCAGcc                                                                                                                                     >  1:113542/1‑143 (MQ=255)
       aatgcaatAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCCGATCAGcc                                                                                                                                     >  1:63729/1‑143 (MQ=255)
       aatgcaatAACCGCAGCCAATGCACTGCACGGACTGGACGTCGACGATACAGTTGGAACACTGAATAATTGCCCCTTCCGCCGGGCACGCTTTCAGGCCGCCTGGCTCGGAACAGTGCATACACCCGTCTTTGCGGATCAggc                                                                                                                                     >  1:183171/1‑141 (MQ=255)
          gcaggcaCCGCCGCCTCTGCACTGGTCGGATTGGAAGTCGACGATACCGTGGGCATACTGAATAATTGCACCTGCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCAtt                                                                                                                                  <  2:840573/137‑1 (MQ=25)
                 ccGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGAATCAACCATTCCAGatt                                                                                                                           >  1:891853/1‑143 (MQ=255)
                 ccGCAGCCAATGCACTGCTCGGACTGAAAGTCGACGATACCGTTGGCATACGGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGGACAGTGCATACAGCCGTCTTTGCGCATCAGCCCTGCCAGttt                                                                                                                           >  1:120971/1‑143 (MQ=37)
                    cggCCAAGGCACTGCTCGGACTGGAAGTCGACGAGACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTc                                                                                                                        <  2:479130/141‑1 (MQ=255)
                    cAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCCGTTTGTc                                                                                                                        >  2:157095/1‑143 (MQ=255)
                     aGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTACATACAGACGTCATTGCGCATCAGCTATACCAGTATATCg                                                                                                                       >  1:436366/1‑143 (MQ=255)
                        cAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCGTCCGCCGGGCACGCGTTCAGGCAGCCTGGATCGGGACAGTGCAGACAGCGGGCTTTGCGGAACAGCCGTGCCGGGTGGTCGtcc                                                                                                                    >  2:364632/1‑141 (MQ=255)
                        cAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTCAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCCTACAGCCGTCTTCGCCGATCCGCCATTCCAGTTTGGCGGTc                                                                                                                    >  2:734445/1‑143 (MQ=255)
                          aTGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCTCGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGGCGTCTTTGCGGGTCAGCACGTCAAGTTTGTGGtatta                                                                                                                  >  2:351880/1‑139 (MQ=255)
                           tGCACTGCTCGGACTGGAAGTAGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCGCGCTTTCTGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCCTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGc                                                                                                                 >  1:991219/1‑143 (MQ=255)
                                tGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCAc                                                                                                            >  1:975793/1‑143 (MQ=255)
                                   tCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTc                                                                                                         <  2:748038/143‑1 (MQ=255)
                                   gcGGACTTGAAGTCTACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTc                                                                                                         <  2:785923/142‑1 (MQ=255)
                                      gACTGTACGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCga                                                                                                      <  2:1059654/143‑1 (MQ=255)
                                      gACTGTACGTCGACGATACCGCTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCGGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCga                                                                                                      <  2:867445/143‑1 (MQ=255)
                                         tGGAAGTCGAGGATACGGGTGGCATACGGAATAGTGGCACCTTCCGCCGGGCACGCTGTGAGGCAGCCCGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGGTTGTCGTTCTGCTCCACTTCCGAGaa                                                                                                   <  2:1043064/143‑1 (MQ=37)
                                          ggAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATGGGAACAGTCCATACAGCCGTCTTTGCGGGCGGCCCACTCCACTCCGTCGTTCTGCTCGACTCCCgagctg                                                                                                  >  2:232949/1‑140 (MQ=255)
                                           gAAGTCGCCGATACCGTTGGCATACTGGATAATTGCCCCTTCCGCCGTGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAgc                                                                                                 <  2:124169/143‑1 (MQ=255)
                                                             ggCATACTGAATGATTGCACCTTCCGCCGGGCACGCGTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCCGATCGGCCATTCCAGTGTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGAtt                                                                               <  2:113542/143‑1 (MQ=37)
                                                             ggCATACTGAATAATGGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCCGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGAtt                                                                               <  2:63729/143‑1 (MQ=255)
                                                               cATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCCCGATTTg                                                                             >  1:237901/1‑143 (MQ=255)
                                                                 tACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTAGATCGGAACAGTGCATACAGCCGTCTTTGAGGAACATCCATTCCAGTTCGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCGCGCTTTGgc                                                                           >  1:705838/1‑143 (MQ=37)
                                                                  acggcataatTGCACGTTCCCCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCGGCCATTCCAGTGTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGCTTTGgcg                                                                          <  2:407134/138‑1 (MQ=37)
                                                                   cTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGCACAGTGCATAGAGCCGCCTCTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCGCCACTTCCGAGGAGCGCATTCGCGTCCCCCGTTGGgcgc                                                                         >  2:134172/1‑143 (MQ=255)
                                                                       ataatTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTCTCATTATACTCCACTTACTAGAAGCGCATTACCGTCAACGATTTGGCGCTTaa                                                                     >  1:651155/1‑143 (MQ=38)
                                                                          atTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTGCAGGATTTGGCGCTTAACTc                                                                  >  1:408898/1‑143 (MQ=255)
                                                                            tGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACATCCGTCTTTGCGGTTCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGACGCGCATTAACATCCACGATTTCGCGCCTAAACCAt                                                                >  2:564006/1‑143 (MQ=255)
                                                                             gCACCTTCCGCCGGGCACGCTTTCAGGTAGCGTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCAtt                                                               <  2:17888/143‑1 (MQ=255)
                                                                              cacgtcccgccgGGCACGCTTCCAGGCGGCCTGGAGCGGGCCGGTGCACACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCGCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTg                                                              <  1:900345/137‑1 (MQ=255)
                                                                                cggtgcgccgGGCACCGTCTCAGGGAGCCTGGATCGGCACACTGCCTACAGCCGCCTTTGCGGATCGTCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTggg                                                            <  2:733668/138‑1 (MQ=255)
                                                                                ccTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTGTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCCTTggg                                                            >  1:887173/1‑143 (MQ=255)
                                                                                 cTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCCTTgggg                                                           >  1:761921/1‑143 (MQ=255)
                                                                                         gggCACGCTCTCGGGCGGCCCGCACCGGCACAGTGCATACAGCCGTCTTTGCGGATCGGCCATTCCAGTTCGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTa                                                   <  2:211555/143‑1 (MQ=38)
                                                                                           gCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTaca                                                 >  2:345366/1‑143 (MQ=255)
                                                                                               gCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACAcccc                                             >  1:502107/1‑143 (MQ=255)
                                                                                                        cAGCCTGGATCGGAACATTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTAtt                                    <  1:960684/143‑1 (MQ=255)
                                                                                                        cAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTAtt                                    <  1:644315/143‑1 (MQ=255)
                                                                                                          gccctgATCGACACAGTACATGCATCCGTGGTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGc                                  <  1:1061202/138‑1 (MQ=255)
                                                                                                            cTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTCGGGTTATCGTACACCCCAATATTATTCCCg                                >  1:3592/1‑143 (MQ=255)
                                                                                                             tGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGACGCGCATTACCGTCCACGATTTGGCGCGTACATCCTTGGTGTTGTCGTACACCCCCTTGGTTTTGCCgc                               >  1:55939/1‑142 (MQ=37)
                                                                                                               gATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATGCCAGTTTGTCGTTCTGCTCCTCTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTTTCGTACACCCCAATGTTATTGCCGACg                             <  2:571881/143‑1 (MQ=255)
                                                                                                                aTCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTTTGCTCCACTTCCGAGAAGCGCATTACCGTCCCCGACATGGCGCTTAAAACAATGACGTGCACGTCCAGCCCAAACGTATTGCCGACgg                            >  2:424044/1‑143 (MQ=25)
                                                                                                                   ggAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACCCCCCCCTGTTATTGCCGACGGTAt                         >  1:969076/1‑143 (MQ=255)
                                                                                                                   ggAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTAt                         >  1:960203/1‑143 (MQ=255)
                                                                                                                     aaCAGTGCATACAGCCGTCTTTGCGGATCAGCCATGCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATcg                       <  1:795685/143‑1 (MQ=255)
                                                                                                                         gTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGCTTTGGCGCCAAAATCATTGGGGATGTCGTACACCCCAATGTTATAGGCGTCGGAACCGCGaa                   >  2:902086/1‑143 (MQ=255)
                                                                                                                         gTGCATACAGCCGTCTGTGCGGATCAGCCATTCCAGTGTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGaa                   <  1:164842/143‑1 (MQ=255)
                                                                                                                          tGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACTGTATCGCGAAt                  >  1:620317/1‑143 (MQ=255)
                                                                                                                          tGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGCAt                  >  1:880160/1‑143 (MQ=255)
                                                                                                                          tGCACACAGCCGGCTGTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCTTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAAt                  <  1:683204/143‑1 (MQ=255)
                                                                                                                            cATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGt                >  1:14072/1‑143 (MQ=255)
                                                                                                                             aTACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGCATCGCGAATGTc               >  1:501014/1‑143 (MQ=255)
                                                                                                                             aTACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGACGCTTAAATCATTGGTGTTGTCGTACACCCCAATGTTATAACAGACGATAACGCGacaatc               >  2:468929/1‑138 (MQ=255)
                                                                                                                              tACAGCCGTCTGTGCGGATCAGCCATCCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCg              <  1:440093/143‑1 (MQ=255)
                                                                                                                              tACAGCCGTCGTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCACTGTTGTTGCCGCCGGTCTCGCGGGTGTCg              >  1:171765/1‑143 (MQ=255)
                                                                                                                                   ccGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGACTGTCGTTGCa         >  1:243929/1‑143 (MQ=255)
                                                                                                                                    cGTCTTGGCGGATCAGCCAGTCCAGTTTGGCGTGCTGCGCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCAc        <  1:115313/143‑1 (MQ=255)
                                                                                                                                         tggCGGATCCGCCGTTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGa   <  2:1027667/141‑1 (MQ=255)
                                                                                                                                          tGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGaa  <  1:157095/143‑1 (MQ=255)
                                                                                                                                          tGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCAGACGGTATCGCGCATGTCGTTCCACTCTGaa  >  1:612324/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                           
AGCCCGCAATGCAATAACCGCAGCCAATGCACTGCTCGGACTGGAAGTCGACGATACCGTTGGCATACTGAATAATTGCACCTTCCGCCGGGCACGCTTTCAGGCAGCCTGGATCGGAACAGTGCATACAGCCGTCTTTGCGGATCAGCCATTCCAGTTTGTCGTTCTGCTCCACTTCCGAGAAGCGCATTACCGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAA  >  NC_000913/4082333‑4082613

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