Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,852,356 T→A 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.9% 311.2 / ‑4.0 88L74Q (CTG→CAG) hypDhydrogenase maturation protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (46/41);  total (47/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CGGTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACTGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACA  >  NC_000913/2852261‑2852452
                                                                                               |                                                                                                
cGGTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCg                                                                                              <  2:1900038/100‑1 (MQ=255)
    naTACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATgg                                                                                          >  2:796659/2‑100 (MQ=255)
    cctACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATgg                                                                                          >  1:1374247/3‑100 (MQ=255)
    cATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATgg                                                                                          >  1:1199785/1‑100 (MQ=255)
    cATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATgg                                                                                          <  1:1849064/100‑1 (MQ=255)
     aTACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATggg                                                                                         <  2:980734/100‑1 (MQ=255)
     aTACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATggg                                                                                         <  2:22914/100‑1 (MQ=255)
     aTACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATggg                                                                                         <  2:1879247/100‑1 (MQ=255)
      tACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGt                                                                                        <  2:1805428/100‑1 (MQ=255)
      tACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGt                                                                                        >  1:1519978/1‑100 (MQ=255)
       aCCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTa                                                                                       >  1:1056865/1‑100 (MQ=255)
        cccACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAg                                                                                      <  2:848808/100‑1 (MQ=255)
        cccACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAg                                                                                      >  1:803550/1‑100 (MQ=255)
        cccACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAg                                                                                      <  2:1678206/100‑1 (MQ=255)
          cACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGaa                                                                                    >  2:1620850/1‑100 (MQ=255)
           aCGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAAt                                                                                   <  1:78895/100‑1 (MQ=255)
            cGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATc                                                                                  >  1:754068/1‑100 (MQ=255)
            cGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATc                                                                                  <  1:337212/100‑1 (MQ=255)
             gCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACTGCCGATGGGTAGAATCg                                                                                 >  2:1046202/1‑100 (MQ=255)
             gCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCg                                                                                 <  2:1173707/100‑1 (MQ=255)
             gCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCg                                                                                 >  2:725861/1‑100 (MQ=255)
                 tCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGacac                                                                             >  1:32263/1‑100 (MQ=255)
                  cTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCTCGGTGTGCGTACAGCCGATGGGTAGAATCGACAcc                                                                            >  1:1696306/1‑100 (MQ=255)
                   tttAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCt                                                                           >  2:1717495/1‑100 (MQ=255)
                    ttAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCtn                                                                          <  2:335841/100‑2 (MQ=255)
                      aaaTTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCg                                                                        <  1:297759/100‑1 (MQ=255)
                          tCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGga                                                                    <  2:1374238/100‑1 (MQ=255)
                           cGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGgag                                                                   <  2:1848021/100‑1 (MQ=255)
                            ggCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGgaga                                                                  <  2:1199779/100‑1 (MQ=255)
                               cTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACAGCTGCGTGGAGATTg                                                               >  2:1036034/1‑100 (MQ=255)
                               cTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTg                                                               <  1:1915599/100‑1 (MQ=255)
                               cTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAAAATCGACACCTGCGTGGAGATTg                                                               >  1:207852/1‑100 (MQ=255)
                                 ctcttccgatcTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgcc                                                             <  1:1088530/90‑1 (MQ=255)
                                  gACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGATATTgcca                                                            <  1:1077663/100‑1 (MQ=255)
                                   aCCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgccag                                                           >  2:1830170/1‑100 (MQ=255)
                                   aCCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgccag                                                           <  2:944135/100‑1 (MQ=255)
                                    ccAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgccagc                                                          <  2:1056860/100‑1 (MQ=255)
                                       gTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCAt                                                       <  2:259220/100‑1 (MQ=255)
                                          aCTGCCGGAAAACGTTGAGATTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCg                                                    <  1:1690299/100‑1 (MQ=255)
                                             gCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACGCCTGCGTGGAGATTGCCAGCCATCCGGaa                                                 >  2:1226797/1‑100 (MQ=255)
                                             gCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgccagatcggaagag                                                 >  2:1088525/1‑90 (MQ=255)
                                             gCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGaa                                                 >  1:1564189/1‑100 (MQ=255)
                                              ccGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAg                                                <  2:1564179/100‑1 (MQ=255)
                                              ccGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAg                                                <  1:1039240/100‑1 (MQ=255)
                                                ggAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTc                                              >  1:915744/1‑100 (MQ=255)
                                                ggAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTc                                              >  2:1252067/1‑100 (MQ=255)
                                                ggAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTc                                              >  1:1290191/1‑100 (MQ=255)
                                                ggAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTc                                              >  1:894800/1‑100 (MQ=255)
                                                  aaaaCGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAt                                            >  1:37576/1‑100 (MQ=255)
                                                  aaaaCGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAt                                            <  2:753780/100‑1 (MQ=255)
                                                  aaaaCGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAt                                            >  2:300376/1‑100 (MQ=255)
                                                    aaCGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAtct                                          >  1:507232/1‑100 (MQ=255)
                                                      cGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAtcttc                                        <  2:536004/100‑1 (MQ=255)
                                                      cGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAtcttc                                        <  1:1442549/100‑1 (MQ=255)
                                                       gTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAtcttct                                       >  2:173733/1‑100 (MQ=255)
                                                       gTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAtcttct                                       >  2:162374/1‑100 (MQ=255)
                                                           aGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTAc                                   >  2:274471/1‑100 (MQ=255)
                                                            gTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTAcc                                  >  1:91295/1‑100 (MQ=255)
                                                              ttATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGACGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCt                                 <  1:1433264/99‑1 (MQ=255)
                                                                 tCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTgg                             >  2:1655767/1‑100 (MQ=255)
                                                                  ccACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTgg                             >  1:1796528/1‑99 (MQ=255)
                                                                  ccACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGc                            >  1:1288099/1‑100 (MQ=255)
                                                                  ccACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGc                            <  2:1473420/100‑1 (MQ=255)
                                                                  ccACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGc                            >  2:1018015/1‑100 (MQ=255)
                                                                  ccACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGc                            >  2:1820587/1‑100 (MQ=255)
                                                                  ccACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGc                            <  2:1775494/100‑1 (MQ=255)
                                                                  ccACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGc                            >  2:1610324/1‑100 (MQ=255)
                                                                  ccACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGc                            >  1:409370/1‑100 (MQ=255)
                                                                   cACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCg                           >  2:1609758/1‑100 (MQ=255)
                                                                   cACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCg                           >  2:198762/1‑100 (MQ=255)
                                                                   cACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCg                           <  2:1273682/100‑1 (MQ=255)
                                                                   cACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCg                           >  2:1473407/1‑100 (MQ=255)
                                                                   cACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCg                           >  2:13292/1‑100 (MQ=255)
                                                                   cACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCg                           <  1:1473416/100‑1 (MQ=255)
                                                                   cACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGCAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCGGTACCTTTGGc                            >  2:1974438/1‑99 (MQ=255)
                                                                      ggTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACg                        <  2:803547/100‑1 (MQ=255)
                                                                      ggTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACg                        <  1:1365075/100‑1 (MQ=255)
                                                                        tCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGcc                      <  2:1519969/100‑1 (MQ=255)
                                                                           ggggTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATg                   >  1:1157143/1‑100 (MQ=255)
                                                                            gggTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATgc                  <  2:91294/100‑1 (MQ=255)
                                                                               tGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGcgcg               >  1:1364536/1‑100 (MQ=255)
                                                                                gCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACAACTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGt              >  1:1443406/1‑100 (MQ=255)
                                                                                   cGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGcc           <  2:1696293/100‑1 (MQ=255)
                                                                                     gtgtGCGGACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCgg         <  1:1310220/100‑1 (MQ=255)
                                                                                      tgtgCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCggg        <  1:173734/100‑1 (MQ=255)
                                                                                       gtgCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCgggg       <  2:32263/100‑1 (MQ=255)
                                                                                        tgCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGa      >  1:679168/1‑100 (MQ=255)
                                                                                            tACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACa  >  2:337036/1‑100 (MQ=255)
                                                                                               |                                                                                                
CGGTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACTGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACA  >  NC_000913/2852261‑2852452

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