Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,852,356 T→A 100% L74Q (CTG→CAG)  hypD → hydrogenase maturation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,852,3560TA98.8% 282.0 / ‑4.5 85L74Q (CTG→CAG) hypDhydrogenase maturation protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (57/27);  total (57/28)
Fisher's exact test for biased strand distribution p-value = 3.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

GCGGTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACTGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACA  >  NC_000913/2852260‑2852452
                                                                                                |                                                                                                
gCGGTCATACCCTCGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACTGcc                                                                                               <  1:89402/100‑1 (MQ=255)
   gTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGAt                                                                                            <  2:2001381/100‑1 (MQ=255)
    tCATACCAACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATg                                                                                           <  2:1243168/100‑1 (MQ=255)
     cATACCCACGCTATCGTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGCGTGCGTACAGCCGATg                                                                                           >  2:2319421/1‑99 (MQ=255)
      aTACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGTCCGGTGTGCGTACAGCCGATggg                                                                                         <  1:128202/100‑1 (MQ=255)
      aTACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATggg                                                                                         <  1:2326396/100‑1 (MQ=255)
      aTACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATggg                                                                                         <  1:467492/100‑1 (MQ=255)
       tACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATggg                                                                                         >  2:2323227/1‑99 (MQ=255)
        aCCCCCGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGGGTTTATCCAGGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTa                                                                                       <  2:798422/100‑1 (MQ=255)
        aCCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTa                                                                                       <  1:1246697/100‑1 (MQ=255)
         cccACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTa                                                                                       <  2:727434/99‑1 (MQ=255)
         cccACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAg                                                                                      >  2:1770042/1‑100 (MQ=255)
         cccACGCTATCTTTAAATTCGGCCTCGACCAGATACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAg                                                                                      >  1:2431034/1‑100 (MQ=255)
             cGCTATCTTTAAATTCGGCCTCGACCAGTCACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGCTGCCCGGTGTGCGTACAGCCGATGGGTAGAATc                                                                                  <  2:1535680/100‑1 (MQ=255)
               cTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGTCCGGTGTGTGTACAGCCGATGGGTAGAATCGa                                                                                >  2:140828/1‑100 (MQ=255)
                tATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGac                                                                               >  1:857269/1‑100 (MQ=255)
                  tCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGATGTGCGTACAGCCGATGGGTAGAATCGacac                                                                             >  2:75246/1‑100 (MQ=255)
                   cTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACAcc                                                                            >  2:630762/1‑100 (MQ=255)
                   cTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACAcc                                                                            <  1:1070756/100‑1 (MQ=255)
                   cTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACAcc                                                                            >  1:2312066/1‑100 (MQ=255)
                   cTTTAAATTCGGCCTCAACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGGGCGTACAGCCGATGGGTAGAATCGACAcc                                                                            >  2:608460/1‑100 (MQ=255)
                    tttAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCAGTGTGCGTACAGCCGATGGGTAGAATCGACACCt                                                                           >  2:1350360/1‑100 (MQ=255)
                    nttAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTANCCACNGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACAcc                                                                            >  1:1223125/2‑99 (MQ=255)
                      tAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGc                                                                         >  2:1502404/1‑100 (MQ=255)
                       aaaTTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCg                                                                        >  1:302051/1‑100 (MQ=255)
                        caTTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGt                                                                       >  1:416569/2‑100 (MQ=255)
                        aaTTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGt                                                                       >  1:1692007/1‑100 (MQ=255)
                         aTTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTg                                                                      >  2:1002334/1‑100 (MQ=255)
                          ttCGGCCTCGATGGGGTACTGCCGGAAAACGGTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTg                                                                      <  2:287466/99‑1 (MQ=255)
                          ttCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTgg                                                                     <  1:315732/100‑1 (MQ=255)
                            cGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGgag                                                                   <  1:624550/100‑1 (MQ=255)
                              gCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGAt                                                                 >  1:903122/1‑100 (MQ=255)
                               ccTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGAtt                                                                >  2:1825181/1‑100 (MQ=255)
                               ccTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGAtt                                                                >  1:2282034/1‑100 (MQ=255)
                                 tCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgc                                                              >  1:2838/1‑100 (MQ=255)
                                  cGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgcc                                                             >  2:108606/1‑100 (MQ=255)
                                   gACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTATAGCCGAAGGGTAGAATCGACACCTGCGTGGAGATTgcca                                                            >  2:1799878/1‑100 (MQ=255)
                                    aCCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgccag                                                           <  1:2061733/100‑1 (MQ=255)
                                     ccAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgccagc                                                          >  1:512271/1‑100 (MQ=255)
                                      cAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgccagcc                                                         >  2:710086/1‑100 (MQ=255)
                                           aCTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCg                                                    >  1:1825859/1‑100 (MQ=255)
                                              gCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATTGGAACCTGCGTGGAGATTGCCAGCCATCCAGaa                                                 >  1:199658/1‑100 (MQ=255)
                                              gCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGaa                                                 >  2:323658/1‑100 (MQ=255)
                                              gCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGa                                                  >  2:328287/1‑99 (MQ=255)
                                               ccGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAg                                                >  2:1495143/1‑100 (MQ=255)
                                                 ggAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTc                                              >  1:880767/1‑100 (MQ=255)
                                                   aaaaCGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAt                                            >  2:1208873/1‑100 (MQ=255)
                                                      aCGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAtctt                                         >  2:1160009/1‑100 (MQ=255)
                                                       cGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAtcttc                                        >  1:1980908/1‑100 (MQ=255)
                                                       cGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAtcttc                                        >  1:1442721/1‑100 (MQ=255)
                                                        gTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCCTCCGGAAGTCAtctt                                         >  2:906286/1‑98 (MQ=255)
                                                         ttGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTg                                      >  1:777833/1‑100 (MQ=255)
                                                           gAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTa                                    <  1:72457/100‑1 (MQ=255)
                                                             gTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTAcc                                  >  1:329264/1‑100 (MQ=255)
                                                             gTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTAcc                                  >  1:1453682/1‑100 (MQ=255)
                                                              tttATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAAGCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCt                                 <  2:230741/100‑1 (MQ=255)
                                                                 aTCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTg                              >  2:2469643/1‑100 (MQ=255)
                                                                 aTCCACGGTCCGGGGTACCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTg                              <  1:1904757/100‑1 (MQ=255)
                                                                  tCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGACTCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCCTCTTCTGTACCttt                               >  2:1093676/1‑98 (MQ=255)
                                                                   ccACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGc                            >  2:956254/1‑100 (MQ=255)
                                                                   ccACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGc                            <  2:186602/100‑1 (MQ=255)
                                                                    cACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACCCCTGCGTGGAGATTGCCAGCCAGCCGGAAGTCATCTTCTGTACCTTTGGCg                           >  2:1768429/1‑100 (MQ=255)
                                                                    cACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCg                           >  1:2149728/1‑100 (MQ=255)
                                                                    cACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCg                           >  2:2535447/1‑100 (MQ=255)
                                                                      cGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGAc                         <  2:1430809/100‑1 (MQ=255)
                                                                       ggTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACg                        >  2:1182018/1‑100 (MQ=255)
                                                                       ggTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACg                        >  2:244496/1‑100 (MQ=255)
                                                                        gTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGc                       >  1:164912/1‑100 (MQ=255)
                                                                        gTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGc                       >  1:1213991/1‑100 (MQ=255)
                                                                         tCCGGGGTGCCCGGTGTGCGGACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGcc                      <  2:256560/100‑1 (MQ=255)
                                                                          ccGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGGGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCa                     >  1:1191961/1‑100 (MQ=255)
                                                                             gggTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATgc                  >  1:950242/1‑100 (MQ=255)
                                                                              ggTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATgcg                 <  1:57428/100‑1 (MQ=255)
                                                                                 gCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGcgcg               <  2:2125301/99‑1 (MQ=255)
                                                                                  cccGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTg             >  2:1616653/1‑100 (MQ=255)
                                                                                  cccGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTg             <  2:2126596/100‑1 (MQ=255)
                                                                                  cccGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTg             <  1:144063/100‑1 (MQ=255)
                                                                                   ccGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGc            <  2:806108/100‑1 (MQ=255)
                                                                                         tgCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACGTTTGGCGACGCCATGGGGGTGCCGGGGa      >  2:467082/1‑100 (MQ=255)
                                                                                         tgCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGa      >  1:306299/1‑100 (MQ=255)
                                                                                         tgCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGa      <  1:257324/100‑1 (MQ=255)
                                                                                          gCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTTTACCTTTGGCGACGCCATGCGCGTGCCGGGGaa     <  2:1670183/100‑1 (MQ=255)
                                                                                          gCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGaa     >  1:1229414/1‑100 (MQ=255)
                                                                                           cGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGaaa    >  1:556349/1‑100 (MQ=255)
                                                                                             tACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCACGTGCCGGGGAAACa  >  1:2525526/1‑100 (MQ=255)
                                                                                                |                                                                                                
GCGGTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACTGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACA  >  NC_000913/2852260‑2852452

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