Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,403,148 A→C N371H (AAC→CAC)  hflK → modulator for HflB protease specific for phage lambda cII repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,403,1480AC99.4% 614.6 / ‑4.2 175N371H (AAC→CAC) hflKmodulator for HflB protease specific for phage lambda cII repressor
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (100/74);  total (100/75)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 175 total aligned reads displayed.

CGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGG  >  NC_000913/4403052‑4403244
                                                                                                |                                                                                                
cGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACg                                                                                               >  1:392284/1‑100 (MQ=255)
 gCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCTGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACgg                                                                                              <  2:828964/100‑1 (MQ=255)
 gCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACgg                                                                                              >  2:2168740/1‑100 (MQ=255)
 gCAAAGTACTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACTAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAATAGCGATCACgg                                                                                              <  2:3706083/100‑1 (MQ=255)
   aaaGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTg                                                                                            <  1:2608849/100‑1 (MQ=255)
      gTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCa                                                                                         <  1:1953961/100‑1 (MQ=255)
        gCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGc                                                                                       >  1:4061030/1‑100 (MQ=255)
         cTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGGGCGATCACGGTGCCAGCa                                                                                      >  2:3119476/1‑100 (MQ=255)
           ggTTAACGATAAAGGTGGCAAACTGGTGGTTCAGCCGATAGCCCAAATGCTGAAAGGTGGTAACGCCCAGGCGTCGGAGTGCGGTCACGGTGCGAGCAAt                                                                                    >  2:1611508/1‑100 (MQ=255)
             ttAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATct                                                                                  >  1:2964333/1‑100 (MQ=255)
             tcaaCGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATct                                                                                  >  1:2224884/3‑100 (MQ=255)
              tAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATctg                                                                                 >  1:3207859/1‑100 (MQ=255)
              tAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATctg                                                                                 >  2:1086194/1‑100 (MQ=255)
              tAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATctg                                                                                 >  2:1543749/1‑100 (MQ=255)
               aaCGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATctgc                                                                                <  1:2469815/100‑1 (MQ=255)
                 cGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATctgctg                                                                              >  1:1616910/1‑100 (MQ=255)
                 cGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATctgctg                                                                              >  1:924943/1‑100 (MQ=255)
                 cGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATctgctg                                                                              >  2:1306632/1‑100 (MQ=255)
                 cGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATctgctc                                                                              >  2:353814/1‑99 (MQ=255)
                  gATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCtgctgc                                                                             <  2:3182542/100‑1 (MQ=255)
                  gATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGAGGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGCTCACGGTGCCAGCAATCtgctgc                                                                             >  2:3289163/1‑100 (MQ=255)
                    ggagagTGGCAAACTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGACGAGCGATCACGGTGCCCGCAATCTGCTGCGt                                                                           >  1:378650/6‑100 (MQ=255)
                      aaGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGCAGAGCGATCACGGTGCCCGCAATCTGCTGCGTCt                                                                         >  2:2703259/1‑100 (MQ=255)
                      aaGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCt                                                                         >  2:2986631/1‑100 (MQ=255)
                       aGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTg                                                                        <  1:2086112/100‑1 (MQ=255)
                       aGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTg                                                                        >  1:999574/1‑100 (MQ=255)
                         gTGGCAACCTGATGGGTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTgcc                                                                      >  2:399814/1‑100 (MQ=255)
                          tGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTgccg                                                                     <  1:2398811/100‑1 (MQ=255)
                            gCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGGCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTgccgcc                                                                   <  2:3622663/100‑1 (MQ=255)
                              aaCCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAg                                                                 <  1:3799917/100‑1 (MQ=255)
                               aCCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGc                                                                <  2:909006/100‑1 (MQ=255)
                                ccTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGcc                                                               <  2:16784/100‑1 (MQ=255)
                                ccTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGcc                                                               >  1:3992949/1‑100 (MQ=255)
                                ccTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGcc                                                               <  1:1528950/100‑1 (MQ=255)
                                   gATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCctct                                                            >  2:2881614/1‑100 (MQ=255)
                                      ggTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTtcc                                                         <  2:2291044/100‑1 (MQ=255)
                                       gTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTtcct                                                        <  2:1282871/100‑1 (MQ=255)
                                        ttCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTtcctc                                                       <  1:473961/100‑1 (MQ=255)
                                         tCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTtcctcc                                                      >  1:2060399/1‑100 (MQ=255)
                                         tCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTtcctcc                                                      <  1:288533/100‑1 (MQ=255)
                                          cTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGTCAGCCTCTTCCTCca                                                     >  2:307407/1‑100 (MQ=255)
                                           tGCCGTTGGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGTAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCcac                                                    <  1:3750909/100‑1 (MQ=255)
                                              cGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAAc                                                 >  1:821909/1‑100 (MQ=255)
                                              cGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAAc                                                 >  1:1026916/1‑100 (MQ=255)
                                              cGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCACCAGCCTCTTCCTCCACAAc                                                 >  1:1401930/1‑100 (MQ=255)
                                                 tAGACCAGATGCGGAAAGGTGGTAACGCCCCGGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAg                                              <  1:31758/100‑1 (MQ=255)
                                                  aGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGt                                             >  2:855614/1‑100 (MQ=255)
                                                   gACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTg                                            >  1:2694097/1‑100 (MQ=255)
                                                   gACCAGATGCAGAGAGGTGGTAACGCTACTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTAACGCCAGCCTCTTCCTCCACAACCAGTg                                            >  2:1344052/1‑100 (MQ=255)
                                                    aCCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTgg                                           >  2:3480872/1‑100 (MQ=255)
                                                     ccAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGa                                          >  1:537042/1‑100 (MQ=255)
                                                     ccAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGa                                          >  2:2690861/1‑100 (MQ=255)
                                                     ccAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGa                                          >  2:1106597/1‑100 (MQ=255)
                                                      cAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGag                                         >  1:962395/1‑100 (MQ=255)
                                                      cAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGag                                         <  1:1749305/100‑1 (MQ=255)
                                                      cAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGag                                         <  2:1424773/100‑1 (MQ=255)
                                                      cAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGag                                         <  2:2146991/100‑1 (MQ=255)
                                                       agagGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGGAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagc                                        <  2:482950/100‑1 (MQ=255)
                                                       aGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagc                                        >  2:1303932/1‑100 (MQ=255)
                                                         aTGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagcaa                                      >  2:2487471/1‑100 (MQ=255)
                                                          tGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagcaag                                     >  2:2228921/1‑100 (MQ=255)
                                                           gCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagcaagc                                    >  1:3281057/1‑100 (MQ=255)
                                                           gCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagcaagc                                    >  1:3436885/1‑100 (MQ=255)
                                                           gCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagcaagc                                    >  2:2609637/1‑100 (MQ=255)
                                                           gCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagcaagc                                    >  1:1900393/1‑100 (MQ=255)
                                                            cTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagcaagca                                   >  2:1856793/1‑100 (MQ=255)
                                                            cTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagcaagca                                   <  2:57242/100‑1 (MQ=255)
                                                            cTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagcaagca                                   <  1:2750420/100‑1 (MQ=255)
                                                             tGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAgcaagcaa                                  >  2:3221342/1‑100 (MQ=255)
                                                              gAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAac                                 >  1:849685/1‑100 (MQ=255)
                                                              gAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAac                                 <  1:1105743/100‑1 (MQ=255)
                                                              gAAACGTGGTAACGCCCCTTCGGCGATGAGCGATCACGGTGCCAGCAATCTTCTGCGTCTGCCGCCAGGCTCTTCCTCCACAACCAGTGGAGCAAGCAac                                 <  2:385158/100‑1 (MQ=255)
                                                                aaGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAacac                               >  2:1952591/1‑100 (MQ=255)
                                                                 aGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAacacn                              <  1:1426136/100‑2 (MQ=255)
                                                                 aGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACAcg                              >  1:2061969/1‑100 (MQ=255)
                                                                  ggtggtAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACAcgt                             <  1:312892/100‑1 (MQ=255)
                                                                   gtggtAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCGGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACAcgtc                            <  2:3151474/100‑1 (MQ=255)
                                                                    tggtAACGCCCCGGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACAcgtcg                           <  2:2268853/100‑1 (MQ=255)
                                                                     ggtAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACAcgtcgt                          >  1:286816/1‑100 (MQ=255)
                                                                     ggtAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACAcgtcgt                          >  1:1271851/1‑100 (MQ=255)
                                                                       tAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTcc                        >  2:493431/1‑100 (MQ=255)
                                                                          ngCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTccacc                     >  1:3260957/2‑100 (MQ=255)
                                                                          cGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTccacc                     >  1:3498160/1‑100 (MQ=255)
                                                                          cGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTccacc                     >  2:342743/1‑100 (MQ=255)
                                                                            ccccTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAg                   >  2:2205715/1‑100 (MQ=255)
                                                                              ccTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTc                 >  2:935547/1‑100 (MQ=255)
                                                                              ccTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTc                 <  1:3714847/100‑1 (MQ=255)
                                                                              ccTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTc                 >  1:201936/1‑100 (MQ=255)
                                                                               cTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCa                >  2:194787/1‑100 (MQ=255)
                                                                                tGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAg               >  2:618984/1‑100 (MQ=255)
                                                                                 gcggcgAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAgg              <  2:2863128/100‑1 (MQ=255)
                                                                                  cggcgAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAggg             <  1:909731/100‑1 (MQ=255)
                                                                                  cggcgAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAggg             >  1:2529299/1‑100 (MQ=255)
                                                                                    gcgAAGAGCGCTCTCGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCTAGCAACACGTCGTCCACCAGTCAGGGCg           <  1:1030385/100‑1 (MQ=255)
                                                                                    gcgAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCg           <  1:2201188/100‑1 (MQ=255)
                                                                                     cgAAGGGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGa          <  2:1070559/100‑1 (MQ=255)
                                                                                           gCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCGCAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGAtattat    >  1:393698/1‑100 (MQ=255)
                                                                                           gCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGAtattat    >  1:3003362/1‑100 (MQ=255)
                                                                                            cGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATCATg   <  1:3138283/100‑1 (MQ=255)
                                                                                             gATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATgg  >  1:2108520/1‑100 (MQ=255)
                                                                                                |                                                                                                
CGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGG  >  NC_000913/4403052‑4403244

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: