Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 550,098 T→C 6.1% W158R (TGG→CGG)  ylbF → putative anaerobic allantoin catabolic oxamate carbamoyltransferase; DUF2877 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913550,0980TC 6.1% 199.1 / 3.2 65W158R (TGG→CGG) ylbFputative anaerobic allantoin catabolic oxamate carbamoyltransferase; DUF2877 family protein
Reads supporting (aligned to +/- strand):  ref base T (31/30);  new base C (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

CAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACA  >  NC_000913/549952‑550244
                                                                                                                                                  |                                                                                                                                                  
cAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCActg                                                                                                                                                   >  8:77268/1‑148 (MQ=255)
cAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCActg                                                                                                                                                   <  7:77268/148‑1 (MQ=255)
cAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCActg                                                                                                                                                   >  4:212232/1‑148 (MQ=255)
cAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTGGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCActg                                                                                                                                                   <  3:212232/148‑1 (MQ=255)
      cATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCActggctgg                                                                                                                                              <  7:139/147‑1 (MQ=255)
        tGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTggg                                                                                                                                             <  5:2185/146‑1 (MQ=255)
        tGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTggg                                                                                                                                             >  6:2185/1‑146 (MQ=255)
          gATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGGTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTaa                                                                                                                                          <  5:69019/147‑1 (MQ=255)
             gCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGc                                                                                                                                      >  3:67762/1‑148 (MQ=255)
               aaCGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGccc                                                                                                                                    >  4:341768/1‑148 (MQ=255)
                    gcgGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGAt                                                                                                                               <  8:70507/148‑1 (MQ=255)
                      ggaggaAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTa                                                                                                                             <  3:341543/148‑1 (MQ=255)
                        aggaAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCTCTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAAc                                                                                                                           >  4:134989/1‑148 (MQ=255)
                            aaCCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGcc                                                                                                                        >  6:218781/1‑147 (MQ=255)
                                  gCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCat                                                                                                                 <  1:361018/148‑1 (MQ=255)
                                   cTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGCCCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCat                                                                                                                 >  5:52305/1‑147 (MQ=255)
                                   cTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCatg                                                                                                                <  3:277083/148‑1 (MQ=255)
                                   attttCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCatg                                                                                                                <  4:67762/147‑1 (MQ=255)
                                     tttCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTTGGTAAAGGCCCCGGATTAACGCCGAGTCattat                                                                                                              >  2:126798/1‑148 (MQ=255)
                                      ttCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGcc                                                                                                                                     <  8:84779/124‑1 (MQ=255)
                                      ttCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGcc                                                                                                                                     >  7:84779/1‑124 (MQ=255)
                                      ttCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCatga                                                                                                               <  8:234813/146‑1 (MQ=255)
                                      ttCGGGCCACTGGCGGTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCatgat                                                                                                              <  5:218781/147‑1 (MQ=255)
                                        cGGTCCACTGGCGGTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGCTTTCCGTCCGCGCTCAATGGCGTTAAGACCGTCTGGCGGCACTGGCTGGGTAACGGCCCCGGATTAACGCCGAGTCAtgatga                                                                                                            <  5:125702/147‑1 (MQ=255)
                                        cGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCAtgatga                                                                                                            <  5:125701/147‑1 (MQ=255)
                                         gggCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGaca                                                                                                          >  3:11793/1‑148 (MQ=255)
                                           gCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCAAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACg                                                                                                        >  5:194370/1‑148 (MQ=255)
                                                 ggCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGc                                                                                                  >  1:110762/1‑148 (MQ=255)
                                                 ggCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGc                                                                                                  >  1:356107/1‑148 (MQ=255)
                                                      tGGCTGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGATTCGTCACTGTTTACAGGCCGCGATCAATGGCGTTAATACCAACTGGCGGCACTGGCTGGGTAAAGTCCCCTGATTACCGCCGAGTAAAGATGACACGCTGACCGGAAt                                                                                             >  6:102883/1‑148 (MQ=255)
                                                       ggcggcAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGa                                                                                               >  6:207502/1‑145 (MQ=255)
                                                       ggcggcAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGa                                                                                               <  5:207502/145‑1 (MQ=255)
                                                              aGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTgg                                                                                     <  3:75611/148‑1 (MQ=255)
                                                                     cGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCt                                                                               <  4:177722/147‑1 (MQ=255)
                                                                       ctgcctgcTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTTGGTAAAGGCCCCGTATTAACTCGGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGt                                                                            >  5:103078/1‑148 (MQ=255)
                                                                           ctgctgAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGtatta                                                                        <  5:291831/148‑1 (MQ=255)
                                                                             gctgAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACCGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCt                                                                               >  3:284806/1‑139 (MQ=255)
                                                                             gctgAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACCGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCt                                                                               <  4:284806/139‑1 (MQ=255)
                                                                                gAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGgcg                                                                   >  2:31922/1‑148 (MQ=255)
                                                                                   ttACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGAACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCtt                                                                >  5:246858/1‑148 (MQ=255)
                                                                                        ccAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGAt                                                           >  2:325689/1‑148 (MQ=255)
                                                                                          aGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATgc                                                         >  8:1547/1‑148 (MQ=255)
                                                                                             ttCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGtatt                                                                         >  8:190855/1‑129 (MQ=255)
                                                                                             ttCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGtatt                                                                         <  7:190855/129‑1 (MQ=255)
                                                                                                gTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCt                                                    <  5:141034/147‑1 (MQ=255)
                                                                                                gTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCt                                                    <  1:325689/147‑1 (MQ=255)
                                                                                                   aCTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGATGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCgg                                                <  1:126798/148‑1 (MQ=255)
                                                                                                    cTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCAGGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGgt                                               >  6:241177/1‑148 (MQ=255)
                                                                                                      gTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGgtcg                                             >  8:55862/1‑148 (MQ=255)
                                                                                                       ttttCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGgtcgt                                            >  1:100967/1‑148 (MQ=255)
                                                                                                        tttCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGATCCGgtcgtc                                           >  2:34495/1‑148 (MQ=255)
                                                                                                              gCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGtttt                                     <  2:100967/148‑1 (MQ=255)
                                                                                                               ccGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGtatt                                                                         <  1:37552/111‑1 (MQ=255)
                                                                                                               ccGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGtatt                                                                         >  2:37552/1‑111 (MQ=255)
                                                                                                                 gcgcTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTg                                  <  7:55862/148‑1 (MQ=255)
                                                                                                                   gcTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGcc                                >  3:134437/1‑148 (MQ=255)
                                                                                                                      cAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGt                             <  3:134989/148‑1 (MQ=255)
                                                                                                                           gCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCg                         <  3:263851/147‑1 (MQ=255)
                                                                                                                              ttAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGt                                                                            >  6:285128/1‑93 (MQ=255)
                                                                                                                              ttAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGt                                                                            <  5:285128/93‑1 (MQ=255)
                                                                                                                               tAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAAt                    >  5:357751/1‑148 (MQ=255)
                                                                                                                               tAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAAt                    >  8:180350/1‑148 (MQ=255)
                                                                                                                                aaGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAAt                    <  5:92045/147‑1 (MQ=255)
                                                                                                                                 aGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCt                  <  1:31922/148‑1 (MQ=255)
                                                                                                                                      gACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGGCCGTTTTTTGCCTGTTCCGACAATCTGCAAc             >  8:311962/1‑148 (MQ=255)
                                                                                                                                       aCTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCGGGGAGTATGGCGCGTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACt            >  2:227124/1‑148 (MQ=255)
                                                                                                                                           gcggcACTGGCCGGGGAAACGCCCCGGATTAACGCCGACTCCTGATGCCACGCCGACCCGAATGGGTCAGGCCGCCCGGTATTATGGCCCCTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGtt        <  4:11793/148‑1 (MQ=255)
                                                                                                                                             ggcACCGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCa               <  1:237430/139‑1 (MQ=255)
                                                                                                                                             ggcACCGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCa               >  2:237430/1‑139 (MQ=255)
                                                                                                                                                aCTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACcac   <  8:218735/148‑1 (MQ=255)
                                                                                                                                                 ctggctggGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACcaca  <  3:312651/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACA  >  NC_000913/549952‑550244

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