Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,077,053 G→C 10.7% intergenic (+15/‑59) yeeT → / → cbeA CP4‑44 prophage; DUF987 domain‑containing protein YeeT/CP4‑44 prophage; cytoskeleton bundling‑enhancing antitoxin CbeA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,077,0530GC10.7% 168.9 / 11.3 56intergenic (+15/‑59)yeeT/cbeACP4‑44 prophage; DUF987 domain‑containing protein YeeT/CP4‑44 prophage; cytoskeleton bundling‑enhancing antitoxin CbeA
Reads supporting (aligned to +/- strand):  ref base G (27/23);  new base C (3/3);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

CGCTGGCACGGCAGTGCTGAAGCGTATACCGGTCGTGAAGTGCAGGATATTCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTGGGGGCTGCCCTGCACCGTGACGCCCTGTTTC  >  NC_000913/2076904‑2077201
                                                                                                                                                     |                                                                                                                                                    
cGCTGGCACGGCAGTGCTGAAGCGTATACCGGTCGTGAAGTGCAGGATATTCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCtgatg                                                                                                                                                      <  1:1862123/150‑1 (MQ=255)
   tGGCACGGCAGTGCTGAAGCGTATACCGGTCGTGAAGTGCAGGATATTCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGc                                                                                                                                                   <  1:1077825/150‑1 (MQ=255)
     gCACGGCAGTGCTGAAGCGTATACCGGTCGTGAAGTGCAGGATATTCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAg                                                                                                                                                 <  1:1190912/150‑1 (MQ=255)
            aGTGCTGAAGCGTATACCGGTCGTGAAGTGCAGGATATTCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGCTCATCAGAAttccc                                                                                                                                          <  1:7600/150‑4 (MQ=255)
            aGTGCTGAAGCGTATACCGGTCGTGAAGTGCAGGATATTCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGCTCATCAGAAttccc                                                                                                                                          <  1:866932/150‑4 (MQ=255)
               gCTGAAGCGTATACCGGTCGTGAAGTGCAGGATATTCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGcc                                                                                                                                       >  2:591693/1‑150 (MQ=255)
                        tataCCGGTCGTGAAGTGCAGGATATTCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGa                                                                                                                              <  2:422096/150‑1 (MQ=255)
                         ataCCGGTCGTGAAGTGCAGGATATTCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGAc                                                                                                                             <  1:591693/150‑1 (MQ=255)
                               gTCGTGAAGTGCAGGATATTCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCt                                                                                                                       <  1:1430328/150‑1 (MQ=255)
                               gTCGTGAAGTGCAGGATATTCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCt                                                                                                                       <  1:1688080/150‑1 (MQ=255)
                                tCGTGAAGTGCAGGATATTCCCGGTGTGCTGGCCGTGTTTCCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCtt                                                                                                                      <  1:1260750/150‑1 (MQ=255)
                                 cGTGAAGTGCAGGATATTCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTa                                                                                                                     >  2:1552451/1‑150 (MQ=255)
                                    gAAGTGCAGGATCTTCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGcc                                                                                                                  <  1:271944/150‑1 (MQ=255)
                                    gAAGTGCAGGATATTCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGcc                                                                                                                  >  1:769147/1‑150 (MQ=255)
                                       gTGCAGGATATTCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGAt                                                                                                               <  1:1073476/150‑1 (MQ=255)
                                            ggATATTCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCTGCTGATGAGCGTCACCCTCAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCt                                                                                                          >  2:160717/1‑150 (MQ=255)
                                            ggATATTCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGGTTTTCt                                                                                                          >  1:170405/1‑150 (MQ=255)
                                            ggATATTCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCt                                                                                                          >  2:1811431/1‑150 (MQ=255)
                                             gATATTCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTg                                                                                                         <  1:777235/150‑1 (MQ=255)
                                                 ttCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGAGAGTTTTGGTCCGTATGTCCGGCTGATGAGCTTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAg                                                                                                     <  1:118437/150‑1 (MQ=255)
                                                 ttCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGGTTTTCTGTAAg                                                                                                     <  2:170405/150‑1 (MQ=255)
                                                    ccGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGAt                                                                                                  <  1:1885753/150‑1 (MQ=255)
                                                    ccGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGAt                                                                                                  >  2:1220178/1‑150 (MQ=255)
                                                      ggTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGAttt                                                                                                >  2:1120753/1‑150 (MQ=255)
                                                           gCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCg                                                                                           <  1:565078/150‑1 (MQ=255)
                                                             tGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCgtg                                                                                         >  1:1739487/1‑150 (MQ=255)
                                                             tGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCgtg                                                                                         <  1:128571/150‑1 (MQ=255)
                                                              ggCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCgtgt                                                                                        >  2:1620852/1‑150 (MQ=255)
                                                               gCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTc                                                                                       >  1:1270895/1‑150 (MQ=255)
                                                                ccGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCa                                                                                      >  2:1071747/1‑150 (MQ=255)
                                                                 cGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAg                                                                                     >  2:512952/1‑150 (MQ=255)
                                                                                 cGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGac                                                                     >  1:751707/1‑150 (MQ=255)
                                                                                  gTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGaca                                                                    <  1:160717/150‑1 (MQ=255)
                                                                                     aGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGacaaca                                                                 >  2:1804623/1‑150 (MQ=255)
                                                                                                      tATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCa                                                <  1:1552451/150‑1 (MQ=255)
                                                                                                       aTGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCAc                                               >  2:1211887/1‑150 (MQ=255)
                                                                                                            cGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCg                                          >  2:223440/1‑150 (MQ=255)
                                                                                                            cGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCg                                          <  1:943189/150‑1 (MQ=255)
                                                                                                                     aGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCtggtg                                 <  2:751707/150‑1 (MQ=255)
                                                                                                                      gcctCACCCAGAACTCAGTCGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCtggtgg                                >  2:80457/4‑150 (MQ=255)
                                                                                                                       cGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTggg                               <  2:769147/150‑1 (MQ=255)
                                                                                                                        gTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCGACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTgggg                              >  2:592007/1‑150 (MQ=255)
                                                                                                                          cACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTGGGGGc                            >  2:1952789/1‑150 (MQ=255)
                                                                                                                             ccTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTGGGGGCTGc                         <  2:968014/150‑1 (MQ=255)
                                                                                                                              cTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTGGGGGCTGCC                        >  2:1905572‑M1/1‑150 (MQ=255)
                                                                                                                                gAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTGGGGGCTGCCCT                      >  1:686841‑M1/1‑150 (MQ=255)
                                                                                                                                      aGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTGGGGGCTGCCCTGCACCg                >  1:1078/1‑150 (MQ=255)
                                                                                                                                       gcgataaTTCTGCTCATCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTGGGGGCTGCCCTGCACCGt               >  1:1451127/6‑150 (MQ=255)
                                                                                                                                        cgataattctGCTCATCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTGGGGGCTGCCCTGCACCctg              >  2:852026‑M1/5‑147 (MQ=255)
                                                                                                                                        cgataaTTCTGCTCATCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTGGGGGCTGCCCTGCACCGTg              >  1:1781914/5‑150 (MQ=255)
                                                                                                                                         gggAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTGGGGGCTGCCCTGCACCGTGa             <  1:735764/150‑1 (MQ=255)
                                                                                                                                         gggAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTGGGGGCTGCCCTGCACCGTGa             >  1:328219/1‑150 (MQ=255)
                                                                                                                                              ttCTGCTCATCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTGGGGGCTGCCCTGCACCGTGACGccc        <  2:1451127/150‑1 (MQ=255)
                                                                                                                                               tCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACCCCCTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTGGGGGCTGCCCTGCACCGTGACGCCCt       >  2:1313196/1‑150 (MQ=255)
                                                                                                                                               tCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTgggg                              <  1:112450/127‑1 (MQ=255)
                                                                                                                                               tCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTgggg                              >  2:112450/1‑127 (MQ=255)
                                                                                                                                                cTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTGGGGGCTGCCCTGCACCGTGACGCCCTg      >  1:432547/1‑150 (MQ=255)
                                                                                                                                                 tgatgaGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTGGGGGCTGCCCTGCACCGTGACGCCCTGt     >  2:284321/1‑150 (MQ=255)
                                                                                                                                                  gatgaGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTcccc                                                                          >  1:714882/1‑80 (MQ=255)
                                                                                                                                                  gatgaGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTcccc                                                                          <  2:714882/80‑1 (MQ=255)
                                                                                                                                                    tgaGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTGGGGGCTGCCCTGCACCGTGACGCCCTGTTTc  >  2:258873/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGCTGGCACGGCAGTGCTGAAGCGTATACCGGTCGTGAAGTGCAGGATATTCCCGGTGTGCTGGCCGTGTTTGCTGAACGCCGTAAGGACAGTTTTGGTCCGTATGTCCGGCTGATGAGCGTCACCCTGAACTGAGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTGGGGGCTGCCCTGCACCGTGACGCCCTGTTTC  >  NC_000913/2076904‑2077201

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